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INTRODUCTION 

In  many  problems  of  practical  importance  one  would  like  to 
identify  unknown  parameters  in  mathematical  models  given  certain 
observations  of  the  underlying  physical  phenomenon  being  modeled. 

A  general  framework  for  approximating  or  estimating  unknown  para¬ 
meters  in  partial  differential  equations,  using  modal  (eigenfunc¬ 
tion)  approximations,  was  presented  in  [11].  The  general  theo¬ 
retical  framework  developed  in  [11]  was  subsequently  applied  to 
spline-based  techniques  in  [  7 ]  to  two  classes  of  second  order 
initial -boundary  value  problems. 

This  thesis  is  devoted  to  developing  numerical  methods  for 
estimating  unknown  constant  parameters  in  certain  fourth  order 
constant  coefficient  partial  differential  equations.  The  approxi¬ 
mation  techniques  follow  the  approach  in  [  7 ]  and  convergence  is 
proved  using  the  theoretical  convergence  framework  developed  in 
[11]  employing  linear  semigroup  theory.  The  necessary  theoreti¬ 
cal  framework  is  summarized  in  Chapter  1. 

We  treat  identification  problems  for  two  specific  equations 
in  one  dimension  which  model  the  transverse  vibration  of  an  elas¬ 
tic  or  viscoelastic  beam,  and  develop  numerical  methods  for  esti¬ 
mating  unknown  parameters. 

We  examine  estimation  techniques  for  the  Euler-Bernoul li 
equation  in  Chapter  2  and  for  the  Timoshenko  equation  in  Chapter 
3.  In  both  cases  we  develop  methods  for  numerically  estimating 
unknown  parameters  and  prove  convergence  of  the  methods.  Numeri¬ 
cal  results  are  provided  to  illustrate  the  theoretical  convergence 
results . 


The  Euler-Bernoulli  equation  which  we  examine  includes  both 
structural  (internal)  and  viscous  damping,  and  various  boundary 
conditions  are  used.  We  introduce  two  methods  for  estimating 
parameters,  one  based  on  quintic  splines  (Section  2)  and  one  based 
on  cubic  spline  approximations  (Section  3) .  For  the  Timoshenko 
equation  (Chapter  3)  we  examine  one  technique  for  estimation  of 
parameters  based  upon  cubic  spline  approximations. 

In  Chapter  4  we  discuss  the  implementation  of  the  approxima¬ 
tion  techniques  into  a  computer  code.  Since  all  of  the  computer 
codes  used  here  for  the  fourth  order  problems  and  in  [  7 ]  for 
certain  second  order  problems  have  the  same  general  structure,  th 
discussion  of  the  implementation  is  made  sufficiently  general  to 
describe  all  of  the  spline-based  methods  which  we  have  developed. 

The  notation  employed  throughout  this  thesis  is  rather  stan¬ 
dard.  For  norms  of  elements  in  Banach  spaces  we  use  | • | ,  whereas 
| | *||  is  used  for  operator  norms.  A  subscripted  norm  | • |  de¬ 
notes  certain  norms  equivalent  to  the  usual  norms  on  the  Sobolev 

spaces  Hm  over  [0,1],  and  specifically  I’Iq  denotes  the  H^1 
2 

(L  )  norm.  Similarly,  inner  products  on  certain  subspaces  of 

the  Sobolev  spaces  Hm  will  be  denoted  by  <•»•>  .  These  will  be 

m 

defined  in  Section  2  of  Chapter  1.  As  we  shall  be  dealing  with 
state  spaces  Z  which  are  products  of  function  spaces,  the  sym¬ 
bol  |*|  with  no  subscript  (and  similarly  <•,•>)  will  be  reserved 
to  denote  the  norm  (or  inner  product)  on  the  state  space  Z. 


CHAPTER  1.  FOUNDATIONS 


Section  1 .  The  Identification  Problem  and  Its  Approximation 

We  begin  by  defining  the  identification  or  estimation  problem 
for  a  process  governed  by  a  partial  differential  equation  and  pro¬ 
ceed  to  its  abstract  formulation.  A  general  framework  for  ap¬ 
proximating  solutions  is  introduced  along  with  tools  necessary  to 
establish  convergence.  The  techniques  for  approximating  solutions 
to  the  identification  problem  follow  those  introduced  in  [11], 
where  modal  (eigenfunction)  state  approximations  were  applied  to 
a  class  of  hyperbolic  and  parabolic  equations,  and  also  used  in 
[  7  ] ,  where  spline-based  state  approximations  were  applied  to  the 
same  class  of  problems. 

This  section  outlines  the  basic  approach  and  theory  which  will 
be  applied  in  later  chapters  to  two  specific  equations,  namely 
the  Euler-Bernoulli  and  Timoshenko  equations  for  the  transverse 
vibration  of  a  beam,  and  to  specific  spline-based  approximations 
for  the  identification  problem. 

We  first  define  the  identification  problem  for  an  initial¬ 
boundary  value  problem.  Suppose  we  have  a  physical  process 
modeled  by  an  initial-boundary  value  problem  with  unknown  para¬ 
meters  q  =  (q^,...,q  )  €  Rp;  the  parameters  q^  may  be  unknown 
constant  coefficients  in  the  partial  differential  equation  or 
parameters  appearing  in  functions  in  the  initial  conditions  or 
non-homogeneous  term.  Suppose  also  that  we  are  given  a  set  of 
output  measurements  from  the  physical  process  which  is  modeled  by 
the  initial -boundary  value  problem.  In  a  sense  which  will  be 
made  more  precise  in  what  follows,  the  identification  (or  estima- 
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tion)  problem  consists  of  finding  the  vector  of  parameters  q  €  R^ 
such  that  the  solution  of  the  initial -boundary  value  problem 
"best  fits"  the  output  measurements  of  the  physical  process. 

he  will  be  interested  in  two  particular  init ial -boundary  value 
problems,  namely  the  Euler-Bernoulli  and  Timoshenko  equations  for 
the  transverse  vibration  of  a  uniform  beam.  The  Euler-Bernoulli 
equation,  including  structural  and  viscous  damping,  has  the  form 

C 1  - 1 )  ytt  *  -q1D4y  -  q2°4yt  "  q3y  +  f(t,x;q),  t  >  0,  x  €  [0,1] 

y(0,x)  =  y  q ( x ; q ) 
yt(0,x)  =  y:(x;q) , 

with  appropriate  homogeneous  boundary  conditions  at  x  =  0  and 
x  =  1,  where  5  8^/8xJ,  y(t,x;q)  is  the  transverse  displace¬ 

ment,  and  f(t,x;q)  is  the  applied  load.  Here,  qj  ,  q2>  and  q^ 
are  unknown  constant  coefficients,  and  q4,...,qp  are  parameters 
appearing  in  the  nonhomogeneous  (load)  term  and  initial  conditions. 
The  Timoshenko  equations  can  be  put  in  the  form 

(1.2)  ytt  =  qxD2y  -  q1Dijj  +  f(t,x;q) 

2  t  >  0,  x  €  [0,1] 

♦tt  =  q5D  *  +  q2(Dy  ' 

y(o.x)  =  y0(x;q) 

yt(0,x)  =  y1(x;q) 

V(0,x)  =  i|>0(x;q) 

*t(0,x)  =  4»xCx;q) 

where  again  we  associate  appropriate  homogeneous  boundary  condi¬ 
tions  with  (1.2)  at  x  =  0  and  x  =  1.  Here  y(t,x;q)  is  trans¬ 
verse  displacent  and  i[/(t,x;q)  is  the  angle  of  rotation  of  a 
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cross-section  of  the  beam. 

More  will  be  said  about  these  equations  in  the  sequel.  The 
above  is  sufficient  at  present  to  set  the  framework  for  the 
identif ication  problem. 

We  shall  be  interested  in  the  ident ii ication  problem  associa- 
ted  with  (1.1)  or  (1.2).  Given  a  set  of  observations  n  = 
where  =  (y  (t  i  ,x1 )  , .  .  .  ,y  (t  i  ,  x^)  )  T  and  yCt^x^)  is  the  ob¬ 

served  displacement  at  t^,  x^  ,  of  a  process  which  we  assume  to 
be  modeled  by  (1.1)  or  (1.2),  find  the  vector  of  parameters 
q  =  (q^,...,q  )  which  minimizes  J(q;y,n),  where  J  is  some 
functional  measuring  the  fit  of  (1.1)  or  (1.2)  to  the  data  n^. 
Specifically,  we  shall  use  a  pointwise  fit-to-data  criterion  of 
the  form 


(1.3)  J (q,y ,n)  =  l  U-  -  h-  I2 

n=  1  1  1 

T 

where  ^  =  (y(t . ,x:) , . . . ,y ( t ± ,x£) )  ,  and  (t,x)  » — =►  y(t,x) 
(respectively,  (t,x)  i — >  (y  (t  ,x)  ,<P(t  ,x) )  is  the  solution  of 
(1.1)  (respectively,  of  (1.2)). 

To  ensure  that  the  initial -boundary  value  problems  are  well- 
posed,  we  shall  assume  hereafter 

1)  Q  is  a  compact  set  in  rP  ,  and 
(HQ)  2)  there  exists  q^  >  0  such  that  q^  >  q^  for  all 

q  €  Q. 

For  the  Timoshenko  equations,  we  also  require,  along  with  the 
above 
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2')  there  exists  >  0  such  that  >_  q^  for  all  q  €  Q. 

Remark .  Note  that  we  have  formulated  the  initial-boundary 
value  problem  for  homogeneous  boundary  conditions.  If  nonhomo- 
geneous  time  dependent  boundary  conditions  occur,  of  the  form 
(DJy) |x=Q  =  h(t,q)  they  will  be  transformed  to  homogeneous  bound¬ 
ary  conditions  by  a  transformation  of  the  form  v(t,x)  =  y(t,x)  - 
g(x)h(t,q),  (see  [29],  where  a  general  method  for  transforming 

(1.1)  with  time-dependent  boundary  conditions  to  (1.1)  with  homo¬ 
geneous  boundary  conditions,  in  the  case  q2  =  0,  q^  =  0,  is  dis¬ 
cussed)  .  Khen  such  a  transformation  is  performed  and  the  output 
measurements  in  the  fit-to-data  criterion  correspond,  as  above, 
to  y,  it  is  appropriate  to  modify  J  in  (1.3)  to 

d.3')  J’(q,v,n)  =  l  |C!  -  n .  1 2 , 

i  =  l  1  1 

where 

Ci  =  (v(ti»xi^  +  g(x1)h(ti,q) , . . . ,v(ti,x£)  +  g(x£)h(ti;q))T. 

Definition  1.1.  The  identification  problem  (ID)  for  (1.1)  (or 

(1.2) )  is  defined  as  the  following:  given  (1.1)  with  unknown 
parameters  q  =  (q^,...,q  )  6  Rp  and  a  family  of  solutions 
(t,x)  > — ^  y(t,x;q)  ((y(t,x;q) ,^(t,x;q))T  for  (1.2)),  and  a  set 
of  output  measurements  (n(ti)}^_1,  find  q  €  Q,  where  Q  is  some 
parameter  set  in  rP  satisfying  (HQ)  ,  such  that 

J(q,y(* » • ;q) ,n)  <  J(q,y(* , • ;q) ,n)  for  ail  q  e  Q, 

where  J  is  the  cost  functional  (1.3). 

Having  formulated  the  identification  problem  for  two  model 
initial-boundary  problems,  we  proceed  to  an  abstract  formulation 
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of  the  problem.  First  we  write  the  initial -boundary  value  problem 
(1.1)  or  (1.2)  as  an  abstract  equation  on  an  appropriate  Hilbert 
space  Z,  in  the  form 

(1.4)  z(t)  =  _a^(q)z(t)  +  F(q,t)  for  t  >  0 
z(0)  =  zQ 

where  z(t)  =  (z^ (t) . zn(t) )  £  Z  for  every  t  £  [0,T],  and 

z-^(t)  corresponds  to  y(t,*)  in  (1.1)  or  (1.2). 

Specific  formulations  of  (1.1)  and  (1.2)  as  an  abstract  equa¬ 
tion  will  be  examined  in  Chapters  2  and  3.  In  each  case,  the  ab¬ 
stract  equation  will  be  formulated  in  such  a  way  that  _o^(q)  gen¬ 
erates  a  CQ  semigroup  T(t;q)  on  Z;  i.e.,  t  • — >  T(t,q)z0 
is  the  solution  to  z(t)  =  _Q^(q)z(t),  z(0)  =  Zq. 

The  notion  of  dissipativeness  will  play  a  central  role  in 
proving  that  ja/(q)  is  generator  of  a  CQ  semigroup  in  the  speci¬ 
fic  examples  to  be  considered.  A  densely  defined  operator  stf  is 
called  dissipative  if  <o£,z>  <  0  for  all  z  €  Vomfaf)  ,  and  maxi¬ 
mal  dissipative  if  its  only  dissipative  extension  is  itself. 

While  there  are  many  simple  conditions  which  guarantee  that  a 
dissipative  operator  generates  a  Cq  semigroup,  we  shall  use  one 
primarily:  a  maximal  dissipative  operator  szf  generates  a  Cq 
semigroup  (T(t)}  of  contraction  on  Z;  i.e.,  || T  ||  <.  1  (cf :  [25, 

p.  88]). 

We  will  be  interested  in  mild  solutions  of  (1.4):  t  i — >  z(t;q) 
is  called  a  mild  solution  on  [0,T]  if  it  satisfies 

(1.5)  z(t;q)  =  T (t ; q) Zq  +  f  T (t -s ; q) F (q , s) ds 

u  Jo 

for  every  t  €  [0,T] .  We  will  also  place  the  following  conditions 
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on  F  which  guarantee  uniqueness  and  existence  of  mild  solutions 
z(-,q)  £  C(0,T;Z)  to  (1.4)  (see  [11,  p.  13]): 

1)  The  map  t  ( — >  F(q,t)  is  measurable, 

(HF)  2)  the  map  q  i — ^  F(q,t)  is  continuous,  and 

3)  there  exists  k(t)  €  L2(0,T)  such  that 
I F (q , t) |  <  k(t). 

We  note  in  passing  that  when  strong  (classical)  solutions  to 
(1.4)  exist  they  will  be  mild  solutions.  In  fact,  if  zQ  € 
Dom(a/(q))  and  t» — >  F(q,t)  is  strongly  continuously  differen¬ 
tiable  in  (0,T)  with  derivative  continuous  in  [0,T],  then, 

[  ,  p.  203] 

i)  z(t)  is  absolutely  continuous  in  (0,T),  where  t  i — >  z(t) 
satisfies  (1.4) 

ii)  z(t)  £  Dom(o/)  for  t  >  0 

iii)  U(t)  -  zQ|  -*■  0  as  t  -*•  0 

iv)  z(t)  satisfies  (1.5);  i.e.,  t  i — =►  z(t)  is  a  mild  solution. 

Also,  the  relation  between  weak  solutions  of  (1.4)  and  mild 

solutions  of  (1.4)  is  given  by  the  following  theorem  [  3,  p.  204], 

If  we  relax  the  assumptions  above  and  only  require  that  F(q,*)  € 

2 

L  (0,T;z),  then  there  exists  a  unique  weak  solution  t  ) — >  z(t;q) 
of  (1.4)  for  0  t  <_  T,  where 

<z(t),5>  is  absolutely  continuous  for  every  £  €  Dom(»*) , 

and 

gY<z(t),C>  =  <  z  (t)  C>  +  <F(q,t),C>  ,  0  <  t  <  T. 

Furthermore,  z(t)  satisfies  (1.5);  i.e.,  t  i — =►  z(t)  is  also  a 
mild  solution. 
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To  define  point  evaluations,  we  note  that  Z  will  always  be 

taken  as  a  function  space  of  Rn-valued  functions  defined  on 

[0,1],  Thus  for  every  q  €  Q,  t  €  [0,T],  we  may  associate  the 

function  u(«,x;-):  [0,1]  Rn  satisfying  u(t,*;q)  =  z(t;q). 

Then,  with  u(t,*;q)  =  (Uj(t,sq) . u  (t,*;q)),  we  define  the 

map  ifQz(t;q)  =  Uj(t,*;q).  In  subsequent  applications,  we  will 

always  take  Z  =  Z,  *  ...  *  Z  where  Z.  consists  of  functions 

1  n  1 

sufficiently  smooth  (in  fact,  subspaces  of  the  Sobolev  spaces  Hm, 
where  m  >  1)  that  point  evaluations  of  u^  will  make  sense. 

Now,  as  we  have  set  up  the  abstract  equation  (1.4)  so  that 
z^(t;q)  corresponds  to  y(t,*;q),  we  may  define  the  cost  func¬ 
tional  J(q , z ( • ,q)  ,n)  by 

(1.6)  J(q,z(* ;q) »n)  =  J(q,*0z(* ,q) ,n) , 

where  J  is  defined  in  (1.3). 

This  leads  us  to  define  in  place  of  (ID),  the  abstract 
identification  problem. 

Definition  1,2.  The  abstract  identification  problem  (IDA)  con¬ 
sists  of  the  following:  find  q  6  Q  <=  rP  such  that 

5(q.z(*  ,q)  ,n)  <  J(q,z(- ,q)  ,n)  for  all  q  €  Q, 

subject  to  z(*;q)  satisfying  (1.5). 

Remark .  This  is,  of  course,  a  reformulation  of  the  original 
identification  problem  in  terms  of  the  abstract  equation,  employ¬ 
ing  mild  solutions.  The  two  coincide  where  classical  solutions 
of  (1.1)  or  (1.2)  exist.  It  is  possible  (see  [11,  pp.  7-8])  to 
formulate  the  cost  functional  J  in  a  more  general  way  so  as  to 
permit  identification  when  data  n  consists  of  measurements  of 
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quantities  other  than  transverse  displacements.  Our  goal  here  is 

to  demonstrate  how  a  numerical  approximation  may  be  constructed, 

based  on  the  form  of  (1.4)  and  the  choice  of  Z,  and  hence  we 

have  restricted  our  attention  to  the  above  problem. 

Since  we  cannot  solve  (IDA)  directly  in  general,  we  describe 

a  general  procedure  for  approximating  solutions  to  (IDA)  following 

the  approach  taken  in  [  7  ]  and  [11]  .  he  take  a  sequence  of  finite 

N  N  N 

dimensional  subspaces  Z  <=  Z  and  define  P  :  Z  ■*  1  to  be  the 
orthogonal  projection  satisfying 

(1.7)  |PNz  -  z  |  |  £  -  z  |  for  all  5  e  ZN,  or  equivalently, 

(1.8)  <PNz  -  z,c>  =  0  for  all  c  €  ZN. 

The  subspaces  ZN  will  be  chosen  so  that  ZN  -►  Z  in  the  sense 

N 

that  P  -►  I  strongly  on  Z ,  We  then  replace  the  abstract  equa- 

v 

tion  (1.4)  on  Z  by  a  sequence  of  approximating  equations  on  Z  : 

(1.9)  zN(t)  =-£/N(q)zN(t)  +  FN (q , t )  for  t  >  0 

N,n.  N,  . 
z  (0)  =  zQ(q) . 


The  approximations  we  use  are 


(q)  =  PV(q)PN 
FN(q,t)  =  PNF(q,t) 

=  pNzo(q)- 

Note  that  this  choice  of  ja^N(q)  requires  that  ZN  c  Dom(o^(q)). 

This  requirement  not  only  imposes  limitations  on  the  smoothness  of 
the  elements  in  ZN  but  also  dictates  that  the  boundary  condi¬ 
tions,  which  appear  in  Dom(j/(q)),  be  satisfied  by  every  zN(t)  €  ZN. 
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Remark .  The  norms  on  the  spaces  Z  used  in  subsequent  ap¬ 


plications  of  this  theory  may  be  parameter  dependent,  in  which 
case  the  projections  p\q)  may  depend  on  the  parameters.  Al¬ 
though  we  do  not  emphasize  this  parameter  dependence  until  we 
encounter  specific  cases,  this  should  be  kept  in  mind. 

Assuming  that  Ji'fq)  is  the  generator  of  a  Cq  semigroup, 
.a/(q)  is  closed,  and  by  the  closed  graph  theorem ,ja^(q)  is  a 
bounded  operator  (see  [  3,  p.  80].  closed  and  B  bounded  im¬ 
plies  is  bounded,  hence  P  j3^(q)P1''  is  a  bounded  operator)  . 


Thus  _o^N(q)  generates  a  semigroup  (T^(t;q)}  given  by 
Moreover,  t  i — >  z^(t;q)  satisfies  (1.9)  on  [0,T]  if  and  only 

if  it  satisfies 


(1.10)  zN(t;q)  =  TN(t;q)PNz0(q)  + 


i: 


TN(t-s;q)PNF(q,s)ds 


Then  standard  Picard  iteration  arguments  imply  that  solutions 
zN(*;q)  of  (1.10)  exist. 

Moreover,  when  _Q^( q)  is  maximal  dissipative  (and  hence  gen¬ 
erates  a  Cq  semigroup  of  contractions),  then  j#'  (q)  will  also 
be  maximal  dissipative,  and  so  ||T^(t;q)||  <  1.  This  follows  from 
the  following: 

«^(q)z,z>  =  <P^(q)PNz,z> 

=  <Mq)PNz,PNz> 

<  0 


and 
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<&/N(q))*z,z>  =  <(PIV(q)PN)*z,z> 

=  <P^  (a^[q) )  *P^z  ,  z> 

=  <c^(q))*PNz,PNz> 

<  0, 

where  we  have  employed  the  self-adjointness  of  the  projection 
operator  PN  and  the  characterization  that  a  dissipative  opera¬ 
tor  is  maximal  dissipative  if  and  only  if  <Q^z,z>  <  0, 

<aAz  ,z>  <_  0  and  s a/  is  closed  [25,  p.  85]. 

We  can  then  formulate  an  approximate  identification  problem. 

Definition  1.3.  The  approximate  identification  problem  (IDN) 
consists  of  the  following:  find  q^'  £  Q  c  rP  such  that 

J(qN,zN(-;qN) ,n)  <  J(q,zN(- ;q) ,n)  for  all  q  £  Q 

subject  to  zN(-;q)  satisfying  (1.10)  on  [ 0 , T] . 

One  of  the  goals  will  be  to  prove  convergence  of  the  solutions 
q^  of  the  approximate  problem  (IDN)  to  a  solution  q  of  the 
identification  problem  (IDA).  In  [11,  pp.  15-17],  it  was  proved 
that  the  map  qi — >  zN(t;q)  is  continuous  for  zN(*;q)  satis¬ 
fying  (1.10).  Since  J  consists  of  point  evaluations  on  the 
first  component  of  zN(-;q),  it  is  easily  argued  that  the  map 
vi — >  J(* ,v, •)  is  continuous  on  C([0,T];z)  (recall  again  that 

we  will  take  1  =  1,  *  ...  *  Z„  such  that  1,  is  a  subset  of  H1 

In  l 

or  H2).  Thus  we  find  that  the  map  q  i — >  J (q , zN ( • ; q) ,n)  is 
continuous,  and  so  for  each  N  there  exists  a  solution  qN  to 
the  approximate  identification  problem  (IDN) .  By  the  compactness 
of  Q,  there  exists  a  convergent  subsequence  (again  denoted  by  q  ) 
such  that  qk  converges  to  some  q  in  Q. 


We  show  that  this  q  is  a  solution  to  the  identification 
problem  (IDA) .  This  will  be  accomplished  if  we  can  show 

(1.11)  lim  jq^-q*|  =  0  implies  1 im | zN (t ;qN) -z^ (t ; q*) |  =  0 

N'veo 

for  every  t  €  [ 0 , T] .  To  see  this,  observe  that  J(q\ zN( • ;q^) ,n) 
J(q,zN(- ;q) ,q)  for  every  q  €  Q,  since  q^  is  a  solution  of 
(IDN) .  Taking  the  limit  as  N  and  applying  (1.11),  we  ob¬ 

tain 

J(q.z(*  ;q)  ,n)  <  J(q,z(- ;q) ,n) 

for  every  q  €  Q.  Thus  q  will  be  a  solution  of  the  identifica¬ 
tion  problem  (IDA)  if  (1.11)  can  be  shown  to  hold.  Since  t  - — > 
z (t ; q)  satisfies 

(1.4)  z ( t ; q )  *  T(t;q)zn(q)  +  [  T(t - s ; q) F(q , s) ds ,  t  €  [0,T], 

u  J  0 

and  t  ■> — »  zN(t;q)  satisfies 

(1.10)  zN(t;q)  -  TN(t;q)PNz0(q)  +  TN (t - s ; q) PNF (q , s) ds , 

t  €  [0,T] , 

an  application  cf  the  Lebesgue  bounded  convergence  theorem  (see 
[11,  p.  20])  yields  a  convenient  criterion  for  showing  when  (1.11) 
holds . 

Proposition  1.4.  Let  t  ■ — >  z(t;q)  and  t  ■ — >  z^(t;q)  be 
solutions  of  (1.4)  and  (1.10)  respectively,  and  assume  (HF)  holds. 
Then 

lim|qN-q*|  =  0  implies  1 im | zN( t ; q^) - z (t ; q*) |  =  0 

JV-MO 


for  every  t  €  [ 0 , T]  if 
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i)  II T^(t ; q^)  ||  <_  Mewt  with  M,a>  independent  of  N,  q. 

N 

ii)  P  -*-1  strongly  in  Z. 

iii)  T^(t;q)  -*■  T(t;q*)  strongly  in  Z  and  uniformly  in 

t  6  [0,T3  when  {q^}  is  any  convergent  sequence  with 


This  proposition,  proved  in  [11,  p.  20]  for  the  more  general  case 
when  F  can  be  mildly  nonlinear,  is  the  fundamental  tool  in  prov¬ 
ing  convergence  of  solutions  of  the  approximate  identification 
problem  to  solutions  of  the  abstract  identification  problem. 
Verification  of  (i)  has  already  been  treated  in  the  case  where 
_3/( q)  is  the  generator  of  a  CQ  semigroup  of  contractions.  In 
that  case,  we  have  already  shown  that  ||T^(t;q)||  <_  1. 

Part  (ii)  of  Proposition  1.4  will  be  verified  for  the  specific 
N 

cases  where  Z  is  the  linear  span  of  certain  cubic  or  quintic 
splines  satisfying  prescribed  boundary  conditions.  This  will  be 
done  in  Section  3. 

Finally,  part  (iii)  of  Proposition  1.1  can  be  established  using 
the  Trotter-Kato  theorem,  vrhich  can  be  viewed  as  a  functional 
analytic  version  of  the  Lax  Equivalence  Theorem  (stability  plus 
consistency  implies  convergence) .  The  version  we  use  is  due  to 
Kurtz  [28] : 

Proposition  1.5.  [28].  Let  (#,  |  •  |  )  and  (#\  |  •  |  )  ,  N  =  1, 
2,...,  be  Banach  spaces  and  let  be  bounded  linear 

operators.  Assume  further  that  y(t)  and  t)  are  linear 
Cq  -  semigroups  on  2  and  3^  with  infinitesimal  generators 
ssf  and  respectively.  If 


¥ 
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i)  lim  |fr^z|v  =  |z|  for  all  z  €  ^, 

N-m»  * 

ii)  there  exist  constants  M,w  independent  of  N  such  that 
||^-N(t)||N  <  Me4"1,  for  t  >  0, 

iii)  there  exists  a  set  <=■  &  such  that  9  c  dom(j/),  9  =  J 

and  (Ag  -£?)3>  =  @  for  some  Aq  >  0  for  which  for  all 

z  €  ^  we  have  lim  |_q?^tt^z  -  =  0, 

N-~°  N 

then  lim  |  J-T^(t)  tt^z  -  z  |  N  =  0  for  all  z  uniformly  in 

N+oo  " 

t  on  compact  subsets  of  [0,»). 

Section  2.  Preliminary  Definitions  of  Spaces  and  Norms 

We  now  introduce  some  of  the  spaces  and  norms  to  be  used  in 
later  chapters  and  collect  some  facts  about  them.  Denote  by  Hm 
the  usual  Sobolev  spaces  over  [0,1)  in  one  dimension.  These 
are 

Definition  1.6. 

Hm  =  ( <t> :  Dm  14>  is  absolutely  continuous  on  [0,1] 
and  Dm<j>  €  L2(0,1)>. 

The  usual  norm  on  Hm  is  1 4>  | 2  =  I  |d^4>|27.  Thus,  H°  =  L2 .  The 

Hm  j  =  0  LZ 

spaces  Hm  with  this  norm  are  Hilbert  spaces  [47,  p.  55;  1,  pp. 

44-47).  Moreover,  the  graph  norm  defined  by 

M2  -  Ul2  *  |d"4>I2 

i/  L 

is  an  equivalent  norm  to  I’l  _  [1  ,  p.  79).  Denote  the  norm  in 

Vi 

H°  by  M0. 

Since  we  will  require  functions  satisfying  certain  prescribed 
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boundary  conditions,  we  will  be  concerned  with  certain  subspaces 
of  Hm.  We  will  also  find  it  convenient  to  weight  the  norm  in 
those  spaces  in  order  to  obtain  dissipative  estimates.  In  parti¬ 
cular,  define  for  a  fixed  a  €  R,  a  >  0 

Hp(a)  =  U  £  Hm :  0(0)  =  4>  ( 1 )  =  0,  (DO)  (0)  =  (DO)  (1)  =  0,..., 

(Dm_10)(0)  =  (Dm*10)(l)  =  0}. 

with  inner  product 

<0,0>  =  <aDm4>  ,Dmip>n  , 

v  r  m ,a  *  O' 

where  <*,->Q  t*ie  H°  CL2)  inner  product.  Denote  the  norm 

induced  by  this  inner  product  as  l*lm  a*  Then  H™(a)  is  a  Hil¬ 
bert  space  and  the  norm  I • I  is  equivalent  to  the  usual  norm 

r  1  1  m ,  a  n 

M  m  on  H™  (cf:  [1,  p.  158]). 

H  2 

We  will  also  need  certain  other  subspaces  of  H  . 

Definition  1.7.  Given  a  >  0,  define 

H2(a)  s  (0  €  H2:  0(0)  =  0(D  =  0} 

H2(a)  =  (0  G  H2:  0(0)  =  0(1)  =  O'  (0)  =  0'd)  =  0} 

H2(a)  =  (0  €  H2:  0(0)  =  0'(0)  =  0} 

2 

where  is  equipped  with  the  inner  product 

<0,0>  2, a  =  <“D24>,D20>0. 

2  2 

Note  that  ( a)  =  Hp(a)  as  defined  above.  When  a  =  1 ,  we 

2  2 

simplify  the  notation  by  writing  h  H^(l). 

2 

Theorem  1.8.  The  spaces  are  Hilbert  spaces  and  the  norm 


» 
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2 

|‘|,  ,  is  equivalent  to  the  usual  norm  |* |  9  on  HT  for 

z ,  i  K 

k  =  1,2,3. 

Proof :  The  case  k  =  2  is  well  known  (see  above)  .  Ke  con¬ 
sider  the  cases  k  =  1  and  k  =  3.  First  we  observe  that  |*|,  . 
u2 

is  a  norm  on  H.  : 

k 

Since  |  -  |  2  ^  is  a  semi-norm  on  H"  [  1,  p.  73],  w'e  only 

2 

need  to  show  J  u ] 2  =  0  implies  u  =  0  for  u  €  H^,  k  =  1,3. 

12  2  ° 

But  | u  j  2  =  0  implies  Iq  (D“u)  *  0,  which  implies  D“u  =  0 

a.e.,  or  u  =  ax  +  b  a.e.;  furthermore,  u,  Du  are  absolutely  con¬ 
tinuous,  u(0)  =  u(l)  =  0  (if  k  =  1)  or  u(0)  =  u'(0)  =  0  (if 
k  =  3) ,  which  implies  a  =  b  =  0. 

2 

1)  k  *  1 .  Ke  first  prov^  is  a  closed  subspace  of  H". 

2 

Suppose  u^  £  is  a  Cauchy  sequence  in  the  | • | 2  norm;  then 

2  7 

u  -*■  u„  €  H.  .  But  u  €  H“  implies  that  u  (0)  =  u  (1)  =  0  for 
n  0  1  n  l  r  n  n 

? 

every  n.  Thus,  ur  -*•  u  in  H  implies  that  un(°)  *  u(0)  ,  and 
u  (1)  -*■  u(l)  ,  and  so  u(0)  =  u(l)  =  0. 

Next  we  prove  |*|2  is  equivalent  to  the  I*Iq  norm  (and 
hence  to  the  j • |  ?  norm. 

2  H  0  2 

Let  u  €  H  Then  there  exists  a  v  €  H  such  that  D  u  =  v 

0.2  2 
and  for  every  v  €  H  there  exists  u  €  such  that  D  u  =  v, 

y  Si  I  J  « 

namely  u(x)  =  /q  /q  v(s)dsds^  -  x  f g1  v(s)dsds^.  Then,  by 
the  Rayleigh-Ritz  inequality  [40,  p.  5] , 

M0  l  ?  lDul0 

=  \  D(  1  vfsjdsdsj  -  x  f  f  1  v(s)dsds1)  Q 
0  J  0  1  J  0  ■*  0 

1  fX  I 1  fsl 

=  -  v(s)ds  -  v(s)dsds1  n 

n  Jo  Jo  Jo  1  u 
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|f  |  1  vCsjasdsJJ. 


v(s)ds  Q  + 


0  '0 


Since  v  €  ,  Schwartz's  inequality  implies  /Qjv(s)jds  <_  |  v  |  ^ . 

Define  g(x)  -  /g|v(s)|ds  <  /J|v(s)jds  <  j v | Q  for  all  x  €  [0,1]. 
Let  c  s  [ v | 0 .  Then  /*  g(s)ds  <  /J  c  =  c.  So,  continuing  the 


chain  of  inequalities,  we  obtain 


Milo  vtsH°  *  Ho  lo1  vtS)dSdSl|o} 
MlDvCs,|Ho  *  ll„  | 1  ivcsi idsdSl 1 0} 


u '  0  i 


o  >  o 

r  1  rS 


|v(s) Ids  0  + 


0  J0 


£  I c I n  +  i c  I  n>  =  b  =  l 


0  J  TT  TT  ,V|(T 


But  v  =  D“u  a . e . ,  so 


ulo  -  ^v|o  =  l|r)2u| 


0  =  IT  1  u  1  2  ’ 


I  D^u| q  <  I D2u | q  +  |u|‘  £  |D4u|p  +  -~|D"u 


2,2  .  4  ,„2  ,2 


-  a  +  -V) i o2ui o» 

TT  ^ 


I  u  I  2  £  I  u  I  G  £  (1  +  — y)  f  u  i  2  * 

7T 

Now,  since  the  graph  norm  |*|r  is  equivalent  to  the  usual 

2  b 
7  7  i  2  2 

H"  norm  |-|  7  =  £  |DJu|n  on  H  ,  we  have  equivalence  of  |-i 

II"  *  _  n  vJ 


H  1  = 
i  to  |  • 


on  HJ 


2)  k  =  3.  The  proof  is  similar  to  the  previous  case  (k  =  1), 

once  we  observe  that  u  €  if  and  only  if  there  exists  a  v  e 

2  x  S1 

such  that  D  u  =  v,  u(0)  =  u’(0)  =  0,  namely  u  =  /Q  / Q  v(s)dsds1. 

The  domains  of  the  operators  we  will  study  will  be  certain 


dense  subsets  of  satisfying  particular  boundary  conditions, 
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Definition  1.9.  We  say  $  satisfies  boundary  conditions  of 
type  k,  k  =  1,2,  or  3,  if 


a) 

o 

V 

II 

o 

& 

=  *(1)  =  4”(1)  =  0 

for  k  =  1, 

b) 

4.(0)  =  ♦,(0) 

=  4>(1)  =  4> '  (1)  -  0 

for  k  =  2 , 

c) 

♦  (0)  =  <D'(0) 

=  *"(1)  =  <T'(1)  = 

0  for  k  =  3. 

Define 

hJ  =  (M  H4 

:  <f>  satisfies  boundary  conditions  of  type  k) 

Note  that  for  the  Euler-Bernoulli  beam,  a)  corresponds  to  the 
boundary  conditions  for  the  simply  supported  beam,  b)  to  boundary 
conditions  for  a  beam  clamped  at  both  ends,  and  c)  to  boundary 
conditions  for  a  cantilevered  beam. 

We  will  also  need  to  use  the  following  fact: 

is  dense  in  H°  (in  the  |*|q  norm). 

This  follows  from  the  fact  that  H™  can  be  characterized  as  the 

closure  of  Cq (0,1)  in  the  Hm  norm  [1,  p.  44;  30,  p.  91], 

and  that  c  H™  c  h£  c  for  m  ^  4;  the  result  follows 
00  *0 

since  is  dense  m  H  [1,  p,  31]. 

Section  3.  Splines  and  Error  Bounds  for  Approximation  by  Splines 

N 

The  approximating  subspaces  Z  c  Z  will  consist  of  sets  of 
cubic  or  quintic  splines  satisfying  prescribed  boundary  condi¬ 
tions.  We  begin  this  section  by  defining  a  particular  class  of 
interpolating  functions,  the  L-splines. 

Given  any  partition  =  (x^J^.q,  0  =  xo  <  X1  <  •  •  •  <  XN  =  1 » 
the  L-splines,  as  defined  by  Schultz  and  Varga  [41],  are  piece- 
wise  continuous  functions  satisfying  certain  interpolating  con¬ 
ditions.  The  functions  themselves  belong  to  the  kernel  of  certain 
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differential  operators  on  each  [x^,x^+^].  To  be  more  precise, 
given  any  differential  operator  L  of  order  m 

m 

(1.12)  Lu(x)  =  l  a  •  (x)  D-’u(x)  ,  m  >  1, 

j  =  0  J 

where  a.  €  C-1  [a,b]  ,  and  given  any  partition  we  defin.  an 

incidence  vector  d  =  (d^,...,d^T  of  positive  integers  with 
1  <  d-  <  m  for  1  <_  i  <_  N-l,  and  define  the  L-spline  space  in 
the  following  way: 

Definition  1.10.  [41].  The  L-spline  space  Sp(L,A^',d)  is 

the  collection  of  all  functions  on  [a,b]  such  that 

L*Ls(x)  =0,  x  €  [a,b]  ^  AN 

D^s(x^)  =  D^s(xp  for  all  0  k  2m-l-d^,  0  <  i  <  N, 

where  L*  is  the  formal  adjoint  of  L. 

We  will  be  interested  in  the  specific  case  where  L  =  Dm  and 
d^  =  1,  0  <  i  <  N.  In  this  case,  the  functions  in  Sp(L,A^,d) 

7  m  -  7 

are  piecewise  polynomials  of  degree  2m-l,  with  C  [0,1]  con¬ 
tinuity.  In  order  to  obtain  error  estimates  for  projections  onto 
spaces  of  splines,  we  use  well-known  error  estimates  for  inter¬ 
polation  by  L-splines.  For  our  purposes,  we  require  the  type-1 
interpolant  of  [41]. 

Definition  1.11.  [41].  Given  f(x)  £  Cm_1[a,b]  and  with  L 

the  differential  operator  of  order  m,  AN,  and  d  as  above,  we 
say  a  function  s(x)  £  Sp(L,AN,d)  is  a  Sp(L ,AN ,d) - interpolant  of 
f(x)  of  type  1  if 
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i)  s(xi)  =  f (xi) ,  0  <  i  <  N 

ii)  (Dks) (x^)  =  (Dk£)(xi),  0  <  k  <  m-1  for  i  =  0  and  i  =  N. 

-I 

It  is  known  141]  that,  given  any  f  €  l  there  exists 

N 

a  unique  s  6  Sp(L,A  ,d)  such  that  s  is  a  type-1  interpolant 
of  f . 

For  each  positive  integer  N,  let  A^'  =  x^  =  i/N; 

in  the  sequel,  we  will  restrict  our  attention  to  such  uniform 
partitions  of  [0,1].  Also  we  take  L  =  Dm  and  d^  =  1,  i  =  1, 
...,N-1.  The  corresponding  spline  subspaces  will  be  denoted 
S^m  ■*'(AN)  =  Sp(Dm,A^,d).  Thus  a  spline  of  degree  2m-l  (order 
2m)  will  be  a  function  s(x)  €  S2m'1(AN)  c  c2m'2  such  that  s(x) 
is  a  polynomial  of  degree  2m-l  (order  2m)  on  each  subinterval 
[Xi,xi+i]  defined  by  the  partition  A^.  In  particular, 

•7  -iw»  ^ 

S"(A1N)  =  (s(x)  6  L  [0,1]:  s  is  a  polynomial  of  degree  3 

on  [xA  ,xi  +  1] ,  i  =  0 , N - 1 } 

is  the  set  of  cubic  splines  (m  =  2)  and 

S^(A^)  =  (s(x)  £  [0,1]  :  s  is  a  polynomial  of  degree  5 

on  [xi,xi+13 ,  i  =  0 ,N-1 } 

is  the  set  of  quintic  splines  (m  =  3).  The  type-1  interpolant  s 
to  f  €  C*(0,1]  from  S^(*^)  satisfies 
i)  s(xi)  =  f (Xj) ,  0  <  i  <  N, 

and 

ii)  s (0)  =  f(0),  s(l)  =  f(l) 

(Ds)(0)  =  (Df ) (0) ,  (Ds) (1)  =  (Df ) (1 ) . 

Given  f  €  C“[0,1],  the  type-1  interpolant  s  to  f  from  S5(A^') 
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satisfies 

i)  s ( x i )  =  f(xi) ,  0  <  i  <  N 

ii)  s (0)  =  f  (O')  ,  s  ( 1 )  =  f  Cl) 

CDs)  (0)  =  (Df )  ( 0)  ,  CDs)  CD  =  CDf)  (1) 

(D2s)(0)  =  C  D 2  f ) CO) ,  (D2s) (1)  =  (D2f) (1) . 

These  finite-dimensional  subspaces  S2m  ^(A^)  of  H2m  ^ 
possess  a  convenient  basis  (the  B-splines)  which  will  be  defined 
at  the  end  of  this  section.  Suffice  it  to  say  for  the  moment  that 
the  B-spline  basis  functions  have  small  (local)  support  which 
make  the  spaces  S2m  ^(A^1)  efficient  for  Galerkin-type  approxi¬ 
mations.  The  second  advantage  of  splines  is  their  approximation 
power . 

Since  we  shall  be  approximating  functions  using  projections 
of  the  form  PN:  Hk  •*>  S2m'1(AN),  k  €  (0, .  .  .  tm-l)  ,  we  shall  need 
to  obtain  error  bounds  for  such  approximations  by  splines.  The 
following  result  for  bounding  the  error  in  spline  interpolation 
will  be  fundamental  in  deriving  error  bounds  for  projections  onto 
subspaces  of  splines.  While  the  results  of  this  theorem  have 
been  proved  for  general  partitions  (satisfying  a  uniformity  condi¬ 
tion  max  (x-.,-x-)/  min  (x- -x-)  <  a,  some  a  >  1)  and 
0<i<N-l  11  1  0<i<N-l  11  1  “ 

for  the  general  L  in  (1.12),  we  state  the  results  only  for  the 
particular  case  we  need. 

Theorem  1.13.  [44].  Given  f  6  H2m,  let  s  be  the  unique 
element  in  S2m  ^(AN)  which  interpolates  f  in  the  sense  of  a 
type-1  interpolant.  Then, 


|DJ(f-s)|0  £  cj,m(|)2m'jlD2mfl0.  j  =  0 . 2m  - 1 , 
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where  c-  are  constants  independent  of  f. 

J  »m  * 

Remark .  This  result  is  proved  for  j  =  0,...,m  in  [41, 
Theorem  9].  For  the  case  where  m  +  1  <_  j  <  2m-l,  see  [44,  Lemma 
3.1]  and  the  references  there.  We  will  also  need  the  Schmidt 
inequality  [40,  p.  7]: 

Lemma  1.14.  If  Pn(x)  is  a  polynomial  of  degree  n  =  1,2, 

...,5,  then  /^[Dp  (x)l2dx  <  c  (b-a)  u  /^[p  (x)]^dx,  where 
a  n  —  n  a 

c1  =  12,  c 2  =  60,  c3  =  2(45  +  /T60T)  =  170.,  c4  =  440,  c&  =  738.8  . 

Now  we  wish  to  obtain  error  bounds  for  projections  onto 
certain  subspaces  of  S3(A^)  and  S3(A^)  consisting  of  splines 
satisfying  prescribed  boundary  conditions.  In  [7],  we  proved 
the  needed  results  for  a  certain  subspace  of  cubic  splines,  and 
we  restate  those  results  here. 

Denote  by  S3(AN)  =  [s  €  S3(AN):  s(0)  *  s(l)  =  0}.  Given  a 
function  z  €  H®,  denote  by  P^'z  its  projection  onto  S3(A^) 
in  the  H°  norm.  Pq  is  the  map  Pq  :  H°  — ^Sq(AN)  satisfying 
<z  -  PqZ,s>q  =  0  for  every  s  €  S3(A^),  or  | PqZ - z I q  = 
inf  {  |  s  -  z  |  q  |  s  €  S0V)}.  In  [7],  we  prove 

Lemma  1.15.  [7,  2.3].  If  z  £  hJ  n  H4  e  H°,  and  pjj  is 

the  projection  P^:  H°  — >  S3(AN),  then 

I  D-*  (PqZ-z)  |  0  <  <0,j(^4'j|D4zi0’  j  =  °’1’2- 

Take  (ot)  to  be,  as  defined  in  Section  2,  equipped  with  the 

|  •  J 1  a  norm  induced  by  the  inner  product  a  =  <aD<)>  ,D\J>>q  . 

1  N  * 

Given  a  function  z  €  H^,  denote  by  P^  z  its  projection  onto  the 
subspace  S3(A^).  P^  is  the  map  P^:  Hq(<x)  — =►  S3(AN) 
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satisfying  <2  -  P^z,s>^  a  =  0  for  every  s  e  Sq(A^),  or 
|P^z-z|^  a  =  inf(|s-z|^  a  |  s  €  S3(AN)}.  In  [  7],  we  also  prove 

Lemma  1.16.  [7,  4.3].  If  z  €  H*  D  H4  and  P^  is  the 

projection  P^:  — =*•  S3(A^)  in  the  |*|^  Q  norm,  then 

|D^  (Pjz-z)  |ljQt  <  /o  k2  j  |D4z|0.  j  =  0,1,2. 

Note  that  P^  and  P^  are  both  orthogonal  projections  onto 
Sq(Aa),  but  the  projections  are  taken  with  respect  to  different 

norms.  Also,  since  H4  is  dense  in  H1  and  in  H°,  the  above 

error  bounds  hold  on  dense  subsets  (H4  fi  Hq)  of  the  domain  of 
definition  of  the  projection  operators. 

We  wish  to  state  and  prove  analogous  results  for  projections 

onto  spaces  of  quintic  splines  satisfying  prescribed  boundary 

5  N 

conditions.  Denote  by  Sk(A  )  the  collection  of  quintic  splines 
satis.*'  j.ng  boundary  conditions  of  type  k,  defined  in  Section  2: 


=  {s  €  S5(AN):  s  satisfies  boundary  conditions  of 
type  k ,  k  =  1,2,3}, 

=  (s  €  C4  [0,1] :  s  is  a  quintic  polynomial  on  each 
Cxi , xi + 1 ] >  i  =  0,...,N-1,  and  s  satisfies  the 
boundary  conditions  of  type  k,  k  =  1,2,3}. 

We  shall  again  require  two  sets  of  projections.  Those  we  need  are 


and 


PN  •  H° 
2,k‘  M 


C5,.N, 
sk(4  ) 


P3,k: 


-5,.N, 
Sk(4  ) 


] 

) 


First,  we  obtain  some  preliminary  bounds  using  the  Schmidt  inequal¬ 
ity. 
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Lemma  1.17.  Let  s  €  S^(AN') 


boundary  conditions  fo  type  1, 
i)  |  s  |  Q  <  CjNj  |sL,  1  < 


ii)  Id2s!2i0  <  =SN 


^  2 


2,0’ 


be  a  quintic  spline  satisfying 
2  or  3.  Then 

j  <  4 


and 

iiij  I s I  2 , a  -  ^  c6N2|s|0- 


Proof :  Apply  the  Schmidt  inequality: 


1 

| Ds (x) | ^dx 

0 


fx.  .  N 

1  |Ds(x) | 2  <  l  c,N 

x,  i=l  5 


1  | S (x) I  2  , 

Jxi-1 


since  s(x)  is  a  polynomial  of  degree  5  on  [x^  j,x.J,  i  =  1, 

This  implies  i)  with  ?1  =  /cs,  j  =  1.  Let  g(x)  =  Ds(x). 
Then  g  is  a  polynomial  of  degree  4  oi  each  tx^^.x^],  i  =  1, 

•  .  .  ,N,  a.rd 


f1  7  7  f1  ,  N 

j  |D  s(x) |  dx  =  | Dg(x) |  ax  <  l 

J0  J0  ~  i  =  l 

=  c.N2  f  |  g(x) | 2dx 
J  0 

■  c.N2  f  | Ds (x) | 2dx 
4  ;0 


1  4 


1  c4i:lN 


rxi  ? 

I  g  Cx) | 2dx 

xi-l 


which  establishes  (i)  with  c2  =  /c4  c5,  j  =  2.  In  other  cases 

(j  =  3,4)  follow  in  a  similar  manner.  Inequality  (ii)  is  derived 

2 

in  the  same  manner.  Let  g(x)  =  D  s(x);  then  g  is  a  cubic  poly¬ 
nomial  on  [x^j.x^I,  i  =  1 , . . . ,N ,  and 
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1 1)2  5  t  2  ,  cx  =  “l^glo  =  ai  f1  D2gfx)dx 

1  =  1  ^  xi - 1 
N  fX . 

=  »  I  1  D(Dg(x))dx 
1=1  ixi-l 

?  N  f x  • 

a  c~Nz  I  1  Dg (x) dx 

i=1  Jxi-1 


.  N  fX . 

a  c9c,  N4  l  I  1  g(x)dx 


i  =  1  J  x  . 

1  1  xi-l 


=  a  c2c3N4|g|2 


=  c2c3N4|s|2>  , 


which  proves  (ii)  with  c,.  =  /c2c3.  Finally 


*i!.«  ■  “i^io  - «.?.  r* 

1=1  ;x. 


D(Ds (x) )dx 


N  rX. 


Ds (x) dx 


<_  a  c.N  £  Ds(x)dx 

i=1  Jxi-1 

<  a  c^CgN4  I  [  1  s  (x)  dx 
1=1  xi-l 


a  ^4^5^  I s I o  > 


which  establishes  (iii)  with  c^  =  /c^c  3 .  □ 


*  5  f  .  h  x 


Finally,  we  bound  the  error  for  projections  onto  S^(A  )  for 
k  =  1,2. 

Lemma  1.18.  Let  P2  ^  be  the  orthogonal  projection  P2  j,: 

H°  — >  sjj(AN)  with  respect  to  the  H°  norm,  for  k  =  1  or  2 ; 
then  if  z  €  (4>  €  H^:  <J>  satisfies  boundary  conditions  of  type  k) 

^  iP2,kz'zto  -  K2,0%^6  >d6zI0 
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ii)  ^(P^z-z)  !0  1  K2,j(F)6’j,I)6z|0’  5  = 

and 

lil-)  ^P2,kZ‘Z^2,a  -  ^  k2,2%^  lD  z^0* 

Proof :  The  first  estimate  (i)  follows  directly  from 

I P2 >kz'z I o  i  I Ikz"z I  0 

where  I^z  is  the  type-1  interpolant  to  z  from  S^(A^).  Since 

z  satisfies  the  boundary  conditions  of  type  k,  and  tvpe-1  inter- 

polants  interpolate  these  conditions,  it  follows  that 
N  5  N 

I^z  £  Sj,(A  ).  Furthermore,  the  basic  spline  interpolation  error 
bounds  (Theorem  1.13)  hold,  and  thus  (i)  follows  directly. 

To  obtain  (ii),  we  write 

(1.13)  |Dj(P^kz-Z)|0  <  |Dj(P^k-I^)l0  .  |Dj(I^-Ol0. 

A  bound  on  the  second  term  in  (1.13)  follows  directly  from  the 
interpolation  error  bounds: 

lDjOk*-Ol0  I  c_1i3(l)6-jlD6z|0. 

To  bound  the  first  term  in  (1.13),  we  again  resort  to  the  Schmidt 
inequality,  since  ^z  -  I^z  €  S^(AN), 

lD^  (p2,kz_Ikz^0  -  lp2,kz"1kz^0’  by  Lemma  1-17» 

i  Vl^-lo 

<_  2cjN^  6  c0  3 1  D6z  |  q,  by  Lemma  1.13, 
which  gives  us  (ii)  with  k~  .  =  2c. cn  ,  +  c.  Finally 
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lP2,kZ'Z^2,a  =  ^P2,kZ'Z^0  -  ^  K2,2fN)  ZI(T 


The  final  projection  error  bounds  which  we  shall  require  are 

2 

those  for  the  projections  in  the  H^Ca)  norm.  Recall  that 
2  2 

=  {$>  £  H  :  <J>  satisfies  boundary  conditions  of  type  k, 

2  2 

k  =  1,2,  or  3},  and  that  Hk(a)  is  equipped  with  the  norm 

2  2 

| • | 2  a  derived  from  the  inner  product  <‘»’>2  a  E  <aP  * ^O' 

We  define  k  to  be  the  orthogonal  projection  Hk(a)  — > 

S5(AN)  satisfying 

|P3,kz-zl2,a  =  inf { | s-z | j ,a  I  s  €  sJ(iN)). 

N  N  2 

Lemma  1.19.  Let  P‘3  k  be  the  projection  P^  Hk(a)  — > 

taken  with  respect  to  the  j*|2  Q  norm.  If  z  e  (<J>  £  H6  : 

<J>  satisfies  boundary  conditions  of  type  k}  and  k  =  1  or  2 , 

then 

^  lP3,kZ‘Z  l2,a  -  ^  k3,0^N-)  I D  z  I  0 
ll)  lD  ^P3,kZ‘Z^  l2,o  -  ^  K3,2(‘N)  zI0 
iii)  | D  (P3>k2-Z)|0  i  <3>2(-N-)  lD  zi(T 

Proof :  Let  I^z  be  the  type-1  interpolant  from  S5(a^)  to 

z.  Then  since  type-1  interpolants  interpolate  the  boundary  condi¬ 
tions  of  type  k  (for  k  *  1  or  2),  I^z  e  S^(A^) ,  and  the  spline 
interpolation  error  bounds  of  Theorem  1.13  hold.  The  inequality 
(i)  follows  directly  from 


iP3,kZ*Z  l2,ct  -  =  (IkZ'Z^0  -  ^  C2 , 3  4  1  d6z  I  0  * 


To  obtain  (ii),  write 


29 


(1.14)  U>V!|,k*-Ol2>0,  i  lI)2(P3,kz-Ik^ 

The  second  term  of  (1.14)  satisfies 


2,a  * 


D2  CI^z-2) 


2,0  *  ^|D4C^-OI0  1  ^  c4>3(i)2|D6z 


from  the  spline  interpolation  error  bounds  (Lemma  1.13). 

To  bound  the  second  term  we  again  use  the  Schmidt  inequality 
and  the  fact  that  ^z,  -  I^z  €  S^(A^)  : 

lD2(P3,kZ‘IkZ)l2,o  -  a5N2|P3,kz-Ikz|2,o-  from  Lemma  1’14' 


-  £SN2(|P3,kz'z|2,o  *  iz'Ikz|2,o) 


<  2ccN  I  z-IjZ  I 


-  5 


k  1  2  ,ot 


=  2  v'ST  c5N2  JD2  (I^z-z)  |Q 

2  &  c5(g)  C-2,1  lD  z  lo 

^  *3,2  W  lD  z  1 0  * 


with  ^3 ^ 2  C4,3  +  2c5c2,3’  0 

Remark .  Projection  error  bounds  for  projections  onto  splines 
satisfying  cantilever-type  boundary  conditions  (type  3)  have 
been  excluded  from  the  above.  The  reason  for  this  is  that  they 
do  not  fit  into  the  above  framework  since  the  type-1  interpolant 
does  not  interpolate  these  conditions.  Furthermore,  the  proof  of 
the  interpolation  error  bounds  requires  that  the  interpolating 
spline  to  f  satisfy  the  first  integral  relation  (cf:  [48]) 


f 1 ( Lf ) 2dx  =  f1 
Jo  Jo 


(L(f-s) }  +  (Ls)  dx , 


2 
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where  L  =  D"3  here.  For  this  to  hold,  it  is  necessary  that 


fl 


(D3s)2dx  =  (-1) 


rl 


s(D  sldx. 


But  the  integration  by  parts  fails  for  type  3  boundary  conditions. 
However,  our  numerical  results  suggest  that  the  projections 


2,k 

and 

3,k 

when  k  = 

3  satisfy  error  bounds  similar  to 

those 

for 

p^ 

2,k 

and  Pi!  , 

3 ,  k 

when  k  =  1,2. 

In  addition  to  their  strong  approximation  properties,  the 
2m -1  N 

spaces  S  (A  )  are  attractive  computationally  because  they 


possess  a  computationally  simple  basis  of  small  (local)  support, 

namely  the  B-splines.  The  B-splines  are  compactly  defined  as 

differences  of  the  truncated  power  basis  functions. 

All  the  B-splines  are  translates  of  one  basic  B-spline. 

Define  B  (x)  =  <5n(x-y)^  1  =  [  (-1) 1  (^)  (x-i)1?  \  as  the  basic 

i  =  0 

B-spline  with  knots  at  0,1 . .  where  6n  is  the  forward 

difference  operator  6f(xQ)  =  ffXj)  -  f  (xQ)  , .  .  .  ,<5k  +  1f  (xQ)  = 
6kf(x^)  -  6^f(xQ),  and  where 

ni  ,  (x-t)n  ,  when  t  <  x 
(x-t);-1  = 

0  ,  when  t  >  x 


is  the  truncated  power  basis  function. 

To  define  the  B-spline  basis  for  S2m'1(AIs'),  take  n  =  2m, 
and  define 

gi,n-l(x)  = 

where  h  =  x^  +  1  -  x^  and  where  we  have  defined  an  extended 
partition  AN  =  xA  =  ih,  of  AN,  and  yi  =  xi  -  mh . 
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N 

Remark.  With  the  B-splines  defined  in  this  manner,  B.  7 

1 , 3 

agrees  with  the  Prenter  definition  of  cubic  splines  [35]  .  Its 

relation  to  the  normalized  B-splines  of  DeBoor  (cf:  [42,  p.  135]") 

is  given  by  B^  =  7 — N,ri(x). 

&  i,n  ( m  - 1 ) !  l v  J 

A 

The  support  of  each  B-spline  B^  is  the  interval 

[x-  ,x.  ],  and  it  is  easily  shown  [42,  p. 116]  that  the  N+m+1 

l-m  l+m  '  r 

r,  ,  •  \N+m-l  r2m-l,.N\ 

B-splines  (B.  7  ,  span  S  (A  ) .  Moreover,  since  all 

l , Zm-1  1=1 -m  1  v 

B-splines  are  obtained  by  scaling  and  translating  one  basic  B- 

spline,  it  is  efficient  to  store  a  function  s  €  ^(A^)  by 

storing  its  coefficients:  s  €  S2m  "'"(A^)  implies  s(x)  = 

N+m-1  N 

7  c.B.  t  (x) .  The  computational  aspects  of  B-splines  will 
i=l -m  1  1»2m-1 
be  explored  in  Chapter  4. 

Finally,  recall  that  we  require  approximating  subspaces  in 
Dom(jtf) .  Thus  we  will  require  our  spxines  to  satisfy  the  boundary 

T  N 

conditions.  For  the  cubic  splines  in  Sq(A  ),  we  require  s(0)  = 
s  C 1 )  =  0.  We  may  take  from  the  N+3  basis  elements  in  S^(A^) 
the  following  N+l  basis  elements  for  Sp(A^): 


rN  -  fiN 
Bi  ,  3  "  Bi,3’ 


,N 


;N 


i  =  2 , . . . ,N- 2 

N 


C1*15)  B0 , 3  =  B0 , 3  *  4B-1,3;  B1 , 3 
BN-1 , 3  =  BN-1,3  '  BN,3/4; 


Bl,3  '  B0,3/4 


,N 


iN 


BN , 3  "  , 3  ‘  4BN+1 , 3 * 


The  same  approach  can  be  used  to  obtain  a  basis  for  quintic  splines 
satisfying  boundary  conditions  of  type  k.  These  are  listed  as 
follows : 

i)  To  obtain  a  basis  for  S^(AN),  take  from  { B^’ 
the  following  N+l  splines  {B^}^_q  (dropping  the  second  sub¬ 
script) 
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Bi 


-v 

Bi> 


3  <  i  <  N-3 


(1.10) 


v  ~\t  aV 
B0  =  B0  '  3B-1  +  1 2B - 2 

bJ  =  B^  -  BN2 

K  /s  \  *  V 

B2  =  B2  ‘  B-2 


with  B N_j(x)  =  B.(l-x),  j  =  0,1,2, 


5,.N, 


ii)  To  obtain  a  basis  for  $2^  ),  take 


(1.17) 


rN  _  pN  33  gN  ,  165  gN 

B0  '  B0  ‘  T  B-1  1“  B  -  2 

N  _  g.N  _  9  gN  ,65  gN 

B1  '  B1  4  B-1  +  T  B -  2 

N  ^ V  l  aV  g  aV 

B2  =  B2  -  F  B-1  +  A  -2 


with  B  ^  =  B?,  3  <  j  <  N-3  and  B^.  (x)  =  bJ’(I-x),  0  <  j  <  2 

iii)  For  the  cantilevered  beam,  we  obtain  a  basis  for 

5  \  \  v 

S3(A‘ )  by  taking  Bq,...,B^_-  as  in  ii)  and 


(1.18) 


BN  =  BN  +  !  BN.l  +  3BN+2 

V  /V  \T  aV 

Vl  ■  BN-1  2BN+2 

rN  _  gN  1  RN 

BN-2  BN-2  2  BN+1 


Having  defined  the  IDA  problem  and  its  approximation  IDN,  and 

having  defined  the  spline  spaces  necessary  to  set  up  the  approxi- 

N 

mating  subspaces  Z  ,  we  now  proceed  to  some  concrete  applications 
of  these  ideas. 


► 
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CHAPTER  2.  APPLICATION'  TO  THE  EULER-BERNOULLI  EQUATION 


Section  1.  The  Euler-Bernoulli  Equation  with  Structural  and 
Viscous  Damping 

Having  introduced  the  basic  ideas  for  the  parameter  estima¬ 
tion  problem,  we  turn  to  the  application  of  these  ideas  to 
specific  equations  arising  in  elasticity.  In  this  chapter,  we 
discuss  the  Euler-Bernoulli  equations.  The  well-known  equations 
for  the  transverse  vibrations  of  a  thin  elastic  beam  are 

El 

3x 

(2.1)  2  2 

3  y  3  LJ( 

m  — l  =  f(t,x;q) 

3t  3x 

where  J(  is  the  bending  moment,  m  is  mass  per  unit  length  and 
f  is  the  applied  load.  Ke  wish  to  include  two  types  of  damping 
in  the  above  model.  The  first  is  velocity-proportional  viscous 
damping  yy  .  The  second  is  structural  damping  arising  from  a 
simple  viscoelastic  model,  the  Kelvin -Voight  model,  where  we  have 
the  constitutive  relationship  o  =  Ee  +  ce,  where  o,e  are  the 
linear  stress  and  strain,  E  is  Young's  modulus,  and  c  is  the 
damping  coefficient.  Thus  we  introduce  damping  proportional  to 
strain  velocity. 

Following  the  usual  development  of  the  Euler-Bernoulli  equa¬ 
tion  (14,  pp.  295-302],  we  obtain 


jK* 


cydA  =  El 

3xZ 


+  cl 


3x  3t 


and  the  equation  (2.1)  becomes 
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(2.2)  >tt  =  -aD4y  -  <5D4yt  -  Yyt  +  f,  t  >  0,  x  £  [0,1], 

where  a>0,  &  >_  0 ,  y  >_  0 . 

In  setting  up  the  operators  and  appropriate  spaces  below,  we 
shall  consider  only  homogeneous  boundary  conditions  of  type  k 
introduced  in  Definition  1.9.  This  results  in  no  loss  of  general¬ 
ity  of  the  methods,  since  any  non -homogeneous  conditions,  includ¬ 
ing  time  dependent  conditions  (which  can  occur  due  to  applied 
moments  and  shears  at  the  boundary)  may  be  transformed  (see  [29]) 
in  a  manner  that  reduces  the  equations  to  ones  with  homogeneous 
boundary  conditions  and  an  additional  "load"  term  included  in  f. 

Within  the  framework  we  have  established,  a  variety  of  ap¬ 
proximations  may  be  used.  Recall  that  we  rewrite  the  boundary- 
initial  value  problem  as  the  abstract  equation 

z(t)  =  ja^(q)z(t)  +  F(t;q)  on  Z 

2(0)  =  ZQ 

and  then  approximate  it  by  a  sequence  of  equations  of  the  form 

zN(t)  =  -3^N(q)zN(t)  +  FN(t;q)  on  ZN 

N,n,  N 
z  (0)  =  zQ. 

One  way  of  obtaining  different  approximations  clearly  is  by  making 

different  choices  for  ja/N.  We  discuss  only  the  choice 
N  N  N 

Srf  =  P  ja^P  here,  but  other  choices  are  possible.  For  example, 
if  o?  2  exists  as  a  differential  operator,  then  p^Q^p^ja/^p^ 
can  be  used  (see  [27]). 

Additionally,  the  choice  of  state  spaces  and  the  iwm  of  the 
evolution  equation  on  these  spaces  leads  to  different  natural 


approximations;  it  is  this  difference  we  wish  to  examine  here. 
In  particular,  for  the  equation  (2.2),  we  shall  write 
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(2.3) 


z(t)  =  V(q)z(t)  +  F(q,t) 
z(0)  =  zn 


-Qf(q) 


-q^ 


■q9D  -q. 


^  in  X  =  H2(a)  x  h° 


or,  (when  y  =  0) 


(2.4) 


r(q)  z(t)  =j/(q)z(t)  *  F (q , t ) 

z(0)  =  zQ 


ss'i  q) 


r(q) 


0  1  0 

0  0  -q1D2 

0  D2  0 

10  0 
0  1  q2D2 

0  0  1 


l  m 


H°  x  h°  x  h° 


where  q:  =  a,  q2  -  6,  q3  =  y. 

As  we  shall  see,  (2.3)  leads  naturally  to  quintic  spline 
approximations  which  are  discussed  in  section  (2)  and  (2.4)  leads 
to  cubic  spline  approximations  which  are  discussed  in  section  (3). 
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Section  2.  An  Approximation  Using  Quintic  Splines 

In  order  to  investigate  and  approximate  solutions  of 


(2.4) 


ytt  ”  ’ayxxxx  *^yxxxxt  +  ^ 


y(0,x)  =  <Hx) 
yt(0,x)  =  iC(x) 


where  y  satisfies  boundary  conditions  of  type  k,  we  write 

2  0 

(2.4)  as  abstract  equation  in  a  subspace  of  H  *  H 


(2.5) 


z(t)  =j*'(q)z(t)  +  F(q,t),  t  >  0 
z (0)  =  zn 


with  z(t) 
and  F(q,t) 


(z,(t) ,z2 (t) )T  =  (y(t,-),yt(t,-))T,  z o 
(0,f(t, • ;q))T,  and 


4q) 


o 

-qxD4 


(*,U0T, 


where  qj  i  a,  q^  =  6.  In  particular,  for  each  of  the  boundary 

2  0 

conditions  of  type  k,  k  =  1,2,3,  we  consider  (2.5)  in  *  H  , 
and  to  be  of  the  form  Ja/,  with  Dom(ja^)  =  *  H4. 

We  derive  the  important  properties  of  these  operators  which 
will  be  used  to  prove  convergence  of  the  approximate  identification 
problem  associated  with  (2.5). 

We  first  consider  the  special  case  where  6=0  (no  damping) 
and  denote  the  corresponding  operators  by  (i.e.,  is 

43^  in  case  q2  =  0).  In  this  case,  we  can  define  the  operators 

/s 

ja^  in  such  a  way  that  maximal  dissipativeness  can  easily  be 

A 

argued  in  the  following  way:  we  will  show  that  J#  ^  and 

are  disspative,  and  stf ^  is  closed;  maximal  dissipativeness  then 


follows  [25,  Theorem  4.4,  p.  87]. 

To  extend  the  results  to  the  case  6  >  0,  we  will  also  need 
to  prove  dense  inclusions  of  certain  subsets  of  Sobolev  spaces. 

To  do  this,  we  will  use  the  following  lemma: 

Lemma  2.1.  Let  a/:  (Dom(^)  <=  X)  ■+•  X  be  linear,  X  a  Hil¬ 
bert  space.  If  Si/  satisfies  the  dissipative  inequality 
<si/x,x>  <_  0  for  every  x  £  Dom(_oO  and  if  R {Si/-  XI)  =  X  for 
some  X  >  0,  then  Dom(j^)  is  dense  in  X. 

Proof :  Suppose  Dom(_oO  is  not  dense.  Then,  there  exists 

a  non-zero  Xq  €  X  such  that  <x^,x>  =  0  for  every  x  £  Dom(ja^). 
Since  by  assumption  R(_of-  XI)  =  X  for  some  X  >  0,  it  follows 
that  x q  =  {Si/-  XI)yQ  for  some  non-zero  yQ  £  Dom(j^)  .  There¬ 
fore, 

0  =  <Xq,x>  for  every  x  £  Dom(iaf) 

=  <XQ,yQ>,  in  particular, 

=  <(jaT-  xi)y0,y0> 

=  <^y0,y0>  -  X  ! y 0  |  2 

<  <^/yQ,yQ>, 

which  contradicts  the  dissipative  inequality  on  Si/.  □ 

Corollary  2.2.  Let  sZ:  (Dom(j^)  c  X)  X  be  linear,  X  a 
Hilbert  space.  If  a/  satisfies  a  dissipative  inequality 
<SiZx,x>  <  0  for  every  x  €  Dom(ja^)  and  if  R(_qO  =  X,  then 
Dom(_o^)  is  dense  in  X. 

Proof :  R(_qO  =  X  implies  0  £  p{aZ)  ,  the  resolvent  set  of 
Si/.  But  p {s/)  is  an  open  set,  and  thus  there  exists  a  X  >  0 
such  that  X  €  p(_aO  and  R{siZ  -  XI)  =  X.  The  above  theorem  now 
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yields  the  result.  o 

Finally,  we  will  obtain  the  desired  results  for  the  case 
6  >  0  by  taking  maximal  dissipative  extensions  of  We  con¬ 

sider  the  six  cases  corresponding  to  ,  k  =  1,2,3,  and  JaK  , 
k  =  1,2,3. 

Case  1.  Consider  the  case  corresponding  to  a  simply  sup- 

2  0 

ported  beam.  Consider  the  operator  in  Z  =  (a)  x  H 
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z1(x)  =  G1(x)  -  xG:(l)  -  (x4/6  -  x/6) Gj  (1)  where  G1(x)  = 


1 

a 


fx  rs, 


(S?  rs2fsi  42 

M  c  1  g (s) dsds. ds?ds _ .  Then  (zltz,)  £  H?  x  H,  = 
JnJnJn  iz  J.  i 


0  J0 


0  J0 


DomC^j)  mi  ^x(^)  *  (.oD4^) 


Lemma 


2.5.  is  densely  defined.  Hence  is  dense  in 


and  H4  is  dense  in  H^(a). 

A 

Proof :  Corollary  2.2  and  R(Ja^)  =  Z  and  dissipativeness 

of  yield  this  result  immediately. 

a  ^ 

Lemma  2.6.  srf  ^  is  skew- ad  j  oint :  Jaf*  =  -ja^. 

Proof :  First,  we  show  =>  Let  2  =  (zi»z2^  G 

Dom(j9^)  ,  and  find  all  w,  v  such  that  <j^^z,w>  =  <z,v>.  Thus 
w,v  €  Z  must  satisfy 


aD2z2D2w1 


0 


4 

aD  z-^w2 


n2  n2 

aD  z^D  Vj 


Z2V2  =  °’ 


where  w  =  (w^,w2),  v  =  (v1,v2).  Or,  integrating  by  parts, 

fl 


(2.6) 


|  z2^D4Wl-v2)  +  a[(Dz2)  (D2w1))  l  -  |^a(D2Zl)  (D2w2  +  D2V;i) 


0 


-  a[D3Zlw2)J]=  0, 


-V 

where  we  have  applied  the  conditions  on  z  €  Dom(_a?^)  ,  w  €  Z. 
Thus,  if 


v2  =  aD  w  4  2 

,  -  1  w  e  H4  x 

D2Vl  =  -D2w2 


and  if  w^'(0)  =  w^'fl)  =  0  and  w2(0)  =  w2(l)  =  0,  then  clearly 

A  A 

(2.6)  clearly  holds,  or  -J#  ^  is  adjoint  toJ3^. 


We  now  show 
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-safy  is  the  (maximal)  adjoint  operator;  i.e.,  ap'j  c  y : 

Let  g  €  Dom(j^)  ,  £  =  J^|g;  we  show  f  =  -_q?  y  g  and 
g  €  Dom(_Q^2).  We  have 


<z,f>  =  <z  ,j^*g >  =  <j£^z,g>  for  every  z  e  Dom(jaP'^) 


or,  with  £  =  (f^f^,  g  =  (g1,g2), 

r  loo  fl  fl 


0 


a(D2Z;i)  (D2f1)  + 


0 


Z2f2 


a 

a 

o 


)(D2z2)(D2g1) 
a(D4z1)g2  =  0. 


for  every  (z^,z2)  €  Dom(j9^)  .  Define  h^  (x)  =  I 
;1 


f2  (s)dsds1 


0 


0 


f2(s)dsds^.  Then  h^ (x)  satisfies  h”  =  f7  a.e., 


h  ( 0 )  =  h1(l)  =  0.  Integrate  the  above  by  parts  to  obtain 
fl  a(D4z1)f1  +  a[(D2z1)Df1]J  -  a  [(D^)  f  j  ]  J  +  |1(D2z2)h1 


0 


+  tz2(Dh1)]J  -  [(Dz2)h1]J  -  J*  a(D2z2)(D2gl)  +  a(D4Zl)g2=0. 

But,  f  6  Z  implies  f-^(O)  =  f^(l)  =  0;  also,  z  €  Dom(j^j)  im¬ 
plies  z2(0)  =  z2(l)  =  (D2z1D (0)  =  (D2z:) (1)  =  Zl(0)  =  z1(l)  =  0, 
and  so  this  becomes 


a|\D4z1)(f1  +  g2)  +  j\ D2z2)(hi  -  aD2g ^ )  =  0. 

Since  the  first  term  is  independent  of  the  second,  this  is  equi¬ 
valent  to  the  pair  of  equations 

{Vzptfl  *  g2)  -  o 
|\d2z2)  Ch1  -  aDZgl)  =  0. 


and 
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2  4  4 

Now,  for  every  v  €  ,  there  exists  z^  €  such  that  D  z^  = v 

2 

so  that  f-^  +  annihilates  all  v  €  H.,  and  so  +  g-,  =  0. 

2 

50  we  have  =  -f^  € 

0  2 

Likewise,  for  every  w  €  H  ,  there  exists  z7  €  such 

2  2  0 
that  D  =  w  and  so  (h^  -  aD  g^)  annihilates  all  w  €  H  . 

2  2  1 

This  implies  h^  -  aD  g^  =  0.  Hence,  D  =  ^  h^  = 

1  /  fx  fS  1  fl  fsi  \  3  \(  fx 

-  (  1  f t (s)dsds,  -  x  f  (s)dsds  ),  D  g.  =  -(  f?(s)ds  - 

a  \J0-,0  z  1  J0J0  z  1/  i  a \jQ  c 

f  1  1  s  i  \  41  3 

j  f2(s)dsds^j,  and  D  g^  =  —  f^  a.e.,  which  implies  D  g^ 

is  absolutely  continuous  and  D^g^  €  and  so  g^  €  H^.  Further¬ 
more,  (D2gl)(0)  =  i  h1(0)  =  0  and  (D2gl)(l)  =  I  h^l)  =  0. 

2  0 

This,  along  with  the  fact  that  g  €  Z  =  *  H  ,  implies 

51  €  H<. 

Thus,  g  =  (g1,g2)  €  x  hJ  =  dom(j/x)  and  (J2^*)g  =  f  = 

--^^g,  or  -Q?  *  c  -jtf  This  proves  the  desired  result,  since 

=>  -vr  ° 

Lemma  2.7.  s£ ^  is  a  closed,  maximal  dissipative  operator, 
and  is  the  generator  of  a  Cq  semi-group  of  contractions  on  Z. 

Proof :  That  ^  is  closed  follows  from 


*)*  =  ( --QT 1 )  *  =  -&*  =J y1. 


We  can  also  establish  maximal  dissipativeness  easily  by  noting 
that 

<Q^  *7-  ,Z>  =  <--Q^Z,Z>  =  0. 

This  implies  that  & ^  is  maximal  dissipative  by  a  theorem  of 
Krein  (25,  p.  873  .  That  *s  a  generator  then  follows  by 
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another  result  [25,  4.5,  p.  88]. 

4  2  1 

Remark .  We  have  shown  that  is  dense  in  H  n  Hq  m 

/v 

the  | * | 2  a  norm  by  showing  that-Q^  is  densely  defined  in  Z. 

In  this  case  the  result  can  also  be  argued  using  the  fact  that 

2  1  0 

A  with  Dom(A)  =  H  n  Hg  is  self-adjoint  in  H  and  the  re¬ 
sults  of  Goldstein  [18,  p.  863. 

Remark.  Since  srf ^  is  skew- ad j oint ,  we  could  also  have  used 
Stone's  theorem  [18,  p.  22]  to  argue  that  stf 2  generates  a  Cq 
group  on  Z.  However  we  emphasize  the  role  of  dissipativeness 
since  we  require  this  in  the  case  6  >  0. 

Case  2.  This  case  corresponds  to  a  beam  clamped  at  both 

2  0 

ends.  Consider  the  operator  in  Z  =  H2(a)  x  H  given  by 


-aD 


on  Dom(j^2)  =  H4  x  H2 


where 


and 


H 4  =  U  e  H4:  <K0)  =  ♦*(0)  =  <f>(D  =  <t>'  (l)  =  0} 
H2  =  {<t>  €  H2:  (K0)  «  <t>'(0)  =  <(>(1)  =  <t>'(l)  =  0} 

2  2 

H2 (a)  =  H2  equipped  with  the  inner  product 
<4>,'P>2,a  =  <aD24),D2^>0. 


Lemma  2,8.  <&2z,z>  0  for  every  z  €  Dom(j^2) 

Proof : 


2Z,Z>  =  <z2,z1>2fOL  *  <-uD  z1,z2>0 
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1 1 


2  2 

aD  z^D  + 


f  1 


-  otD  z^2 


aD2z2d2z1  -  a[(D3Zl)z2]Q  +  a  [  (D^)  (Dz2)  ]  J 

-  f  aD2z1I 
J  0  1 

A 

=  0  for  every  z  €  Dom(_Q^2). 


D2z. 


Lemma  2.9.  R(_q/2)  =  Z 


take  z2  =  f,  and 


Proof :  We  show  that  for  every  (f,g)  €  Z,  there  exists  a 

^ 2 1 » z 2 )  £  Dom(J^?)  such  that  a/ =  (g): 
z1(x)  =  GjCx)  -  ( 3x2  -  2x3)G1(1)  -  (x3-x2)GJ(l)  where 

1  g(s)dsds1ds2ds3 .  Then,  (z1,z2)  € 


G1(x) 


1 

a 


fx 

S3 

'S2 

0- 

0  • 

0  ■ 

rS] 

0 


H2  *  H2  '  Bomti,)  and  S?2Q‘)  -  (.^J  *  (*)• 

Lemma  2.10.  S& 2  is  skew-adjoint:  -2^*  =  '-^2' 

Proof :  First  we  show  z>  -3 2 '  Let  z  =  (z1»z2^  €  Dom(J^2) 


and  w  =  (Wj,w2)  €  Dom(_a/2)  .  Then  we  show  <&7z,v;>  =  <z,-J^2w>: 


<fl^2z,w>  -  <  z  ,  -_£^2  w>  = 


fl  2  2  f* 

aD  z^D  w^  + 


0 


0 


•a (D  z1)w2 


fl 


a(D2Zl) (-D2w2) 


z2(aD  w:) 


=  0. 

Thus,  -J*2  is  adjoint  to  _a?2.  Now  we  must  show  <=  -Jp’  . 

Let  g  €  Dom(-V^),  f  =  -2^*g;  we  show  f  =  -J&  g.  We  have 
<z,f>  =  <z,jrf*g>  =  <s^  2  z ,  g>  for  every  z  e  DomCj^),  or, 
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2  2 

aD  zxD  f:  + 


f  1 


2  2f  2 


1  2  2  r1  4 

D  z/g  +  a  D  z  g 

0  21  J0  12 


for  every  z  =  (z^,z2)  £  Dom(j^?).  Denoting  (x) 


rs. 


f7(s)dsds  +  ax  +  b ,  where  a,b  are  arbitrary,  a  £  R,  b  £  R, 


fX 


and  h7(s)  =  f?(s)ds  +  a  we  integrate  by  parts  to  obtain 
2  J  0  2 

-i  /•  1 

a(D4z1)f1  +  a[(D2Zl)  (Df^]  J  -  a  [  (D3Z;l)  f  ^  J  +  J  (D^)!^ 

+  [z2h2]J  -  [(Dz2)h1]J  -  af1(D2z2)  (D2gl)  +  f^t  (D^)  g2  =  0 . 


J  0 

2  0 

But,  f  €  Z  implies  €  H2  and  f2  €  H  ,  so  the  first  two 

A 

bracketed  terms  vanish;  likewise  z  £  Dom(ja^?)  implies  the  last 
two  terms  in  brackets  vanish.  Thus,  we  have 

rl  „  rl 


0 


a(D  z1) (fx  +  g2) 


(D2z2)(h1  -  aD2gl)  =  0. 


Since  these  two  terms  are  independent,  this  is  equivalent  to  the 
pair 

rl  4 

f  n  7  ^  r  4>  r, 

>2' 


and 


0 


(D  z. ) (f  +  g7)  =  0 


(D2z2)(h1  -  ctD2gl)  =  0 


for  every  z  €  Dom(j/2)  =  H2  *  H2 .  Thus  +  §2  =  0  imPlies 


®2 


fl  €  H2 ’ 


°2gl 


1  .  _  1 
a  hl  ”  a 


r  S 


0 


1  f 2 (s ) ds  +  ax  +  b,  D3g1  =  ^  j  f2(s)ds  +  a 


is  absolutely  continuous,  D4g^  *  ^  f2  £  H°  implies  g.^  £  H4  and 
g  €  Z  implies  €  H2,  so  that  g^  £  H4  n  H2  =  H4 .  Hence 


! 
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g  £  Dom(-2^2^  with  -£/*g  =  f  =  or  *  c  -j£'2 .  Thus 

jy*  =  &  2  • 

A 

Lemma  2,11.  _a^2  is  a  closed,  maximal  dissipative  operator 
and  is  the  generator  of  a  semi-group  of  contractions  on  Z. 

Proof :  Identical  to  the  proof  of  Lemma  6. 

4  2 

Remark .  The  statement  H7  is  dense  in  (in  the  usual 

2  OQ  QQ 

n  norm)  follows  directly  from  the  fact  that  c”  =  (<j>  €  C°°: 

2 

<P  has  compact  support  ir.  (0,1)}  is  dense  in  H?  and 

C0  C  H2  C  H2  (cf :  [30  *  p-  91] ) ‘ 

Case  3.  Finally  consider  the  case  corresponding  to  a  canti- 

■>  0 

levered  beam.  Consider  the  operator  in  Hj(ot)  *  h  given  by 


1  ’ 
0 


on  Dom(^3)  =  x 


where 

H ^  =  (d>  €  H4:  4(0)  =  ♦,(0)  =  4"(1)  =  4"’  (1)  =  0} 
H3  =  (4  €  H2:  4(0)  =  4>  ’  (0)  =  0} 

2  2 

and  H  ^  C a )  is  H3  equipped  with  the  inner  product  <*,*>7 

2  2  * 

«*D  • , D4 • >0 . 

Lemma  2.12.  <-^^z,z>  <_  0  for  every  z  €  Dom(J^3)  . 

Proof : 


<^3z,2>  =  <z2,Zl>2>a  +  <-aD4Zl,z2>( 


=  f1  a (D2z_) (D2z, )  +  f1  -  a(D4z,)z 
JO  L  1  j  0  1 
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g  €  Dom(^->)  with  -£*g  =  f  =  --^g,  or  *  c  -jp'2 .  Thus 

V*  =  jaP'7  • 

_  i- 

Lemma  2.11.  is  a  closed,  maximal  dissipative  operator 

and  is  the  generator  of  a  Cq  semi-group  of  contractions  on  Z. 

Proof :  Identical  to  the  proof  of  Lemma  6. 

4  7 

Remark .  The  statement  H,  is  dense  in  H"  (in  the  usual 

H“  norm]  follows  directly  from  the  fact  that  =  {0  6  C°°: 

2 

<P  has  compact  support  in  (0,1)}  is  dense  in  H7  and 

C0  C  H2  C  H2  (cf :  [30’  p-  91]) * 

Case  3.‘  Finally  consider  the  case  corresponding  to  a  canti- 

?  0 

levered  beam.  Consider  the  operator  in  H“(a)  *  H  given  by 


'  0 
-aD4 


1 

0 


on  Dom(i^3)  =  ^3  *  H2 


where 

H4  =  { 4>  €  H4:  $(0)  =  4>'(0)  =  <f>"(l)  =  <J>"  *  (1)  =  0} 

H2  =  { <J>  €  H2 :  4> (0)  =  4>-  (0)  =  0} 

2  2 

and  H,(a)  is  H,  equipped  with  the  inner  product  <•,•>., 

J  J  to  )  Ot 

<aD2. ,D2*>0. 

Lemma  2.12.  <-^^z,z>  <_  0  for  every  z  €  Dom(J^3). 

Proof : 

<J*3Z’Z>  =  <z2,z\>  2  ,a  +  <‘aI)4zl’22>0 

=  f1  a (D2z2) (D2z, )  +  f1  -  a(D4Zl)z7 

Jo  2  }  0 
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a(D2z2)  (D2Z;i)  -  ct  I  CD3Z;l  )  z2  )  J  ♦  a[(D2Z;L)  (Dz2)]  J 


-  f1  ct(D2z1 

1  n 


)(D  z2) 


0  for  every  z  €  Dom(S/ . 


Lemma  2,13.  -  Z. 

Proof :  We  show  that  for  every  (f,g)  £  Z,  there  exists  a 

(z-pZ-J  £  DomCJ^j)  such  that  take  z2  =  f,  and 

.2  3  2 2 

2j(x)  =  GjCx)  -  G^’(l)  -  (^-  -  ^_)G'"(1)  where 

■XfS  rs  fs  4  2 

G, (x)  =  ^  g  (s) dsds , ds 7ds , .  Then,  (z,,z7)  £  H,  x  H,  = 

1  '  O'  0  '  0  '  0  ,  1  2  21  i  i  $  ■> 

Dom(^3)  and  V3(.*)  -  (,oP4^)  -  (*)• 

/V  A  A 

Lemma  2.14.  ^  is  skew-adjoint:  ■&  *  = 

A  A  A 

Proof :  First  we  show  J3^*  ^  -~sa^».  bet  2  =  £  DomC^j) 

and  w  =  (w^,w2)  €  Dom(ja/3).  Then  we  show  <jS73z,w>  =  <z,-ja/jW>: 


<-3^z,w>  -  <z,--Q^3w> 


■  j'  *  j‘  -C”4M)"-2 

2  2 

a(D2z,) C-D2w?)  -  f  z-(aD4w .) 
n  1  2  J  n  2  1 


Thus,  --S^j  is  adjoint  to  Now  we  must  showr -Qh' *  c  -J3^  • 

Let  g  €  Dom(Q^|),  f  =  a^g;  we  show  f  =  -Ja^g.  W'e  have 

A  A  A 

<z,f>  =  <z,ja/?g>  =  <-Q/,z,g>  for  every  z  €  Dom(J3^_),  or, 


C  “CdShD2^)  -  {'  *2*2  -  cpy^HD2^)  *  o|V4z1 


)g7  =  0 


for  every  z  =  ( z  ,  z  2 )  £  Dom(s^j).  Denote  by  h(x)  the  function 


■X  fS  (-1  fl  ,s 

h(x)  =  f?(s)dsds1  -  x  f7(d)ds  -  f?(s)dsds, 

■Vo  -  1  Jo  Vo 

1 

+  f 2 (s) ds 

0  z 

? 

satisfying  D“h  =  f  7  ,  h ( 1 )  =  h'(l)  =  0.  Integrate  by  parts  to  ob¬ 
tain 

^  a(D4z1)f1  +  a[(D2z1)(Df1)]J  -  a[(D3z1)f1]J  +  j*(D2z2)h 

+  [z2(Dh)]J  -  [(Dz  2 )  h  ]  -  aj\D2z2)  (D2gl)  +  (D4Z;i)g2  =  0. 

But  f  =  (f^.f-J  €  Z  implies  f-^CO)  =  f^(0);  likewise, 
z  €  Dom(ji'3)  implies  that  ZjfO)  =  z|(0)  =  2  ”  Cl)  =  zj"(l)  =  0 
and  z2(0)  =  ^(O)  =  0.  Thus  the  above  becomes 

jo  a(D4Zl)(f1+g2)  +  |\d2  22)  (h  -  aD2§1)  =  0. 

Since  the  first  term  is  independent  of  the  second,  this  is  equi¬ 
valent  to  the  pair 

{Vzpcq  -  i2)  ■  o 

and 

jo(D2z2)(h  -  <xD2gl)  =  0. 

2 

Thus,  +  §2  =  0  which  implies  g2  =  -f1  €  .  Furthermore, 

h  -  aD2g1  =  0  implies  D2g;L  =  ^  h,  D3gl  =  i  (j  f2(s)ds  - 

|  f2(s)ds^,  and  D4g^  =  ^  f2 ,  so  that  D3g^  is  absolutely  con¬ 
tinuous  and  D4gl  €  H°  which  implies  g^  e  H4 .  Also  (D2g^)(l)  = 
^  h(l)  =  0  and  (D3g^ ) (1)  =  ~  h'(l)  =  0.  This,  together  with  the 
fact  that  g  =  (gl,g2)  €  Z  =  H2  x  H° ,  implies  that  g  e  H4  x  H2  = 
DomCj^j),  with  (J^*)g  *  f  *  -J^g,  or  -^3  <=  "-^3*  This  Proves 
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A  A 

the  desired  result,  since  we  already  have  stf  *  r>  □ 

Lemma  2.15.  S&  ^  is  a  closed,  maximal  dissipative  operator 
and  is  the  generator  of  a  Cq  semigroup  of  contractions  on  Z. 

Proof :  Identical  to  the  proof  of  Lemma  6. 

Case  4.  Now  consider  the  general  case  corresponding  to  a 
simply  supported  beam  with  internal  damping.  Define  the  operator 
in  H2(ct)  x  h°  =  Z  by 


on  Domfja^)  =  H4  x  H4 

where  a€R,  a>0,  6€R,  6>0, 

H2  =  (4>  £  H2:  0(0)  -  0(1)  =  0} 

H4  =  {0  €  H4:  0(0)  *  0" ( 0)  =  0(1)  =  0"(1)  =  0} 

2  2  2  2 
and  H^a)  is  with  the  inner  product  <-,->2  a  ~  <ctD  *  »D  •  >0 

By  the  results  of  Section  1.2,  H4  is  dense  in  H°  and  by  Lemma  2.5 

4  2  , 

above,  is  dense  in  H^(a);  thus  is  densely  defined. 

Lemma  2.16.  is  dissipative. 

Proof:  Let  z  =  (z^,z2)  £  Domf-Q^).  Then, 


<  Jafz ,  z  > 


<z2,zl>2,a  +  <_aD  Z1,Z2>0  +  <'6D  Z2’Z2>0 
a|^(D2z2) (D2z1)  +  |^(-aD4z1)z2  +  j  (-6D4z2)z2 

-a([(D3z1)z2]  l  -  [(D2Zl)(Dz2)]J)  -  «|\d2z?)  2 

-  6([(D3z2)z2]  l  -  [(D2z2)(Dz2)]J)  =  -6  |\d2z2) 


<  0  for  6  >  0. 
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2  2 

Remark .  Note  that  we  required  (D  z^)  (0)  =  (D  z7)(l)  =  0 
to  establish  the  dissipative  inequality,  and  this  is  one  reason 
for  including  these  conditions  in  Dom(-9^). 

-9^  itself  does  not  generate  a  CQ  semigroup  on  Z,  since 
R(_9^  -  XI)  ^  Z.  However,  does  possess  a  maximal  dissipative 

extension  9^  which  is  closed  (cf:  [25,  p.  87]).  In  the  case 

A 

6  =  0,  we  have  seen  that  the  maximal  dissipative  extension  is  sd  ^ . 
Since  Dom(-9^)  is  dense  in  Z,  the  following  dense  inclusions 
hold : 

Dom(  9^)  e  Dom(ja^)  c  Z, 

and  the  maximal  dissipativeness  of  Sd  ^  is  sufficient  to  ensure 
that  it  is  the  generator  of  a  CQ  semigroup  of  contractions  on  Z. 

2  0 

Case  5.  Consider  next  the  operator  in  H^fa)  *  H  given  by 


^2  " 


-aD4  -6D4 


on  Dorn  ( -9^)  =  H4  x  H4  ,  a,  6  >  0 


This  corresponds  to  the  beam  clamped  at  both  ends  as  in  Case  2, 
but  with  internal  damping  included. 


Lemma  2.17.  sd^  is  dissipative. 


Proof : 


See  Lemma 
<  -9/  z  ,  z  >  = 


2.16. 
fl 


-  6 


0 


The  same  integration  by  parts  yields 

(D2z2)2  £  0 


for  every  z  =  ( z  ^  ,  z  2 )  e  DomC-Q^)-  D 

Since  R(-9^2  ’  /  Z,  we  will  again  need  to  take  a  maximal 

dissipative  extension  -9^2  where 

Dorn  (  sd 2 )  c  Dorn  (  9^ )  c  Z 


are  dense  inclusions. 


Lemma  2.18  and  a  theorem  of  Krein  [25,  4.3, 


p.  87]  are  sufficient  to  guarantee  the  existence  of-Q^*  1°  t^ie 
case  6  =  0,  we  have  seen  that  =  & 2‘  another  standard  re¬ 
sult  ([25,  4.4,  p.  87])  £^2  is  closed,  and  generates  a  CQ 

semigroup  of  contractions  on  Z. 


Case  6.  For  the  final  case  we  wish  to  consider,  we  take  the 
2  0 

operator  in  Z  =  x  h  defined  by 


^ . 


•  aD 


-5D 


on  Dom(i^3)  =  H4  x  H3 


where 


and 


H4  «  £  H4:  4>(0)  =  4>'(0)  =  <J>"(1)  -  ♦,,,(1)  =  0), 

H2  =  (f  £  H2:  <|>(0)  =  <p(l)  =  0}, 

2  2 

HjCa)  =  equipped  with  the  inner  product 


<  • 


> 


0* 


Lemma  2.18.  is  dissipative. 


Proof : 
<-Q^2,Z> 


See  Lemma  2.16.  The  same  integration  by  parts  yields 

f  1  2  2 

=  -6j  (D  Z2)  0  for  every  z  =  (z^^)  €  Dom(Ja^). 


0^3  has  a  maximal  dissipative  extension  by  the  theorem  of 

Krein  [25,  4.3,  p.  87],  and  since  Dom(  is  dense  in  Z, 

Domf-Q^)  c  DomCJa^j)  c  Z  where  all  inclusions  are  dense.  Thus 
[25,  p.  88]  -^3  generates  a  Cq  semigroup  of  contractions  on  Z. 
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We  shall  require  certain  dense  subsets  of  smooth  functions 
in  Dom(j^^)  to  apply  Proposition  1.5  (Trotter-Kato  theorem). 
To  construct  these  dense  subsets,  we  will  use  eigenfunctions  of 
B4  on  Dom(D4)  =  H4. 

2 

We  first  note  that  if  u  £  H^,  k  =  1,2,  or  3,  then 
max  |u(x) |  <  K|u| 2  » 

and 

|u(xx)  -  u(x2)|  <  K  /x:-x2  | u | 2 . 

These  inequalities  follow  from  an  application  of  the  Cauchy- 

Schwartz  inequality  and  from  the  fact  that  for  u  e  ,  lu^  < 

K | u | 2  for  some  constant  K  (this  is  true  because  |*|  and 

1  2 

| * | 2  are  equivalent  to  the  usual  H  and  H  norms  respect- 

2  ? 
ively  on  Hj"  and  |u|  £  |u|  9).  In  particular,  for  u  e  IC , 

k  Hi  k 

fx2 

x1,x2  £  [0,1],  u(x2)  =  u  (x^ )  +  u'(£)d£  implies 

jxl 

|u(x2)  -uCXj)  |2  <  |  2  u *  (  C) d C|2 

Jxl 

2 

<  | x 2 - |  |u'|q,  by  the  Schwartz 

inequal i ty , 

=  lx2  "xi  I  lu  li  I  1(2  lx2  X1  I  lu  if  * 

Then  we  obtain  from  this  (with  x^  =  0)  that 

|u(x2)  |  <  K/x^"  |u  1 2  for  every  x2  £  [0,1] 

which  implies  the  first  inequality.  Thus  an  application  of  the 

Ascoli  Theorem  (see  [38,  pp.  249-250])  implies  that  from  every 

2 

sequence  (un(x)}  which  is  bounded  in  the  norm  of  H^,  it  is 

possible  to  select  a  uniformly  convergent  subsequence.  Addi- 

2 

tionally,  for  any  two  elements  u,v  €  H^, 
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i <u , v> 0 I  <  maxju|  max | v (  £ K2 | u | 2 | v | 2 . 

4  2 

Now  let  Bj,,  k  =  1,2,  or  3  be  the  operator  D  in  with 

Dom(B^)  =  H^.  Bj^  exists  (see  proofs  of  Lemmas  2.4,  2.9,  and 
2.13) ,  and 


- 1  2 
<Bj,  u,v>2  =  <u,v>0  for  every  u,v  €  . 

B^  is  clearly  symmetric  and  furthermore  it  is  a  compact  opera- 

2  2 
tor  on  H^:  given  a  bounded  sequence  (un(x)}  in  ,  we  can 

select  a  uniformly  convergent  subsequence  (vn(x)},  and 


- 1  2 
and  thus  [38,  p.  206]  B^  is  a  compact  operator  from  X  = 

into  X. 

For  the  moment,  fix  k  =  1,2,  or  3.  Since  B’1  is  symmetric 

and  compact,  standard  results  [45,  p.  343]  imply  that  B^1  has  a 

complete  orthonormal  set  of  eigenfunctions  \p  and  associated 

00  ^ 

eigenvalues  X  ,  and  f  =  £  <f,\p  >-\p  for  every  f  €  Dom(Bv)  =  h5. 

^  n  c  n  k  k 

We  have  B,  ip.  =  n  f  ,  where  p  =  1  /  X  . 

k  .7  n  Ti  n 

0  _ 

We  can  construct  a  dense  set  in  H  by  taking  <f>n  =  /p^  . 

4 

Then,  since  for  every  f  £  H^, 


<f’V2  =  <I)2£>d2Vo  =  <f’I)4*n>0  =  <£*Wo  =  <£’Vo  /V 


m 


we  have  |£  <f >4'n>2^n  ’  fl2 
m  £ 


0  implies  that 
l  <£ »<fn>o<fn  "  f|2  — >  °>  which  implies,  since  the  norm 


is 
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m 

stronger  than  the  |*|0-norm,  that  |£  <f,$n>4»n  '  fl0  — *  0  for 

1 

every  f  €  H? .  As  a  consequence,  the  set  span({<J>  })  is  dense 
x  n 

in  H°,  since  is  |*|0-dense  in  H°.  Also,  B^<|>n  =  for 

n  1,2,...  . 

The  smoothness  of  the  functions  <J>n  (we  shall  only  re- 

6  2  2 
quire  H°  smoothness)  follows  from  the  fact  that  B]<^n  = 

2  8 

implies  that  €  Dom(B^)  c  H  ,  and  similarily  for  4>n . 

2  0 

Thus  for  our  dense  subsets  of  smooth  functions  m  Z  =  * H  , 

we  are  led  to  define 


f 

f  o 

Tu  = 

J 

.  0  . 

>Y2,j  = 

, 

l 

Let  3 N  «  spaniy,  .  U  spanly,  .  }^.  ,  and  U  ^N.  Clearly, 

1  >  J  J  ”  1  ^  9  J  J~  x  1 

by  the  above,  ^  is  a  dense  subset  Z,  and  c  Dom(-Q^J  <= 

Dom(ja^)  c:  Z. 


Remark.  The  set  (4>^ }  associated  with  B,1  and  the  set  & 

j  K 

constructed  from  the  will  clearly  be  a  different  set  of  func¬ 

tions  for  each  k.  Since  their  properties  are  the  same,  we  de¬ 
note  these  by  the  generic  symbols  {<)>.}  and  where  it  is 

•  2 

understood  by  context  that  these  represent  the  functions  in 
2  0 

and  x  H  respectively. 

We  require  one  further  property  of  @  to  apply  Proposition 
1.5  (Trotter-Kato) .  We  must  show 


(2.7) 


(ja^.  -  AqI)£0  =  Z  for  some  Aq  >  0. 


First  observe  that  -  AqI  is  invertible  for  Aq  >  0,  since 
generates  a  contraction  semigroup  on  Z,  and  further  that 


when  restricted  to  Dom(ja^).  Then  since  ££  c  Dom(-a^), 
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cvk  -  xon^  =  c  ^k  -  \Qi)& . 

In  fact,  we  show  that  (Ja^k  '  Agl)^  =  5^  and  the  result  (2.7) 
follows  from  the  denseness  of  Q>  in  Z. 


Lemma  2.19.  (o/  k  -  A^I)^  for  some  A^  >  0. 

Proof :  a)  (jai^  '  A^I)^c^  is  trivial, 

b)  (fl/k  ~  AqI)^  =  ^:  Since  (_ci^k  ‘  AQI)  *  exists,  it 
suffices  to  show  that  for  a  typical  basis  element  this  holds.  We 
do  the  typical  basis  elements  y1  .  and  y~  ■  separately.  We 
show  that 


i)  there  exists  a  z  =  (z^,z2)  €  such  that 


f-^k  '  X0I)z  •  Y2,  j  *  (*j)- 
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there 

exists  a 

(^1 

-  = 

iff 

(J^l  -  *0 

iff 

{ _iV  * 

1  'aD  2i 

iff 

-(a+5AQ)D 

'2  0  2 

,2  _  ,„4. 


with  z ^  =  a^.,  this  can  be  solved  for  (2^,2.,)  €  <^  (using 


D 


iff  -(o  +6A0)y^a  -  A“a  =  6y .  +  a( 


or 


a  *  - 


tu]  *  xo 


<“  *SAo)lij  *  'o 


Since  all  quantities  in  the  denominator  are  positive,  this  solu¬ 
tion  is  possible  (a  finite)  for  any  AQ  >  0.  □ 

Thus  for  each  R  =  1,2,  or  3  we  have  a  set  @  such  that 


<=■  Dom(^k)  «=  Dom(.a^k)  <=  Z, 
where  all  inclusions  are  dense,  and 

(j^k  -  aqi =  (Ja^k  -  A0l)S?  -  z 

for  any  Aq  >  0. 

We  can  sununari2e  the  previous  results.  The  operators  -^k  > 
k  =  1,2,3  defined  in  Hk(ct)  x  H°  with  Dom(-^k)  =  ^k  *  Hk 
have  a  maximal  dissipative  extension  J^k  on  Dom(-Q^k)  where 
Dom(.Q^k)  c  Dom(J^k)  c  Z  densely.  Furthermore,  for  each  operator 
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k  =  1,2,3,  there  exists  a  set  52  consisting  of  smooth  func¬ 
tions  such  that  52=  Z  and  (jai^  '  XI)52  =  Z. 

We  now  consider  the  approximate  problem.  Define  for  each 
k,  k  =  1,2,3,  the  spaces  Z  and  Z^  by 

ZN  =  Sjj(AN)  x  Sjj(AN),  Z  =  Z(q)  =  Z(a)  =  h£(«)  x  H° 

N 

(recall  a  =  q  ),  and  let  P‘  (q)  be  the  orthogonal  projection 
N  N 

P  (q) :  Z  — >  Z  .  The  approximate  problem  corresponds  to  solu¬ 
tions  of 

i\t)  =  ^N(q)zN(t)  +  FN(q,t) 
zN(0)  =  zjjfq) 

where  ^N(qj  =  PN(q)  J^(q)  P%)  ,  FN(q,t)  =  PK(q)F(q,t),  and 

Zo’(q)  =  PN(q)z0(q)  • 

We  give  a  concrete  realization  of  these  approximating  equa¬ 
tions  at  the  end  of  this  section,  along  with  numerical  results  in 
Section  4.  Our  first  goal  is  to  prove  convergence  of  qN  — >  q 
in  the  context  of  Chapter  1.  That  is,  we  prove  that  solutions  of 
the  approximate  identification  problem  (IDN)  converge  to  solutions 
of  the  full  identification  problem  (IDA),  using  Proposition  1.4. 

We  first  prove  the  approximating  subspaces  converge  in  the 
appropriate  sense. 

N 

Lemma  2.20.  Assuming  (HQ),  the  projections  P  (q)  converge 
strongly  to  the  identity  I  in  Z  as  N  “,  for  k  =  1  or  2. 

Proof :  Let  z  =  (z-^.z^  6  52  where  52  is  the  dense  subset 
in  Z  defined  in  Lemma  2.19  for  k  =  1,2  respectively.  Then 
PN(q)z  =  (zj,z2)  =  (PjZ1,P2z2) ,  where  P ^  is  the  projection  of 
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the  first  coordinate  onto  S^(A^)  in  the  l*L  norm  and  P ^ 

k  2 ,  a  2 

is  the  projection  of  the  second  coordinate  in  the  I 'ip  norm* 

Then , 

|PNCq)z-z|2  =  I P321_z1 I  2 ,  +  lP2Z2'22lo 

-  ^  k3,0  lD  2I0)  +  ^  k2,0*'N')  zU')’ 

where  we  have  applied  Lemmas  1.18  and  1.19,  which  implies 
|PN(q)z-z|  0  as  N  00  for  z  €  2.  But  since  is  dense  in 

Z  and  the  projections  P^(q)  are  bounded  operators,  it  follows 
that  |P^(q)z-z|  — >  0  for  any  z  £  Z.  □ 

We  may  now  state  and  prove  the  main  theorem  of  this  section. 

Theorem  2.21.  Let  hold.  Then  the  semigroups  T(t;q) 

and  TN(t;q)  generated  by  J^^(q)  and  respectively, 

k  =  1  or  2,  satisfy  ||T(t;q)||<  1  and  ||  TN(  t ;  q)  ||  <  1.  More¬ 
over,  for  any  sequence  {qN}  converging  to  q*  in  Q,  we  have 

I  T^  (t ;  qNj  z  -  T(t;q*)z|  — >  0  uniformly  on  [  0  ,  T] 

for  each  z  €  2 . 

Proof:  The  bound  |J  T  ( t ;  q)  {}  <_  1  follows  because  gen¬ 

erates  a  CQ  semigroup  of  contractions  on  Z.  Also  by  the  re¬ 
marks  in  Chapter  1,  Section  1,  it  follows  that  ||T  it;q)|]  <  1.  To 
establish  the  convergence  results,  we  apply  the  Trotter-Kato 
theorem  (Proposition  1.5).  Let  Z(q*),.#N  =  ZN(qN)  ;  also 

take  -2/=  -^k(q*)  and  _a^N  =  J^k(qN)  .  Let  ttN:  Z(q*)  — >  Z(qN) 
be  the  cononical  isomorphism  between  Z(q*)  and  Z(qN) .  Then 
q  — >  q*  implies  |  7tNz  |  — >  |zj,  verifying  hypothesis  i)  of 
Trotter-Kato.  Thus  it  remains  to  verify  iii) .  We  have  already 


Ki-* 
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defined  for  each  k  a  set  such  that  (jj^Cq*)  -  >1)5^  is 

dense  in  To  establish  convergence,  we  see  that  for  each 

z  =  ( z j » z 2 )  (suppressing  the  notation  irK) 

(2.8)  |^N(qN)z  -  J^(q*)z|  =  |  PN(qN)  J^(qN)  PN(qN)  -  J^(q*)z| 

<  |  (_a/(qN)  -  ^(q*))PN(qN)z| 

+  |  V(q*)(PN(qN)z-z)  |  +  |(PN(qK)-I)^(q*)z|  . 


We  bound  each  of  these  terms  separately. 

The  second  term  can  be  written  explicitly  as 


Mq*)  (PN(qN)  Z-Z)  = 


■q*D 


-q^D 


(PN(qN)z-z) 


ZN  .  _ 

2  2  z2 


-q*D4(z^'Z1)  -  q*D4(z*2-z2) 

where  P^(q^)  =  (z^.z^),  and  Zj  =  P^z.^,  z^  =  ?^z2>  ''here  P^ ,  P^ 
are  the  projection  operators  of  Lemmas  1.19  and  1.18.  Thus, 


|j^(q*)(PK(qN)z-z) I2 


N 

z-,- 


'4(z?- 

,_N 


2  ~  2 2  ^  2  ,  q*  +  lqlD  (zi'21)+<?2d4Cz2-z2^  0 


-  lZ2"Z2^2,q*  +  {q*lD  Iq 

+  I  q  *  I  I  D4  ( z  ^  -  z  2 )  |Q}‘ 


-  ql  K2,2(^N')  Z2  ^  0  +  ^ql  <3,2^N')  21*0 

!->  2  i  ,.6. 


+  lq2l  *2,4^  I  D  z2  I  0^ 


by  an  application  of  Lemmas  1.18  and  1.19.  Thus, 
I  ^(q*)  (PN(qN)z-z)  |  ■*>  0  as  N  -*•». 
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\’ 

For  the  last  term  in  (2.8),  we  have  P‘  (q)  -*•  I  strongly 
on  Z.  Convergence  of  the  first  term  follows  once  we  write  out 
explicitly 


(J^(qN)  -  ^(q*))PX(qN')z  = 


and  note  that  D4z^  -*•  D4z^  and 
q^  -*■  q*.  Thus  Trotter-Kato  (Prop 
semigroups  |T^(t;q^)z  -  T(t;q*)z 
each  z  €  Z.  □ 


0 

-(q^-q*)D4Zj  -  (q2-q*)D4Z2 

D4z‘^"  -*■  D4Z2  in  and 

.  1.5)  yields  convergence 
|  0  uniformly  on  ( 0  , T ] 


J 

of  the 
for 


Remark .  Note  the  role  that  the  dense  subsets  Q)  played  in 

the  proof,  in  addition  to  possessing  the  required  smoothness  to 

apply  Lemmas  1.18  and  1.19.  Since  c  Dom(ja^)  c  Dom(j^)  ,  we 

were  able  to  use  the  form  of  -0^  explicitly  in  the  proof  since 

S&  ■  restricted  to  Dom(-Q^.)  is  itself.  Thus  the  Trotter- 

l  v  i  i 

Kato  approach  does  not  require  that  we  know  explicitly.  □ 

A  direct  application  of  Proposition  1.4  along  with  the  com¬ 
pactness  assumption  (HQ)  yields  convergence  of  a  subsequence  of 
— N 

the  solutions  q  of  the  approximate  identification  problem  to  a 
solution  q  of  the  original  identification  problem  (IDA) . 

The  case  where  y  >  0  in  (2.2)  follows  easily  from  pertur¬ 
bation  theory.  In  this  case  in  (2.5)  becomes 


\  - 


-ctD 


-6D4-y 


Consider  =  -3^,  +  B,  where 


■aD 


-6D 


as  before, 
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and  ^ ^ts  maxi-mal  dissipative  extension,  and  B  = 

is  a  bounded  perturbation  of  since  |Bz|  <_  Y I  2  2  ^  0  ^or 

every  z  =  ( z  ^ ,  z  2 )  €  Z.  So  standard  perturbation  results  [18, 
p.  33]  imply  that  generates  a  CQ  semigroup  on  Z.  It  is 

easily  verified  that  (  -  AI)£2  =  £2.  The  proof  of  Theorem  2.21 

remains  essentially  unchanged  by  the  inclusion  of  this  new  term. 

We  now  turn  to  a  concrete  realization  of  the  approximate 

N 

identification  problem.  Let  B^  be  the  quintic  B-spline  satis¬ 
fying  the  appropriate  boundary  conditions,  as  defined  in  (1.16)- 
(1.18),  and  recall  that  {B^}^_o  are  a  basis  for  S^(A^).  As  a 
basis  for  ZN  =  sjj(AN)  *  Sjj(AN),  then  we  take  where 

eN  ,  |  0  i  1  i  N 

1  1  N*1  5.  i  5 


Our  approximating  equation 


zK(t)  =  -0/  z  (t)  +  PNF  (t) 

is  then  defined  by  _q/N  =  PN£*TPN  (with  F  =  (0,f)T),  where  PN 
is  the  orthogonal  projection  of  Z  onto  ZN.  To  obtain  a  reali¬ 
zation  of  this  method,  we  seek  a  z^(t)  €  Z^  satisfying 

<zN(t)  -  ^NzN(t)  -  PNF(t),£>  =  0  for  all  c  €  ZN. 

Since  zN(t)  and  ?  €  ZN,  they  have  representations  z^(t)  = 

£  Wj(t)Bj  and  C  =  The  condition  above  thus  reduces  to 


<1  *j(t)8j  -  I  wjjtfp'j1  '  =  0 


for  i  =  0 , 1 , . . . , 2N+1 ,  which  in  turn  is  equivalent  to  the  vector 
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system  of  ordinary  differential  equations 


(2.9) 


QNwN(t)  =  KNwN(t)  +  RNF(t) 


QN'wN(0) 


N 

R1  'F 


where  wN(t)  =  (wjjct)  , .  . .  twJN+1  (t)  )  ,  Q^-  =  < e i  »  B j  >  *  Kij  = 

(R^F)j  =  <B^,F>,  and  V  =  C 4> » 'l')  is  the  vector  of  ini¬ 
tial  conditions  in  (2.4).  For  the  examples  under  consideration 
(as  well  as  others  presented  below) ,  the  matrix  structure  in  the 
approximating  equations  facilitate  computations.  In  particular, 
one  finds 


qn  = 


>N 


0 


„ZnN 


•W 


where  (Qj)^  =  a<D"B”  ,D"B^>0 ,  (Qp^  =  <B‘.\B:?> 


KN  - 


K 


N 


,N 


kN 

K3 
,4„N  „N 


and 


with  =  Q‘p  (Kp.j  =  <  -aD  BV,Bj  >q  =  -  (Qp  ij  ,  and 
(K^ij  -  <  -6D4B^7 ,Bj  >  -  <YB^,Bj!>0  =  -  |(Qj,)ij  -  yCQ^)  i-j  »  where  one 
uses  integration  by  parts  and  the  boundary  conditions  to  establish 
these  identities.  Thus  the  equation  in  (2.9)  reduces  to 


(2.10) 

where 


wN(t)  =  GNwN(t)  +  FN (t) 


,N 


-(Q?)-1Q? 


I 


Y  I 


and 
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fn  =  cqn)_1rnf  = 


0 

(Q2 ) "lfN 


Kith  fN  =  <f ,  B  j  >Q  ,  j  =  0,1 . N. 


N, 


The  approximate  displacement  y  (t,x),  used  in  the  cost  func¬ 
tional  J,  is  given  by 

y  (t ,  * )  =  z , (t)  =  l  w. (t)  B .  . 

1  i=0  1  1 

N  N 

The  matrices  and  Q7  have  a  banded  structure  and  can  be  used 

efficiently  in  solving  (2.10).  More  will  be  said  about  this  in 
Chapter  4. 


Section  3.  An  Approximation  Using  Cubic  Splines 

In  Section  2,  we  solved  the  approximate  identification  prob¬ 
lem  (IDA)  using  state  approximations  based  upon  quintic  splines. 

In  this  section,  we  lower  the  smoothness  requirement  for  the 
basis  elements  by  rewriting  (2.2)  as  an  abstract  equation  which 
permits  the  use  of  cubic  splines.  To  do  this,  we  consider  here 
only  the  case  w'here  6  and  y  are  zero  (no  damping),  with  bound¬ 
ary  conditions  of  type  1,  corresponding  to  a  simply  supported 
beam.  Then  we  may  rewrite  (2.2)  as  an  abstract  equation  in 
Z  =  x  x  of  the  form 


(2.11) 


t(t)  =  Mq)z(t)  +  F(q,t),  t  >  0, 
2(0)  =  zQ 


where  z(t)  =  (z2 (t) ,z2 (t) ,z3(t) )  =  (y (t , • ) ,yt (t , • ) ,yxx (t ,  • ) ) , 
2 q  =  (4>,^.4>xx)  ,  and 
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</(q) 


0  1  0 
0  0  -  /a  D2 

0  /a  D2  0 


on  Dom(a/(q))  =  x  n  hJ  x  H2  (1  Hq.  Here  E  a. 

We  first  wish  to  prove  that  ja/(q)  generates  a  CQ  semigroup 
on  Z.  The  following  lemmas  will  do  this,  starting  with  dissi¬ 
pativeness  of  .£/( q)  . 


Lemma  2.22. 
large . 

Proof :  Let 

<  C  ja/(q)  -<ui)v,v> 


(j^(q)  -  (joI)  is  dissipative  for  oj  sufficiently 


v  =  (v1,v2,v3)  €  Dom(  a^(q) ) .  Then, 


= 

<V2>V0 

i 

V 

+ 

D2 

v3’V  + 

2 

</a  D  v2,v3> 

,2 

i  1 2 

i  1 2 

-  w|v1 

lo 

-  v2 1 0 

-  “1 v3( 0 

<v2>Vo 

1 

i  2 

i  i  2 

i  i  2 

-  WlVj 

l0 

-  u)|v2l0 

-  v3< 0 

1  . 

2  i  i 

2  ,|2 

< 

lV2l0lVli 

o  '  “1 

Vll 

0  '  v2> 

o  -  “lv5lo 

i  i  2  i 

i  2 

1  v2  1  0 

vi  1 0 

„  i  2  | 

i  2 

< 

2  +  ~ 

~~2 — 

“1 

vi  1  o  "  “1 

V2I0  -  WIV3I 

= 

( j  ~  u) 1 v 

i  2 

iU  + 

<7 

-  «)|v2|J 

-  a>|v3lo 

< 

0  for  a) 

>  1/2 

□ 

To  establish  that  Ja^(q)  is  a  generator  of  a  Cq  semigroup,  we 
will  show  that  ja^(q)  is  closed  and  R(ja^(q)  -  XI)  =  Z.  For 
simplicity,  we  first  partition  the  operator  _o/(q)  as 
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’01  0 
Mq)  =  o 
[  0 

where 

0  -/cT  D2 

/a  D2  0 

2  12  1 

on  Dom(B(q))  =  H  n  Hq  x  H  fl  Hg,  and  derive  some  results  for 
B(q)  which  will  be  used  to  show  that  q)  is  a  generator  of  a 
CQ  semigroup. 

Lemma  2.23.  B(q)  is  a  dissipative,  skew-adjoint  operator 
and  generates  a  Cg  semigroup  of  contractions  on  Y  =  H°  x  h°. 

Proof:  Dissipativeness  follows  easily  from  an  integration 

by  parts.  Let  v  =  (v^v^  €  Dom(B(q));  then,  <B(q)v,v>  = 

<-/a  D2v2v1>0  +  </a  D2v1,v2>q  =  0. 

It  is  easily  verified  that  -B(q)  is  adjoint  to  B(q) 

(B*  o  -B) : 

<B(q)v,w>  -  <V2-B(q)w>  =  <-/a  D2v2,w1>0  +  </a  D^.w  >Q 

-  <Vl,^  D2w2>q  -  <v2,-/o  D2w1>q 

=  0  for  all  w  =  (w1,w2)  €  Dom(B(q)). 

We  show  B*  c  -B.  Let  g  =  (g1,g2)  €  Dom(B(q)*),  f  =  (fj,^)  = 
B(q)*g.  We  show  f  =  -B(q)g  and  g  £  Dom(B(q)).  We  have 

<v,f>  =  <v,B(q)*g>  =  <B (q) v , g>  for  every  v  €  Dom(B(q)), 

or, 

j„  vlfl  *  j0  v2f2  *  I"  -  jo  <W>2Vl)g2  *  0 


65 


rx 


rs. 


for  every  v  6  Dom(B(q)).  Let  h.,  (x)  = 

i  1  >  n; 0 

f1  fs  i  fx  fs  i 

x  f, (s)dsds.  and  h. (x)  s  f0(s)dsds.  - 

JqJO  1  1  1  J  o J  0  L  1 


fj (s) dsds^ 


0J  0 
1 1 1  s. 


'fSl 

f?(s)dsds,.  Then  integrating  the  terms  in  the  above  equa- 
0J0  z  1 

tion  by  parts  yields 


j0(I)2vl)hl  *  j0tI)2v 2)h2  *  l0  /“(D2v2)Sl  ■  f0  ^<d2v1>S2  "  °- 

2  1 

where  we  have  applied  the  conditions  v^»v2  ^  H  n  Hq  and 

h^(0)  =  h^(l)  =  0,  for  i  =  1,2.  Thus  we  obtain  the  pair  of  equa¬ 

tions 

f1  2 

o(DZv1)(h1  -  /Zg2)  =  0 

and 

r1  2 

(DZv2)  (h2  +  /agj.)  -  0. 

Thus  g_  =  —  h.  ,  Dg_  =  —  Dh  ,  D2g  =  ~  f  €  H°  implies  g?  €  H2, 

✓a  -1  /a  /a  1  ^ 

and  g2(0)  =  —  h.  (0)  =  0,  g7(l)  =  —  h,  (1)  =  0  implies 

2  1  2  1  2 
g2  e  H  0  HU’  Similarly,  we  obtain  g1  £  H  fl  Hq  and  D  gj  = 

-  ~  f2*  and  30  (fl’f2)T  =  B<^*g  =  D2g2>  ‘  ^D2gi)T  =  -B(q)g. 

Therefore,  B(q)  is  skew-adjoint.  Since  <B(q)*v,v>  = 

<-B(q)v,v>  =  0,  B(q)  is  maximal  dissipative  [25,  p.  87],  and  so  it 

generates  a  CQ  semigroup  of  contractions  on  Y  [25,  p.  88].  o 

We  may  now  state  and  prove  the  result  for 


Lemma  2,24.  _Q^(q)  is  the  generator  of  a  CQ  semigroup 
T(t;q)  on  Z  =  H°  x  h°  x  H°. 


Proof :  Let  z  =  (z1,z2,z3)  €  Dom(  Ja/(q) )  .  Since  B(q)  gen¬ 
erates  a  semigroup  on  H°  x  H°,  we  have  [25  ,  p.  87]  R(B(q) -  I)  = 
HB  x  for  any  X  >  0  and  B(q)  is  a  closed  operator.  But 
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since 


(-af(q)  -  XI)2  = 


z2  ’  *Z1 


.  CB(q>-»I)(^) 


R(B(q)  -  AI)  =  x  implies  that  given  f9,f_  €  H®,  there 

exists  z2»z3  e  Dom(B(q))  such  that  (B(q)  -  AI)^2^  =  (f^)’ 

,0 


then,  given  €  Hv 


j(z2  '  e  and  z  solves 


(J*(q)  -  AI)z  =  f,  £  =  which  implies  R(_af(q)-AI)  =  Z. 

Also  the  fact  that  B(q)  is  a  closed  operator  yields  immediately 
that  ja^(q)  is  a  closed  operator,  and  these  two  conditions  along 
with  the  dissipativeness  of  ja^(q)  -  <ul  are  sufficient  [25,  pp .  87- 
88]  to  ensure  that  ja/(q)  generates  a  semigroup  on  Z .  □ 

We  again  require  a  dense  subset  of  smooth  functions 
Dom(a^(q)).  Note  that  for  <f>j  (x)  =  /I  sin(jTrx),  forms 


a  complete  orthonormal  set  for  H°.  Define 


♦j 


f 

0 

'  0  ' 

*J  ! 

[  0  J 

Y2,j  = 

♦j 

,  and 

Y3,j  ■ 

0 

.  Then  define 

0 

.  . 

>N 


«  span{YlJ}JT=1  U  span{Y2J>j  =  1  U  span  {y^  >  .  }*  =  1 , 

00 

2>-  U  .  then  Q)  =  Z.  For  A  >  0,  it  is  easily  verified  that 

N=  1  - 

(AI  =  Q>,  and  so  (AI  -  srf{q))Q>  =  Z. 


i  N 


tN 


and 


For  the  approximate  identification  problem,  define 


tN 


.3,aN 


.3 ,  A  N 


■  3 ,  aN 


N 


N. 


,N 


In  this  case,  for  z  =  (z^,z2,Zj)  €  Z, 


Z  =  Sq(A  )  x  Sq(A  )  x  Sq(A  ),  and  let  P  be  the  orthogonal  pro¬ 
jection  P 
N  N  N  N  N 

P  z  =  (PqZj ,PqZ2 ,poz3)  where  Pp  is  the  orthogonal  projection  of 
H°  onto  S03(AN).  The  approximate  identification  problem  then 
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3Svm 


corresponds  to  solutions  of 

(2.12) 


z  (t)  = 


N(q)zN(t)  +  PNF(q,t),  t  >  0 


zN(0)  =  zjj 

where  ja^N(q)  =  PI'',Q^'(q)PN  and  z^  =  PNZq. 

We  defer  obtaining  a  concrete  realization  of  this  problem  to 
the  end  of  this  section,  after  we  have  proved  convergence  of  a 
subsequence  of  the  solutions  q‘  of  the  approximate  identifica¬ 
tion  problem  converge  to  a  solution  q  of  the  original  identifi¬ 
cation  problem. 

Lemma  2.25.  Assuming  (HQ),  the  projections  PN  converge 
strongly  to  the  identity  I  in  Z. 

Proof:  Let  (z1,z2,z3)  €  where  &  is  the  dense  subset 

@cz  Dom(_£/(q) )  c:  Z  as  defined  above.  Then,  P^z  =  (PqZj, 

N  N  \T  \7  XT  \T 

P0Z2,P0Z3-*  ”  (zi»z2»z3^>  where  PQ  is  the  |  •  |  Q -pro  j  ection  of  Z 
onto  Sq(A^).  Then, 


„N  ,2  r  | dN  ,2 

P  z’zl  *  il-L  lP0Zi'2i  1  0 


£  .^[*0,0^  I  c4zi  I  0  ]2  »  by  Lemma  1-5- 


Therefore  P 


N 


l  on  2$,  and  the  boundedness  of  the  projection 
_  N 

operator  implies,  since  Q  =  Z,  that  P  — >  I  strongly  on  Z.  o 
We  now  state  and  prove  the  convergence  result  for  cubic  spline 


approximations . 


Theorem  2.26.  Let  (HQ)  hold.  Then  the  semigroups  T(t;q) 
and  TN(t;q)  generated  by  ja^(q)  and  ja^(q)  satisfy 
1 1 T  ( t ;  q )  |(  <  Mewt  and  ||TN(t;q)||  <_  Mewt: .  Moreover,  for  any 


» 
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sequence 


{qN}  converging  to  q*  in  Q,  we  have 
|TN(t;q^)z  -  T(t;q*)z|  — >  0  uniformly  on 


[0  ,T] 


for  each  z  6  Z. 

Proof :  The  bound  ||T(t;q)||  <  Meut  follows  because 

generates  a  Cp  semigroup  on  Z  [31,  p.  10].  Since 

q)z,z>  =  <  Ja^(q)  P^z  ,  P^z>  =  w|P‘^z|  <_  cu  |  z  |  ,  by  the  dissipative¬ 
ness  of  -  oil,  for  a)  sufficiently  large,  we  obtain 
II  TN(t ;  q)  ||  <  Meut. 

We  refer  again  to  the  Trotter-Kato  theorem  (Proposition  1.5) 
to  establish  convergence.  Let  30  =  Z,  30^  =  z\  and  take 
to  be  the  identity.  All  that  remains  to  be  verified  is  condi¬ 
tion  (iii)  of  Trotter-Kato. 

Taking  &  as  above,  we  see  that  for  each  z  =  (z1,z2>z.5)  6  , 

|j/N(qN)z--^(q*)z  I  *  |PN-o^(qN)PKz- j/(q*)z  | 

(2.13)  i  I  (J^(qN)-Mq*))PNz| 

+  |^(q*)(PNz-z)|  +  |  (PN-I)ja^(q*)z|  . 


We  bound  each  of  the  terms  above.  The  second  term  may  be  written 
explicitly  as 


J*( q*)  (PNz-z) 


0 

0 

0 


1  0 
0  -  /a  D2 

✓o  D2 


,„N 

(P  z 


-/5’D2(z^-z3) 

✓aD2(z^-z2) 


z) 


0 
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where  q*  =  a  and  z ^  =  P^z^,  i  =  1,2,3.  Thus, 

i  Mq*)(PNz-z)  |2  =  I  z  2  "  2  2  I  0  +  «Id2(z3-z3)Io  +  «  |  D2  (zj-z^  | 2, 
which  implies,  by  Lemma  1.15, 

|ja/(q*)(P  z-z)(  <  [<0,0(N)  Z2^o]  +  a[<0,2(-N-*  Z3^o] 

+  a[K0,2  I  D  z2  I  o]  » 


and  so  |  ja/(q*)  (P  z -z)  |  -*■  0  as  N  -*•  <*>. 

For  the  last  term  of  (2.13),  we  have 


I  strongly  on 


Z.  Convergence  of  the  first  term  follows  explicitly  from  the 
form  of  J3^(q)  : 


(^(qVa'(q'>))Pt'’z 


-(/q?  -  ^)D2z3N  . 

-(/qf  -  •'  q$Wz* 


Convergence  follows  since  D  z, 

* 

in  H  ,  and  by  hypothesis  q^  - 


D2z2  and  D2Zj 


This  theorem,  together  with  Proposition  1.4,  is  sufficient 
to  ensure  that 

lim|qN-q*j  =  0  implies  lim| zN(t ;qN) -z (t ;q*) |  =  0. 

N-*®  N-*-® 


We  now  turn  to  a  concrete  realization  of  the  approximate  prob- 

N 

lem  using  cubic  splines.  Let  {C^}  be  the  cubic  splines  defined 
in  (1.15)  which  span  Sq(AN).  The  approximation  subspace  ZN  is 
Sq(AN)  *  Sq(AN)  x  Sq(AN).  We  follow  the  procedure  outlined  at  the 
end  of  Section  2  where  we  obtained  a  concrete  realization  using 


quintic  splines. 


We  form  projections  :  2  -*•  2^  and  seek  approximate  solu- 

N 

tions  of  the  form  z^(t)  =  £  w^(t)C^  where  is  given 

i=0  1  1  1 

by 


i  =  0  ,  . . . ,N, 


*i-< 


(C^,0,0)T 

(0’Ci-(N+l) ’°)T  1  =  N+1»*-*>2N+1> 

(0,0,C^.(2N+2))T  i  =  2N+2 » . . . , 3N+2 . 


We  are  thus  led  to  a  system  of  3N+3  differential  equations  for 
N 

the  w^  (compare  with  the  2N+2  system  in  the  quintic  formula¬ 
tion.') 


(2.14) 

where 


.  N  ,  _N  N.  A  _N n 
w  (t)  =  G  w  (t)  +  F  (t) 


,N 


0 

0 

0 


0 

rflN, -l.N 
(A2)  a2 


„,.N\  -l.N 
-a(A2)  A1 


with 


and 


^i,j  -  <ci*cj>o*  ^ij  -  <K>»ch 

rnN.  j-  r  .N  .N 
(Q  )  =  diag(A2,A2 ,A2) 


(QN) _1rnf 


(RNF)i  - 


0 

<f  cN 

r,Li-(N+l)  2 


0  <  i  <  N,  2N+2  <  i  <  3N+2 
N+l  <  i  <  2N+1 . 


Finally,  we  present  an  approximate  method  using  cubic  splines 
in  the  case  where  structural  damping  is  included.  While  the  con¬ 
vergence  proof  falls  outside  the  framework  above,  we  present 
numerical  results  in  Section  4  which  demonstrate  the  effectiveness 
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of  this  method.  We  next  consider  (2.2)  with  Y  =  0,  a  >  0,  6  >  0 
and  boundary  conditions  of  type  1  (simply  supported)  .  The  initial 
conditions  are  the  same  as  above.  Equation  (2.2)  is  rewritten  as 

y  =  v 

6D2u  +  v  =  -ctD2u 
n2 

u  =  D  v 

so  that  for  z  =  (y,v,u)T  =  (y.y^y^)1  in  Z  =  H°  x  H°  *  H°,  we 
obtain  the  abstract  equation 

Tz(t)  =  j*z(t) 

(2.15)  T 

z(0)  =  (♦,*,♦"). 

Here  is  given  by 

’  0  1  0 

S?  =  0  0  -aD2 

0  D2  0 

on  Dom(_aO=  (H2  n  hJ)  x  (H2  n  hJ)  x  (H2  n  hJ),  and 

1  0  0 

r  =  0  1  6  D2 

.0  0  1 

Note  that  only  if  <5  =  0  (no  structural  damping)  does  this  reduce 
to  the  form  (2.11)  and  fall  within  the  framework  of  the  theory  out¬ 
lined  for  (2.11)  and  (2.12).  Nonetheless,  the  general  ideas  dis¬ 
cussed  above  lead  to  efficient  computational  schemes  for  (2.11). 
Taking  ZN  =  Sq(AN)  x  Sq(AN)  x  Sq(AN)  and  PN  the  orthogonal 

projection  of  Z  onto  ZN,  we  define  =  P^Ja^PN  and 

rN  „NrnN  . 

T  =  P  TP  and  use 
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(2.16)  rNzN(t)  *  JtfNzN(t) 

as  the  approximating  equation  for  (2.14). 

Section  4.  Numerical  Results 

We  provide  a  few  examples  of  numerical  experiments  illustra¬ 
ting  the  methods  presented  in  this  chapter.  Many  examples  have 
already  been  reported  in  [  6  ]  which  compare  the  method  based  upon 
quintic  splines  to  the  method  based  upon  cubic  splines  for  the 
Euler -Bernoull i  equation  with  structural  damping.  Included  were 
examples  with  time-dependent  boundary  conditions  which  were  trans¬ 
formed  to  homogeneous  boundary  conditions.  We  refer  the  inter¬ 
ested  reader  to  those  examples  and  do  not  repeat  them  here. 

We  provide  some  new  examples  not  included  in  [  6 ] .  The 
first  example  concerns  a  simply  supported  beam  with  both  struc¬ 
tural  and  viscous  damping.  The  others  illustrate  the  convergence 
properties  for  the  important  case  of  the  cantilever  beam.  All 
of  these  examples  use  the  quintic  spline  approximation. 

The  numerical  experiments  consisted  of  taking  as  data  the 
values  of  a  solution  of  a  model  equation  of  the  form  (2.1)  whose 
parameters  were  known  and  then  seeking  a  solution  q^  of  the 
approximate  identification  problem  for  different  values  of  N. 

The  "data"  {y^},  i  =  l,...,r;  j  =  1 , . . .  ,  Jl  were  generated  by 
solving  (2.1)  approximately  by  a  Crank-Nicolson  scheme  and  using 
y  —  =  y(t^,Xj)  as  the  observations  or  input  to  the  approximate 
identification  packages.  For  the  experiments  listed  below,  we 
used  l  =  3,  r  =  10  and  generated  the  "observations"  at 
=  j/4>  3  =  1,...,3  for  times  ti  =  i/10,  i  =  1,...,10. 


The  optimization  algorithm  (see  Chapter  4)  requires  an  ini¬ 


tial  guess  for  the  parameters;  these  are  referred  to  as  the 
"start  up"  values  in  the  tables.  The  values  to  which  the  optimi- 

_ v 

ration  algorithm  converged  for  a  fixed  N  are  denoted  by  q1  , 

— N  _v 

and  J(q')  denotes  the  cost  functional  evaluated  at  q'  . 

Example  2.1.  We  consider  an  example  where  we  hav'e  included 
both  structural  and  viscous  damping  in  the  model.  We  consider 
the  model  equation 


=  -q,y  -  q,v  „  -  q_v-  +  ex  sin  t 

•  tt  ll'  xxxx  ~  xxxxt  M37t 

yCO.x)  =  yt(o,x)  =  n 

y(t,0)  =  >’xx (t » 0)  =  y(t , l)  -  >'xx(t,l)  =  o 

Table  2.1  summarizes  the  results  for  the  parameter  estimation 
problem  for  this  equation  using  the  quintic  spline  approximation. 


N 

N 

qj_ 

Table  2.1 

N 

N 

J(qN) 

2 

. 999938 

.009875 

.032420 

.2001  x  10~5 

3 

.999940 

.009974 

.022713 

.8231  x  io*7 

4 

.999936 

.010004 

.019860 

.3324  x  10'7 

TRUE  VALUE 

1.0 

.01 

.02 

START  UP 

.65 

.005 

.005 

Exampl 

e  2.2. 

Next  we  consider  an  example  for 

a  cantilevered 

beam.  With 

4>(x) 

=  sin  ax  -  sinh  ax 

+  K(cosh  ax  - 

cos  ax)  where 

K  =  (sin  a 

+  sinh 

a)/(cos  a  +  cosh  a) 

and  a  -  7 . 

8547  (corres- 

ponding  to  the  third  mode  for  an  undamped  cantilevered  beam)  ,  we 


consider  the  equation 


^tt  ^l^xxxx  ^2^xxxxt 
y  (0  ,x)  =  <(»  Cx)  ,  yt(0,x)  =  0 

y(t,o)  =  y x  C  t , o )  =  yxx(t,i)  =  yxxx(t,i)  =  o. 
The  numerical  results  are  summarized  in  Table  2.2. 


N 

ii 

Table  2.2 

s 

J(qN) 

2 

.  3821 

.0067 

.025  x  10'1 

3 

.4960 

.  0098 

.735  x  10'4 

4 

.4996 

.010020 

.116  x  io'6 

6 

.4997 

.009995 

.768  x  io'7 

TRUE  VALUE 

.5 

.01 

START  UP 

.35 

.005 

Example  2.3.  A  second  example  for  a  cantilevered  beam  is 
provided  for  the  model  equation 

ytt  •  -Vxxxx  •  V'xxxxt  *  <>10  «P(-20(l-x)]exp[-q9t),  t  >  0 
y(o,x)  =  yt(o,x)  =  o 

y(t,o)  =  yx(t,o)  =  yxx(t,i)  =  yxxx(t,D  =  o. 


Table  2.3  summarizes  the  results  for  the  estimation  problem  corres¬ 
ponding  to  this  equation  using  the  quintic  spline  approximation. 


Table  2.3 


N  qj 

^2 

N 

4)  J^N) 

2  .497863 

.009923 

2.02101 

10.0458  .110  x  io‘5 

3  .498527 

.009995 

2.01284 

10.0255  .693  x  io'7 

4  .499284 

.009996 

2.00602 

10.0143  .449  x  10'8 

8  .499775 

.  009982 

2.00183 

10.0075  .356  x  10'9 

TRUE  VALUE  .5 

.01 

2.0 

10. 

START  UP  .35 

.005 

1.5 

o 

CO 

When  q^  and  q^g 

were  treated  as  known 

values  and  optimization 

was  performed  for 

and  q2 

only,  the 

following  results  were 

obtained . 

Table  2. 

4 

N 

n  N 
^1 

N 

2i 

J(qN) 

2 

. 500394 

.009445 

.146  x  10'5 

3 

. 500045 

.009763 

.141  x  10'6 

4 

.500006 

.009817 

.340  x  10*7 

TRUE  VALUE 

.5 

.01 

START  UP 

.35 

.005 

We  note  that 

in  all  of  the 

examples 

considered  for  the  Euler 

Bernoulli  equation  rapid  convergence  of  to  the  true  valued 

occurred . 

Fc  the  optimization  algorithm  (see  Chapter  4) ,  we  specified 
small  convergence  tolerances  to  obtain  best  possible  results. 

For  a  typical  problem,  convergence  was  obtained  in  3-5  iterations 
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of  the  Levenberg-Marquardt  algorithm. 

A  variable  step/variable  order  method  (DGEAR)  was  used  to 
integrate  the  approximating  system  of  ordinary  differential  equa¬ 
tions.  For  example  2.2  with  a  requested  local  error  tolerance 
of  1.  x  10*6,  DGEAR  took  211  steps  (last  stepsize  was  .018)  and 
used  fifth  order  methods  to  integrate  the  system  when  stiff 
methods  were  used.  By  comparison,  when  non-stiff  methods  were 
chosen,  it  took  5766  step  (last  stepsize  was  .213  *  10  to 
achieve  the  same  local  error  and  it  primarily  used  order  2 
methods,  indicating  a  moderate  degree  of  stiffness  due  to  the 
damping  term. 

As  was  the  case  with  the  second  order  equations  considered 
in  [  7 ] ,  the  approximations  for  linear  constant  coeffici¬ 

ent  problems  were  extremely  accurate,  even  for  small  N.  It  is 
expected  that  the  power  of  these  approximations  will  be  more  im¬ 
portant  for  other  problems  to  be  considered  in  the  future, 
particularly  the  case  of  spatially  varying  coefficients. 


vrii'  ~ , 


CHAPTER  3.  APPLICATION  TO  THE  TIMOSHENKO  EQUATION 


Section  1.  The  Timoshenko  Equation 


The  second  example  of  the  application  of  these  techniques 
that  we  wish  to  consider  is  the  identification  problem  associated 
with  the  Timoshenko  equations  for  transverse  vibrations  of  beams. 
The  Timoshenko  theory  extends  the  Euler-Bernoulli  theory  by  tak¬ 
ing  into  account  the  effects  of  rotary  intertia  and  shear  dis¬ 
tortion.  These  effects  play  a  significant  role  when  the  depth 
of  the  beam  is  large  when  compared  to  its  span  and  when  high  fre¬ 
quency  oscillations  must  be  considered. 

While  the  Timoshenko  formulation  may  be  written  as  a  single 
fourth  order  partial  differential  equation,  it  is  easier  to 
handle  boundary  conditions  when  it  is  written  as  a  system  of  two 
partial  differential  equations  in  the  transverse  displacement 
y(t,x)  and  angle  of  rotation  ^(t»x)  of  the  beam  cross-section 
from  its  original  vertical  position  [14,  p.  300]: 


(3.1) 


ytt  ■  at>,xx  -  V  * 

■  c\x  *  btyx  -  *>• 


t  >  o,  x  e  [0,i] 


where  a  =  k'AG/m,  b  =  Aa/',  c  =  El/m,  with  A  =  cross  sectional 
area  of  beam,  E  =  Young's  modulus,  G  =  shear  modulus,  I  =  moment 
of  inertia,  and  k'  =  shear  coefficient  (cf.  [15]),  and  where 
f(t,x;q)  is  the  "load"  or  applied  force. 

The  initial  conditions  are 


t 
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y  (o 

,x)  = 

(x) 

yt(o 

,x)  = 

(x) 

iKO 

,x)  = 

^0 

(x) 

*t(0 

,x)  = 

♦l 

(x)  , 

some 

of 

the 

common 

.  boundary 

conditions  are 

i) 

y(t 

,0) 

=  y(t. 

1)  = 

*(t. 

0) 

=  *(t,l)  =  0 

ii) 

y(t 

,0) 

=  y(t. 

1)  = 

2i(t 

3xl 

,0) 

-  H1*'15  ■  ° 

iii) 

y(t 

,0) 

=  <Kt, 

0)  = 

*£rt 

,D 

=  iKt.l) 


These  correspond  to  the  boundary  conditions  for  a  beam  fixed  at 
both  ends  (i) ,  a  simple  supported  beam  (ii)  ,  and  a  cantilevered 
beam  (iii)  ,  comparable  to  the  boundary  conditions  of  type  k 
for  the  Euler-Bernoulli  equation  (cf.  [13,  p.  97]). 

For  the  ease  of  exposition,  we  shall  limit  our  discussion 
to  boundary  conditions  (i)  ,  although  the  theory  is  generally 
applicable  to  all  three. 


Section  2.  An  Approximation  Using  Cubic  Splines 

Following  again  the  approach  outlined  in  Chapter  1,  we  write 
(3.1)  as  an  abstract  equation 


(3.2) 


i(t)  =  ^q)z(t)  +  F(q,t) 
2(0)  =  zQ(q) 


in  a  Hilbert  space  Z,  where  q^  =  a,  q2  =  b,  q3  =  c,  and 
q4,...,qp  are  parameters  in  the  load  term  F  and  initial  func¬ 


tion  Zq. 


■ju.  'A •  -  VuiVaa 
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Again,  there  are  many  choices  of  (3.2)  and  the  space  Z 
which  lead  to  well-posed  problems.  The  forms  that  we  choose  for 
(3.2)  lead  to  specific  state  approximations.  We  shall  discuss 
two  such  possibilities  and  discuss  the  approximate  identification 
problem  of  the  first  in  detail. 

The  first  such  choice  of  the  abstract  equation  will  be 
(3.2)  with  z(t)  =  (z1(t) ,z2(t) ,z3(t) ,z4(t))T  =  (y (t , • ) ,yt (t , • ) , 

^  (t ,  • )  ,ipt  (t ,  •  )  )  ,  in  Z  =  V(a)  *  H°  *  V(c)  x  H^1,  and  with 


0 

1 

0 

0 

aD2 

0 

-aD 

0 

(3.3) 

Mq)  =  o 

0 

0 

1 

bD 

0 

cD2  -b 

0 

on  Dom(j^)  =  H2  n  hJ  *  hJ  *  H2  n  hJ  *  hJ,  where  V(a)  =  hJ(cx) 
equipped  with  the  inner  product  <u,v>vj-a^  s  <u,v>1  ^  =  <aDu,Dv>Q. 
We  denote  the  corresponding  norm  on  V(a)  as  | • | .  ,  and  take 

F(q,t)  =  (0,f(t,-;q),0,0)T. 

Theorem  3.1.  The  operator  s/(q)  defined  by  (3.3)  with 

^1  =  a*  ^2  =  q3  =  c  *s  tbe  generator  semigroup 

T ( t ;  q )  on  Z  satisfying  ||  T ( t ;  q)  ||  <_  ewt  for  some  a>  >  0  . 


Proof :  Define  the  operator  on  Z  by 


where 


'0  1  0  0 

<(q)  -  aD2  0  0  0  -  f  ®l  ° 

0  0  0  1  0  B2 

0  0  CD2  0 

0  1  2  11 
B,  =  ,  on  Dom(B. )  »  H  n  H«  x  H«  and 

1  aD^  0  1  u  u 

0  1  2  11 
7  on  Dom(B7)  =  H  D  Hn  *  Hn.  B.  and  B-.  are 

CD  0  u  u  i  l 
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the  wave  operators  considered  in  [7,  p.  29],  where  they  were 
shown  to  be  generators  of  Cq  semigroups  of  contractions  in 
V(q^)  x  H°  for  i  =  1  and  i  =  3  respectively.  Thus  9^  is 
the  generator  of  a  Cq  semigroup  of  contractions  Tp(t;q)  in 
Z,  9^  is  maximal  dissipative,  and  ||  TQ  (t ;  q)  ||  <_  1. 

Now  we  apply  perturbations  results  to  obtain  the  conclusion 
of  the  theorem.  We  first  show  that  9^  is  a  bounded  perturbation 
of  the  operator  9^.  Note  that  9f=9^  +  sif  where 


Thus  srf  is  a  bounded  perturbation  of  9^,  and  we  can  apply  the 
perturbation  results  of  [31,  p.  80]  to  infer  that  srf  generates 
a  Cq  semigroup  on  Z,  satisfying 


II  T(t ;  q)  ||  <  e  ||  9/  ||  t  = 


e 
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Corollary  3.2.  The  operator  {srf  -  u)I)  is  dissipative  for 
all  a)  sufficiently  large. 

Proof :  Let  z  €  Dom(ja^)  .  Then 

<S&z,z>  =  <(-&0  +&)z,z> 

=  <  q/qZ,z>  +  <-Q^z,z> 

<_  <J*z,z>,  since  ^  is  dissipative, 

I  |-*z|  I  z  I 

<_  |  z  |  by  the  boundedness  of  srf  . 

Thus  to  I )  2  ,  z  >  <_  0  for  all  oj  K^. 

The  proof  of  the  above  theorem  is  interesting  in  that  it 
states  that  the  Timoshenko  equations  can  be  viewed  as  a  pertur¬ 
bation  of  wo  simple  wave  equations.  Because  of  this  fact,  we 
were  able  to  apply  the  results  of  (  7 ]  where  the  wave  equation 
involving  an  operator  of  the  form  B-^  was  treated. 

Corollary  3.3.  Let  =  Dom(j^)  n  (Hm  x  Hm  x  Hm  x  Hm)  , 
m  ^  3.  Then  -  Aot)j^  =  Z  for  some  AQ  sufficiently  large. 

Proof :  Since  s&  generates  a  Cq  semigroup  on  Z, 

R(-Q^-  AqI)  =  Z  for  Aq  sufficiently  large.  So  for  such  AQ 
in  tho  resolvent  set  p(_»0,  Srf -  AQI  is  invertible. 

Consider  Y  c  Z.  Y  =  (Hm  x  Hm'2  x  Hm  x  Hm*2)  n  Z.  Clearly 
Si 4-  AqI  is  invertible  on  Y,  so  take  f  =  (f  ^  , f  2  , f ^  ,  f 4)  e  Y. 

We  show  there  exists  a  z  6  such  that 

(3.4)  (j/-  AQI) z  =  f. 

We  know  by  the  invertibility  of  of-  AQI  that  there  exists  a 


z  £  Dom(j*0  =  Hq  n  H  x  Hq  x  Hq  fl  \r  x  hJ  such  that  (3.4)  holds 
We  need  only  show  z  €  But, 

"X02l  +  2 2 

-A0z2  +  aD2z1  -  aDz3 
'Vs  +  2  4 

2 

-A-z.  +  bDz.  +  cD  z,  -  z, 

0  4  1  3  3 

So,  (j/-  AqI)z  =  f  if  and  only  if 

(i)  -AoZl  +  z2  -  f2 
2 

(ii)  ‘AqZ2  +  aD  z1  -  aDz3  =  ^7 

(iii)  -  AqZ^  +  z4  =  f3 

(iv)  -AQz4  +  bDZj^  +  cD2z 3  -  cz3  =  f 4 . 

First  consider  m  =  3.  Then,  f  €  Y  =  (H3  *  H1  x  h3  *  H1)  n  Z 
implies  f2  €  H1,  and  z  €  Dom(ja^)  implies  z^  £  HX, 
z3  €  H2  n  Hq  (hence  Dz3  €  H1)  .  So  g  =  aDz3  +  ^q22  +  ^2  e  ^ 
and  (ii)  implies  D  z -^  =  g/a  €  H  and  €  H  .  Then  (i)  im¬ 
plies  z2  £  H3.  Similarly,  f  6  Y  implies  f^  £  and 
z  £  Dom(_o/)  implies  z^  £  H2  D  H^,  z^  £  Hq  and  so  g  =  (*0Z4  ' 

bDz^  +  cz  +  f^)  £  H  and  so  D  z3  =  g/c  £  H  ,  which  implies 

3  7 

z3  £  H  .  And  then  (iii)  implies  z^  £  H  . 

Next  consider  m  =  4.  We  have  z  £  (H3  x  H3  x  H3  x  H3)  fl 

Dom(j^)  from  above,  and  we  assume  f  £  (H4  x  H2  x  H4  x  H2)  n  Z. 

2  2 

Repeating  the  arguments  in  the  last  paragraph  yields  D  £  H 
which  implies  z ^  £  H4  and  z2  £  H4  (from  (i)).  Likewise, 

D2z3  €  H2  implies  z3  €  H4  and  z2  £  H4 . 


{srf-  AQI) z  = 
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A  simple  induction  argument  yields  the  result  for  arbitrary 
m.  [Note:  we  actually  only  use  m  =  4]  .  So,  {£/-  AQI)^=  Y, 
and  ( &  -  X Q I =  Y  =  Z  since  Hra  is  dense  in  Hm  , 
j  -  1 , . . . ,m.  o 

We  may  now  discuss  the  state  approximations  which  we  use  to 

solv  the  approximate  identification  problem.  We  take 
N  3N  3N  3N  3N 

Z  =  S Q C A  )  x  Sp C A  )  x  Sq )  x  Sp(Aw),  where  we  recall  that 
S^(A^)  is  the  set  of  cubic  splines  s(x)  with  knots  on 
satisfying  s(0)  =  s(l)  =  0.  Define  P^(ot)  the  projection 
P^(a) :  V(a)  — >  Sq(A^)  in  the  | • | ^  a  norm,  and  Pq  to  be 
the  projection  P^:  H^1  — >Sq(A^)  in  the  norm.  Then  ;t 

follows  immediately  from  Lemmas  1.15  and  1.16  that  with  P^ 
defined  as  the  projection  P^:  Z  — >  Z^  given  by  P^z  = 

^•P1Z1  ,P0Z2  ,P1Z3  ,P0Z4^  »  tliat  (see  also  ^  7  >  P-  33H  the  following 
holds . 

N 

Lemma  3.4.  P  — =►  I  strongly  on  Z. 

With  and  P^1  defined  as  above,  we  solve  the  approximate 

identification  problem  associated  with 


(3.4) 


zN(t)  =  _0^N (q)  zN (t)  +  FN(q,t),  t  >  0 

Nrn-,  _  dN 
z  (0)  -  P  zQ 


in 


where  again  we  take 
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a  Cg-semigroup  T^(t;q)  =  ^  0n  Z.  By  Theorem  3.1,  J^(q) 

is  the  generator  of  a  Cq  semigroup  satisfying  ||T(t;q)||  <_  ewt, 
so  that  ja/(q)  -  wl  is  the  generator  of  a  semigroup  of  contrac¬ 
tions  and  is  maximal  dissipative  [25,  p.  90].  Thus  by  the  re¬ 
marks  in  Chapter  1,  Section  2,  #^(q)  -  wl  generates  a  Cq  semi 
group  of  contractions,  or  _o^(q)  generates  a  Cq  semigroup 
TN(t;q)  satisfying  ||  TN(t ;  q)  ||  <  eut . 

We  can  now  state  and  prove  the  convergence  result  that  the 
solutions  of  the  approximate  identification  problem  using  state 
approximations  (3.4)  converge  to  the  solution  of  the  identifica¬ 
tion  problem  corresponding  to  (3.2  -  3.3). 

Theorem  3.5.  Assume  (HQ)  holds.  Then  ja^(q)  and  _a/^(q) 
are  generators  of  CQ  semigroups  satisfying  ||  T ( t ;  q)  ||  <_  ewt  and 
HT^(f,q)||  5  eajt  for  some  m  >  0 .  Moreover,  for  any  sequence 
{ qN)  converging  to  q*  in  Q,  we  have 

|T^(t;q^)z  -  T(t;q*)z|  — >  0  uniformly  on  [0,T] 
for  every  z  e  Z. 

Proof :  We  have  already  proved  the  first  assertion.  We  in¬ 
voke  the  Trotter-Kato  "’heorem  once  more  to  prove  convergence. 

Let  3  =  Z(q*),j0N  =■-  ZN(qN),  and  let  ttN:  Z(q*)  — ^  Z(qN)  be 

the  canonical  isomorphism  between  Z(q*)  and  Z(qN) .  Then 
N  i  N  i 

q  — >  q*  implies  |  tt  z]  — >  z  ,  so  that  the  first  hypothesis 
of  Trotter-Kato  is  satisfied.  Thus,  it  remains  to  verify  hypothe 
sis  (iii)  .  We  have  already  defined  the  set  such  that 
(ja/(q*)  -  AI)£2  is  dense  in  <2. 
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To  prove  that  the  last  hypothesis  holds,  let  z  = 

XT 

Zj»z^)  €  Again,  we  shall  suppress  the  notation  it  ,  and  ob¬ 


tain 


|j^N(qN)z  -  q*)z|  =  |PN(qN)  J3f(qN)PN(qN)  -  J^(q*)z| 

£  |(j^(qN)  -  J*(q*))PN(qN)zi 

+  |^(q*)(PN(qN)z-z)  | 

+  1  (PN(qN)-I)  Qfi q * )  z  |  • 


(3.5) 


We  estimate  each  of  these  terms  separately.  The  second  term  in 
(3.5)  can  be  written  explicitly  as 


o/(q*)  (PN(qN)z-z)  = 


0 

l 

0 

0 

q*D2 

0 

-q*D 

0 

0 

0 

0 

1 

q*D 

0 

q§D2-q* 

0 

(PN(qN)z-z) 


N 

z2  ~z2 

qlD2(zl'Zl)  '  qlD(z3'z3^ 

Z4  ‘  z4 

q*D(z1-z1)  q*D  (Zj-Zj)  -  <\2(Z5~Z3^ 


where  we  have  denoted  PN(q^)z  =  (z^jZ^.z^.z^) .  Thus, 


M;q*)(PN(qN)z-z) |2  =  IzS- 


2 ' Z 2  ^ 1 , q*  +  lqlD  (zi'zi)  "qiDf^z3'z3^  I  0 


+  I  Z4"Z4  1 1  ,q|  +  |q2D(zrzl)+q3I,2(z3'Z3)’q2(z3'z3)l 
i  lZ2-Z2il,q*  +  ^  qf|D(z^-z3)|0)2 

+  I  z4-z4  ll,q*  +  {  1^2  I  I  D('Zl'Zl^  0  +  q3  I  °2  ^z3*z3)  I  0 

+  |qf I  I z3'z3 I  0* 
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Now  since  q  £  Q  and  Q  is  compact  by  hypothesis,  and  since 
z1,z2»z3»z4  are  by  the  choice  of  2$,  we  have 

I  a/(q*) (PN(qN)z-z)  I  2  =  0(1/N2)  for  z  6  0, 


guaranteeing  convergence  of  the  second  term  in  (3.5). 

V 

For  the  last  term  in  (3.5),  we  have  the  convergence  of  P 
to  I  strongly  on  Z.  Convergence  of  the  first  term  follows 
explicitly  from  the  form  of  the  operator, 

0 

rK  rri2  N  „  N. 

(qx-qi)(D  Z2-DZ3) 

(j^(qiN)-^'(q*))P1S(qiS)z  0 

_  (q3-q§)(D2z^VCq2-qf)CDz5I-2^) 

since  qN  — >  q*,  D2zJ  — >  D2zlt  Dz3  —  ^  Dz3,  D2Zj  — ^  D2z3, 

Dz^  — >  z i  and  z3  — >  z^. 

Thus  the  hypotheses  of  the  Trotter-Kato  hold  and  we  may  con¬ 
clude  by  it  that 

|T^(t;q^)z  -  T(t;q*)z|  — >  0  uniformly  on  [0,T] 


for  each  z  €  Z. 


□ 


AD-AU-9  050  AIR  FORCE  INST  OF  TECH  HRIGHT-PATTERSON  AFB  OH  P/E  u/, 

junE82C#j  mEcro2le?F  p<raneter  ioentifkation  for  PROBLEMS  APIS — ETC t U) 

UNCLASSIFIED  AFIT/CI/NR/82-330 

2  op  2 
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N 

Now,  since  P‘  — >  I  strongly  on  Z  (Lemma  3.4),  Theorem 
3.5  along  with  Proposition  1.4  guarantee  the  convergence  of  a  sub- 
sequence  of  solutions  q-  of  the  approximate  identification  prob¬ 
lem  to  a  solution  q  of  the  original  identification  problem  (IDA). 

We  now  briefly  describe  the  concrete  form  this  approximation 
takes  on  ZX  =  S3(AK)  x  S^(AN)  x  Sq(AN)  x  Sq(AN).  With  C^1  e  B?  3 
as  the  basis  elements  for  S3(A^)  defined  in  (1.15),  we  obtain 
as  a  basis  for  ZlN,T  the  set  where 

(C?,  0,0,0),  i  =  0 . N 

(°’Ci-(N+l) ,0’°) »  1  =  N+1 . 2K+1 

(0’°’Ci-(2N+2) *  1  =  2N+2 . 3N+2 

(0,0,0,Cj.(3N+3)),  i  «  3N+3 . 4N+3 . 

The  usual  Ritz -Galerkin  formulation  leads  to  a  4N+4 -dimensional 

\T 

matrix  system  for  the  coefficients  w  'v  ( t )  in  the  expansion  for 

N  v 

z  (t)  relative  to  the  basis  for  Z  .  In  particular,  one  finds 

(see  the  end  of  Section  2.2)  that 


M  4N  +  3  V 

Z  (t)  =  I  w-(t)^ 
i  =  0  1  1 


where  wN(t)  =  (wJJ(t)  , .  .  . , wJN+J (t) )  satisfies 


nN.N. . .  „N  N , ^  dN- 
Q  w  (t)  =  k  w  (t)  +  R  F 


with  (QN)ij  =  <Bi,Bj>,  (KN)i;j  =  <Bi ,  > ,  (RNF)i  =  <6i,F>.  This 

becomes 


(3.6) 


wN(t)  =  GNwN(t)  +  FN 


with 
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G 


N 


with 


(A?)-  • 

2Ji,] 


<C^  DC^> 


0 

I 

0 

0  ' 

aCA*)-^ 

0 

■a (A^) -1Aj 

0 

0 

0 

0 

I 

b(A2)  A3 

0 

c  (A^ ) '  1A^'-bI 

0 

N  N 
<  C  C  > 
i* Lj  O’ 

(a?).  . 

l'l.J 

=  -<DCN,DCN>n, 
i  *  j  0’ 

(A?)  . .  - 
v  3^  1J 

fj,  i  =  N+l . 2N+1 

0  ,  otherwise 


N  N  N 

where  f  =  (f0,...,f  )  corresponds  to  the  "load"  £  in  (3.1) 
and  is  given  by  fN  =  (A^)_1(RNf)  with  (ftNf) .  =  <C^',f>n. 

Also,  since  z^(t)  =  y(t,-)>  the  approximation  to  the  dis¬ 
placement  y  is  given  by 

(3.7)  yN(t,  •)  «  zj(t)  =  l  w^(t) c^T. 

1  i=0  1  1 

As  we  did  with  the  Euler-Bernoulli  equation,  we  may  rewrite 
(3.1),  uusing  a  simple  change  of  variables,  as  an  abstract  equa¬ 
tion  in  a  space  which  permits  the  use  of  lower  order  splines  for 
N 

Z  .  Again,  this  leads  to  a  different  approximation. 

Consider  the  following  change  of  variables: 


Vj  =  /a  Dy  -  /a  ip 

vz  =  yt 

Vj  =  /c  Di p 
v4  ■  *f 


Then,  (3.1)  becomes 


(3.7) 


V1 

0 

/a 

0 

° 

V2 

/a 

0 

0 

0  S 

V3 

0 

0 

0 

/c  3x 

V4. 

.  o 

0 

/c 

o 

0  0  -/a 


0  0 


0  I  v  . 


0  v. 


b//a  0  0  0  j  v. 


=  C-^Dv  +  C2v 

where  v  =  (vi >v2 » V3 »v4^ '  This  is  the  standard  form  for  a  hyper¬ 
bolic  first-order  system  of  partial  differential  equations  [37, 
p.  108]  ,  since  is  symmetric. 

This  may  be  put  into  the  previous  formulation 


z(t)  =  ja/(q)  Z  +  F(q,t)  ,  t  >  0 
z (0)  =  z0 

with  z(t)  now  taken  to  be  z(t)  =  (z^  (t)  , . .  .  ,  z4  (t) )  =  (v^ft), 

• • • ,v4(t))  in  H°  x  h°  x  H°  x  H°.  In  the  case  of  boundary  condi¬ 
tions  corresponding  to  a  fixed  beam  (y(t,0)  =  y(t,l)  *  ip  ( t ,  0)  = 
iKt,l)  =  0).  we  obtain  v2(0)  =  v2(l)  =  0  and  v4(0)  =  v4(l)  =  0. 
This  formulation  in  a  product  of  L  spaces  was  used  to  generate 
the  "data"  for  numerical  experiments  using  the  method  of  lines 
package  M0L1D  [22] . 

Since  this  formulation  does  not  yield  y(t,*)  directly,  an 
auxiliary  equation  z^  =  z2  is  included  and  integrated  along 
with  the  system  above,  yielding  y ( t , • )  =  z^(t). 
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Section  3.  Numerical  Results 

We  summarize  in  this  section  results  of  some  of  our  numeri¬ 
cal  experiments  using  the  approximate  scheme  presented  in  the 
previous  section,  based  upon  cubic  splines.  A  description  of  the 
package  which  implements  this  method  is  the  subject  of  Chapter  4. 

As  in  the  Euler-Bernoulli  case,  we  took  as  data  the  values 
of  the  solution  of  a  model  of  the  form  (3.1)  whose  parameters 
were  known  and  sought  a  solution  q^  of  the  approximate  identifi 
cation  problem  for  different  values  of  N.  The  "data"  {y^}, 

i  =  l,...,r;  j  =  1,...,£  were  generated  in  this  case  using  the 
2 

L  formulation  of  the  Timoshenko  equations  (3.7)  and  using  a 

general  purpose  computer  code  (MOL1D)  to  solve  this  system  of 

first-order  hyperbolic  equations  to  obtain  displacements  y(t,x). 

Then  (y,.},  i  =  l,...,r;  j  =  with  y .  .  =  y(t.,x.)  were 

^  J  J  1  J 

used  as  the  observations  or  input  to  the  approximate  identifica¬ 
tion  package.  For  the  experiments  listed  below,  we  used  £  =  9, 
r  =  10  and  generated  the  "observations"  at  x^  =  j/10, 
j  =  1,...,9  for  times  t^  =  i/10,  i  =  1,...,10. 

The  optimization  algorithm  requires  an  initial  guess  for  the 
N  0 

parameters  q  *  ,  which  is  referred  to  as  the  "start  up"  values 
in  the  tables  below.  The  values  of  the  parameters  to  which  the 
optimization  algorithm  converged  for  a  given  N  are  denoted  by 
qN,  and  J(qN)  is  the  cost  functional  for  those  values  of  the 
parameters . 

Example  3.1.  We  consider  the  motion  of  a  beam  initially  at 
rest  with  fixed  ends  and  described  by  the  system 


■*'9*»*S**&k  ■■  tfW  ♦  -  V. 
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tt 

■  Vxx 

•  c«l*x  * 

10e 

tt 

=  Vxx 

*  q2Cyx  ■ 

-  40 

ytt.O)  =  y(t  ,1)  =  ^ (t , 0)  =  *(t,l)  =  0 
y(0,x)  =  yt(0,x)  =  4/(0, x)  =  4/t(0,x)  =  0. 

The  numerical  results  are  given  in  Table  3.1. 


N 

-N 

^1 

Table  3.1 

3 

3 

J(qN) 

6 

.9882 

726.19 

3.6864 

.526  x  10'4 

8 

.9969 

781.00 

3. 9218 

.735  x  10‘5 

10 

1.0009 

794.29 

3.9684 

.108  x  10'5 

12 

1.00036 

794.90 

3.9732 

.165  x  10'6 

16 

1.C0033 

797.85 

3.9883 

.256  x  io'7 

TRUE  VALUE 

1.0 

800. 

4.0 

START  UP 

.9 

1000. 

3.9 

When  only  one 

parameter 

was  treated 

as  unknown 

and  the  other 

were  held  fixed  at  their  "true  values"  in  Table  3.1,  the  follow¬ 
ing  results  were  obtained. 


..-If-  «*  ~ 


Table  3.2 


N 

J(qN) 

3 

.9 3922 

.122  x 

10’1 

4 

.98452 

.105  x 

ID’2 

5 

.99533 

.211  x 

10"3 

6 

.99824 

.581  x 

10'4 

8 

.99975 

.  770  x 

10'5 

TRUE  VALUE 

1.0 

START  UP 

.9 

Table 

3.3 

N 

-K 

^2 

J(qN) 

4 

814.85 

.109  x 

10'2 

5 

804.07 

.215  x 

10'3 

6 

801.47 

.  587  x 

io*4 

8 

800.21 

.  772  x 

10'5 

TRUE  VALUE 

800. 

START  UP 

1200. 

Table 

3.4 

N 

J(qN) 

3 

3.77091 

.130  x 

10'1 

4 

3.91936 

.108  x 

io-2 

5 

3.97700 

.214  x 

10'3 

6 

3.99157 

.  58  6  x 

10'4 

8 

3.99917 

.772  x 

10'5 

TRUE  VALUE 

4.0 

START  UP 

3.8 
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Example  3.2.  We  consider  a  clamped  beam  deformed  to  the 
shape  <t>(x)  =  cos  Ax  +  cosh  Ax  -  K(sin  Ax  +  sinh  Ax)  with 
A  =  4.730,  K  =  (sin  A  +  sinh  A)/(cos  A  +  cosh  A),  then  allowed  to 
vibrate  freely,  which  can  be  described  by  the  system 

ytt  =  Vxx  '  Vx 
*tt  "  ^3^xx  *  M^x  ■  « 

y (o ,x)  =  4>(x) ,  yt(o,x)  =  o,  i^(o, x)  =  <j>'(x),  c o , x )  =  o, 
y(t,0)  =  y(t  ,1)  =  ip(t,0)  =  ip  Ct ,  1)  =  0. 

The  numerical  results  are  given  in  Table  3.5. 


Table  3.5 


N 

!? 

Icr 

J(qN) 

4 

1.1812 

1325.3 

3.6437 

.613 

X 

10'2 

8 

.  9480 

1125.1 

3.993 

.167 

X 

10'2 

10 

.9938 

1222.3 

4.1074 

.477 

X 

10'3 

12 

.9510 

1134.7 

4.1053 

.468 

X 

10'3 

16 

.9908 

1152.3 

3.8875 

.242 

X 

10*3 

20 

.9959 

1181.7 

3.9593 

.145 

X 

10"3 

24 

1.0009 

1193.7 

3.9771 

.122 

X 

10‘3 

TRUE  VALUES 

1.0 

1200. 

4.0 

START  UP 


.9 


800. 


3.8 
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CHAPTER  4.  IMPLEMENTATION  OF  THE  APPROXIMATE 
IDENTIFICATION  PROBLEM 


Section  1.  General  Description  of  the  Codes 


The  final  chapter  describes  how  the  previously  discussed 
methods  were  implemented  into  computer  codes.  While  we  have  a 
different  computer  program  for  each  of  the  methods  discussed  in 
Chapters  2  and  3,  the  basic  structure  and  much  of  the  code  is 
the  same  for  each  program.  In  fact,  the  computer  codes  used  for 
the  identification  problems  discussed  here  were  developed  from 
codes  written  for  the  identification  problem  associated  with  a 
wave  equation  model  examined  in  [  7 )  . 

The  computer  codes  were  written  with  flexibility  as  the  prim¬ 
ary  guiding  principle.  When  more  than  one  approach  could  be  used, 
the  one  that  would  most  easily  be  adapted  to  handle  extensions 
or  modifications  of  the  current  problem  was  chosen. 

We  describe  below  the  algorithms  in  a  general  setting. 

Most  of  what  follows  is  also  applicable  to  codes  used  for  identifi¬ 
cation  problems  associated  with  the  one-dimensional  hyperbolic 
and  parabolic  equations  in  [  7  ]  and  with  convection-diffusion 
equations  discussed  in  [  5 ] .  All  of  the  codes  were  written  in 
FORTRAN  and  implemented  on  the  IBM  370  at  Brown  University. 

While  this  chapter  is  intended  as  a  description  of  algori¬ 
thms  used  and  not  as  documentation  for  the  computer  codes,  refer¬ 
ence  is  made  to  the  specific  subroutines  which  implement  the 
various  algorithms. 


We  begin  by  outlining  the  general  program  for  estimating  the 
solution  to  the  identification  problem  (IDA).  We  have  seen  two 
formulations  of  the  Euler-Bernoulli  equation  and  one  of  the 
Timoshenko  equations  for  transverse  vibrations  of  a  beam.  All 
gave  rise  to  an  abstract  equation  of  the  form 


(4.1) 


z(t)  =  Mq)z(t)  +  F(q,t),  t  >  0 
z(0)  =  zQ(q), 


where  z(t)  =  (z1(t)  , . . .  ,zg(t))  €  Z,  zQ(q)  =  Y(x;q).  In  every 
case,  we  have  y(t,*;q)  =  z^(t;q)  where  y(t,x;q)  is  the  trans¬ 
verse  displacement  of  the  beam  at  x  6  [0,1]  and  at  time  t 

p 

associated  with  the  parameters  q  €  Q  a.  R  . 

In  order  to  discuss  the  common  computational  features  of 
each  of  the  approximate  identification  problems  discussed  in 


Chapters  2  and  3,  we  note  that  in  each  case  we  took  Z 


N 


N  N  N  M  N 

Z1  x  ...  x  Zj  where  Z ^  possessed  a  basis  aPPro" 


priately  modified  B-splines;  the  basis  for  Z^  then  was 


{6i}i=o>  °  2  S(N+1)  •  where  8a-i)(N+i)+j(x)  =  Bj(x)ei’ 


i  =  i,...,s;  j  =  0,...,N,  with  e^  the  usual  basis  for  R' 


In  each  case,  the  choice  of  =  p^_afp^  led  to  an  approxi¬ 


mation  of  y(t,x)  of  the  form 


N 


(4.2) 


yN(t,x)  =  l  w^ (t)B^(x) 
i  =  0  1  1 


where  wN(t)  =  (w\  (t) , . . . ,w^(t) )  satisfies 
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(4.3) 

with 


QNwN(t)  =  KNw ( t )  +  RNF(q,t),  t  >  0, 
wN(0)  =  PNV , 


N  N  N  N  N  N 

(Q  )Aj  =  (kn) i;j  = 

V  N  N  - 1  N 

(R  F).  =  <6 i , F> ,  P  f  =  (Q1  )  R1^). 


This  formulation  includes  the  following  three  problems  from 
the  previous  chapters: 

1.  The  Euler-Bernoulli  equation  with  boundary  conditions  of 

type  k,  with  approximations  based  upon  quintic  splines.  In  this 

N  N 

case,  s  =  2,  and  B ^  =  B ^  ^  satisfying  boundary  conditions  of 

*  2  0 

type  k.  The  inner  products  are  those  in  Z  =  H^Ca)  *  H  ,  namely 
<u,v>  =  <ulfVl>2>a  ♦  <u2,v2>0. 

2.  The  Euler-Bernoulli  equation  with  boundary  conditions  of 


type  1  (simply  supported),  with  approximations  based  upon  cubic 

N  N  N 

splines.  In  this  case,  s  =  3,  and  =  3  with  B^(0)  = 

B^(l)  =  0.  The  inner  products  are  those  in  Z  =  H°  *  H°  *  H° 

namely  <u,v>  =  <u1>v1>0  +  <u2,v2>0  +  <U3'V3>(T 

3.  The  Timoshenko  equations  with  fixed  end  conditions. 


with  approximations  based  upon  cubic  splines.  In  this  case,  s  =  4 
and  Bi  E  Bi  3  with  B^(0)  =  B^(l)  =  0.  The  inner  products  are 
those  in  Hg(a)  *  x  hq(c)  *  namely  <u,v>  =  <u1»vj>j  a  + 

<U2'V0  +  <U3-Vl,c  +  ^’Vo- 

Then,  given  data  {y^.},  i  =  l,...,r;  j  =  1,...,H  correspond¬ 
ing  to  displacements  of  a  beam  cross  section  located  at  Xj  G  [0,1] 

at  time  t^,  the  approximate  identification  problem  is  to  find 
— N 

q  which  minimizes 
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(4.4)  Aq)  5  J(q,yN(t,-;q),n) 

=  I  I  |yK(ti,x. ; q )  -  y..  |2, 
i=l  j-1  1  3  1J 

N 

subject  to  t  » — >y  (t,*)  satisfying  (4.2)  and  (4.3),  where 
Hi  =  (y(ti,x1) , . . .  ,y(ti,xjl)) . 

So  any  implementation  of  the  approximate  identification  prob¬ 
lem  requires  an  iterative  procedure  to  minimize  $N(q)  and  a 
method  for  approximating  solutions  to  the  system  of  ordinary 
differential  equations  (4.3). 

We  discuss  the  common  features  of  the  computer  program  for 
each  of  the  above  problems. 

As  mentioned  above  there  are  two  major  tasks  involved  in  the 

approximate  identification  program.  The  first  is  an  unconstrained 

minimization  problem  for  the  sum  of  squares  $N(q)  = 

.1  ei(q)Tei(q),  where  (q)  =  yN(ti,xj;q)  -  y^.  ,  ei  = 

(e^  ,  .  .  .  ,  eim)  .  and  yN(t^,*;q)  is  the  solution  of  (4.2)  - 

(4.3).  This  task  is  solved  efficiently  by  the  Levenberg- 

Marquardt  algorithm,  and  IMSL's  version  ZXSSQ  has  been  used  for 

this  purpose.  The  second  major  task,  required  each  time  we  evalu- 
N 

ate  4>  (q) ,  is  to  solve  numerically  the  system  of  ordinary  dif¬ 
ferential  equations  (4.3).  Since  a  wide  variety  of  problems  have 
been  solved,  some  of  which  were  stiff  (for  example,  the  Euler- 
Bernoulli  model  with  structural  damping  coefficient  6  >  0) ,  a 
general-purpose  variable- step ,  variable -order  method  was  required 
which  would  handle  both  stiff  and  non-stiff  equations  efficiently. 
Gear's  algorithm,  with  a  switch  between  Adams-type  and  backwards- 
difference  methods  was  suited  to  this  task;  IMSL's  version  DGEAR 
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was  used.  Both  of  these  algorithms  will  be  described  below. 

In  addition,  integration  of  the  system  (4.3)  requires  an 
efficient  method  for  performing  each  of  the  following  subtasks: 

i.)  Compute  the  basis  elements.  This  requires  an  effici¬ 
ent  method  for  evaluating  the  modified  B-splines 

N  N  N 

B.(x)  and  their  derivatives,  where  (B.(x)).  n  are 

splines  (linear  combinations  of  the  standard  B-splines) 

which  satisfy  prescribed  boundary  conditions. 

N  v 

ii.)  Compute  the  matrices  Q  and  K'  in  (4.3).  This  in¬ 
volves  evaluation  of  inner  products  of  the  form 
a^  =  <D^B^ ,DVB^ >p ,  with  y,v  €  (0,1,2),  and  storing 
the  resulting  matrix  (a^)  efficiently. 

iii.)  Computing  projections.  hre  need  to  project  the  initial 

’  a  .. 

function  P^T  =  ][  b.g‘.  ,  with  b  =  (b.  ,  .  .  .  ,b  )  given 

i=0  11  1  0 

by  b  =  (QN)  1R^(4'),  and  similarly  we  require  P^F(t) 

where  F  is  the  non-homogeneous  term  in  (4.1). 

N  N  v 

iv.)  Compute  the  "spline  series”  l  w.B.(x). 

i  =  0  1  1 

v.)  Evaluate  the  right-hand  side  of  the  system  (4.3),  namely 
f(w,t)  =  (QN) " 1 (KNw(t)  +  RNF(q,t)). 

Subtasks  i.)-iv.)  comprise  the  "spline  package",  discussed  in  Sec¬ 
tion  5.  The  subtask  v.)  is  discussed  in  Section  6,  where  alterna¬ 
tive  ways  of  performing  the  computations  are  discussed.  First, 
we  give  an  overview  of  the  program  structure. 
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Section  2.  Program  Structure 


We  give  an  overview  of  the  program  structure.  We  have  noted 
previously  that  macroscropical ly  the  program  may  be  described  as 
an  optimization  algorithm  for  minimizing  (q) ,  combined  with 
an  approximate  method  for  integrating  the  system  of  ordinary  dif¬ 
ferential  equations  (4.3)  in  order  to  evaluate  4>^'(q).  Since  both 
of  these  tasks  are  solved  by  iterative  methods  and  can  require 
considerable  computation,  it  is  useful  when  analyzing  efficiency 
to  divide  the  program  into  a  structure  on  three  different  levels: 

Level  1.  Preprocessing,  or  computations  which  can  be  done 
prior  to  iteration.  These  include  reading  data 

\T 

{ y ^ j  } ,  computing  and  storing  inner  products  in  Q‘ 
and  k\  and  computing  projections  which  do  not 
depend  on  the  parameter  vector  q. 

Level  2.  Computations  inside  the  iterative  loop  for  minimiz¬ 
ing  $^(q) .  These  include  projections  of  the  ini¬ 
tial  values  (which  may  depend  on  q)  ,  solving  the 
system  (4.3),  and  computing 

v  N  v 

y  (x,t)  =  l  (t) B.  (x) 
i  =  0  1  1 


Level  3.  Computations  which  are  performed  at  each  step  of 
the  algorithm  for  integrating  the  system  (4.3); 
these  computations  are  those  required  to  evaluate 
the  right-hand  side  of  (4.3). 

The  levels  are  connected  in  the  following  manner: 
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1 

_L 

2XSSQ 

T 


DGEAR 

I 

3 

v 

where  ZXSSQ  is  the  routine  to  minimize  4  (q)  ,  and  DGEAR  numeri¬ 
cally  integrates  system  (4.3).  For  a  particular  N  (dimension 
of  approximating  subspaces),  computations  at  level  1  need  only  be 
performed  once.  Computations  at  level  2  must  be  done  at  each 
iterative  step  of  the  optimization  routine  ZXSSQ.  Since  each 
iterative  step  of  the  optimization  algorithm  requires  that  the 
system  (4.3)  be  solved  (at  least  once  per  step),  and  the  numeri¬ 
cal  algorithm  for  integrating  (4.3)  requires  many  evaluations  of 
the  right-hand  side  of  (4.3),  computations  at  level  3  are  done 
most  frequently. 

In  terms  of  routines  in  the  packages,  the  routines  at  level  1 
are  those  called  from  the  main  routine.  These  include  READ,  which 
reads  input  data  and  program  parameters,  and  SETVL  which  sets 
the  knots  and  computes  inner  products  and  projections  which  do 
not  depend  on  q  by  calls  to  the  spline  package. 

The  computations  at  level  2  include  those  performed  in  sub¬ 
routine  FUNC  and  the  routines  called  from  FUNC,  such  as  INIT, 


« «  ►c**  1  W«W~  » 
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which  computes  the  projections  of  the  initial  data,  and  CNVRT 

N 

which  computes  Iw^B^(x)  by  calls  to  the  spline  package. 

The  computations  at  level  3  are  those  done  exclusively  in 
routine  DEV. 

Section  3.  Optimization  by  the  Levenberg -Marquardt  Method 

We  briefly  describe  the  Levenberg -Marquardt  algorithm.  We 
wish  to  find  q*  which  minimizes 

$(q)  =  j  eT(q)e(q) 

where  e(q)  is  the  vector  whose  components 

e(i-l)£+j(q)  =  y(ti’xj;q)  '  yij ’  1  =  1 . r;  j  =  1 . Z 

are  the  pointwise  errors  (residuals)  stored  as  a  vector.  We 
construct  a  sequence  of  iterates  q^  such  that  q^^  — >  q* 
by  using  an  iteration  scheme  of  the  form 

(4.8)  q(j+1)  =  q(j)  +  d(:i) 

where  d^^  is  a  solution  of 


(4.9) 


AU)d(J)  =  . 


V$(qCj))  , 


where  V  denotes  the  gradient  (w— , .  . .  )  and 

dql  9qp 


>0) 


A*--"  is  a 

matrix  characteristic  of  the  method. 

If  =  I,  the  identity  matrix,  then  the  method  is  steep¬ 

est  descent.  If  A^  =  V2$(q^),  where  V24>  denotes  the 


Hessian  matrix  with  entries 


82$ 


then  the  method  is  Newton's 


SqTTq" 

method  obtained  by  truncating  the  Taylor  series  expansion  after 


■feifeiiii 
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the  quadratic  term.  While  the  former  method  is  a  strict  descent 
method  (<S>(q^+^)  <_  $(q^)  and  is  theoretically  globally  con¬ 
vergent,  in  practice  convergence  is  both  slow  and  unreliable 
[17,  p.  18],  On  the  other  hand,  Newton's  method  converges  quad- 
ratically  in  a  neighborhood  of  q*  [17,  p.  35],  but  it  requires 
evaluation  of  the  Hessian  matrix,  which  makes  it  impractical  for 
the  problem  at  hand. 

However,  since 


(4.10) 


V$(q)  =  G(q)TG(q) 


where  G(q) 


is  the  Jacobian  of  e 


(4.11)  V2$(q)  -  Iei(q)V2e.(q)  +  G(q)TG(q), 

2  T 

we  may  approximate  V  4(a)  by  G(q)  G(q)  for  small  residuals. 
Thus  (4.8)  and  (4.9)  with 


(4.12) 


=  G(q^^)TG(q^^) 


is  the  Gauss-Newton  method.  Convergence  of  this  method  is  at 
least  linear  in  a  neighborhood  of  q*,  and  quadratic  if  4(q*)  =  0 
[17,  p.  94].  However,  this  method  experiences  problems  if  G  is 
ill-conditioned  (there  is  no  guarantee  that  G  even  be  full  rank 
away  from  q*) .  The  Levenberg-Marquardt  approach  is  to  replace 

(4.12)  with 

(4.13)  A ^  =  y^D  +  G(q^))TG(q^))  , 

where  D  is  a  positive  diagonal  matrix,  with  y^  0  chosen 
sufficiently  large  to  ensure  A^  is  positive  definite  and  to 
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ensure  strict  descent  (4>(q^+^)  <  4>(q^^)).  In  Marquardt  and 

Levenberg's  original  scheme,  D  =  I  was  used.  The  IMSL  imple- 

T 

mentation  ZXSSQ  uses  =  (G  G)^. 

Note  that  as  u  ^ becomes  large,  the  direction  of  search 

d^  approaches  that  of  steepest  descent,  while  p  ^  =  0  is 

the  Gauss-Newton  method.  The  method  is  locally  convergent,  and 

the  rate  of  convergence  is  at  least  linear  [17,  p.  96  ]  . 

3e-(q) 

Since  G--  =  — 4 -  is  not  generally  available  analytically, 

ij  aq^ 

it  is  computed  numerically.  An  option  in  ZXSSQ  permits  one  to 
have  the  Jacobian  computed  numerically  by 


where  <5 ^  =  max ( j I , . 1) */u,  where  u  is  the  relative  precision 
of  floating  point  computations.  While  this  avoids  the  problem 
of  having  to  provide  an  analytic  Jacobian,  it  costs  p  addi¬ 
tional  evaluations  of  4>(q)  to  numerically  compute  G.  Thus  when 
G  is  fully  evaluated,  p+1  evaluations  of  ^(q)  are  required 
per  iteration. 

If  the  Jacobian  is  not  changing  too  rapidly,  it  is  possible 
to  approximate  G  using  the  information  in  the  direction  of  the 
most  recent  step  to  update  the  Jacobian.  ZXSSQ  uses,  when  ap¬ 
propriate,  a  rank-one  update  of  the  form 

(4.15)  G  ^  +1^  =  G(j)  +  — i— [e  (q^ +1^ )  -  e(q(:i))  -  G^aJa1 

a  a 

where  a  =  q^+^  -  q^  .  When  this  approximation  to  the  Jacobian 
is  used,  no  additional  evaluatons  of  4>(q)  are  required,  and  so 
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<t(q)  will  be  evaluated  only  once  in  such  an  iteration  step. 

The  iteration  (4 . 8) - (4 . 9) - (4 . 13)  is  repeated  until  one  of 
the  following  conditions  is  met: 


1. 

|V$(q(J)) | 

<  6 

2. 

|*(q^+1^) 

-$(q(;)’ 

^  )/4(q^  ^ )  |  <  e 

3. 

<5>(q^  ) 

and 

4>(q^)  agree  to  NSIG  significant 

digits . 

Section  4.  Integrating  the  System  of  Ordinary  Differential 
Equations 

We  next  describe  very  briefly  the  method  used  for  integrating 
the  system  (4.3),  namely  IMSL's  routine  DGEAR  (23).  For  a 
description  see  [21]  .  We  only  include  those  details  which  afffect 
the  choice  of  parameters  for  the  ID  package  and  provide  some  ap¬ 
preciation  for  the  amount  of  computation  required. 

To  abbreviate  the  notation,  write  the  system  (4.3)  in  the 

form : 


(4.17) 


*  *  f (w,t) ,  w(0)  =  w0. 


The  methods  in  DGEAR  are  based  upon  difference  approximations  of 


the  form 


(4.18) 


K  K 

Kk  ’  Tj  Vk-j  *  h  $  eA-r  6o  > 


where  ,  8^ 

and  where  w. 

k 


are  constants  associated  with  a  particular  method 


is  an  approximation  to  w(t^) ,  w^ 


f(wk»tk)  is 


an  approximation  to  and  h  is  a  constant  stepsize 

(tjs+^  =  t^  +  h)  .  The  Adams  methods  of  order  m  correspond  to 
the  values  *  1,  =  m-1,  and  the  backwards  difference  for¬ 

mulas  (BDF)  correspond  to  the  values  *  m,  ^  -  0.  If  (4.18) 
is  solved  with  all  past  values  exact,  then  w,  -  w(t.)  =  0(hm+*) 
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for  small  h,  and  hence  the  method  is  said  to  be  of  order  m. 

As  (4.18)  is  generally  implicit,  it  must  be  solved  itera¬ 
tively.  If  we  write 

K  K 

(4.19)  g(wk)  -  wk  -  h90f(Vtk)  -  I1  ajwk_j  -  h  I2  S^.j- 


then  (4.18)  is  the  implicit  equation  g(w^)  =  0.  DGEAR  provides 
a  variety  of  methods  which  the  user  may  choose  to  solve 
g(wk)  =  0. 

First  a  predicted  value  of  w£^  is  obtained  by  an  expli¬ 
cit  method.  This  is  equivalent  to  (4.18)  with  BQ  *  0,  but  DGEAR 
uses  the  Nordsicck  formulation  (see  [21)).  Then  g(w^)  =  0  is 
solved  iteratively  by  one  of  a  class  of  methods  of  the  form 


(4.20) 

Note  if  P 


w'i+1)  -  .»>  -  (p(iV1g(w(i)) 


(i) 


-  tjf) 


"  1  '  h60(W) 


(W(i))»  then  (4.20) 


is  Newton's  method.  If  p£i)  =  then  (4.20)  is  functional 
(fixed  point)  iteration. 

Between  these  two  extreme  lie  various  approximations  to 

(4&)  ,  (i),  which  are  more  cheaply  evaluated.  One  choice  is 
°w  *-wk 

Pk*^  *  p£^  ;  i.e.,  do  not  re-evaluate  the  Jacobian  at  each  itera¬ 
tion.  The  so-called  chord  method  uses  P^  =  P^ ,  for  some 
k’  <_  k,  corresponding  to  the  parameter  MITER  =  1  or  MITER  =  2. 

A  still  less  costly  method,  particularly  where  large  systems  are 

d  f 

concerned,  is  one  in  which  is  approximated  by  a  diagonal 

matrix  whose  entries  are  forward  difference  approximations  to  the 
diagonal  of  the  Jacobian  (MITER  =3).  If  is  not  avail¬ 

able  analytically,  it  may  be  computed  numerically  (MITER  =  2). 
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For  most  of  the  examples  considered  here,  the  Adams  methods 
(METH  =  1)  with  functional  iteration  (MITER  =  0)  were  most  effici 

V  V 

ent.  However,  it  should  be  noted  that  the  P  methods  can 

lead  to  moderately  stiff  ODEs,  as  in  the  case  where  structural 
damping  was  included,  in  which  case  the  BDF  methods  (METH  =  2) 
and  numerically  computed  Jacobian  (MITER  =  2)  were  used. 

While  the  methods  (4.18)  are  based  upon  constant  stepsize  h 
DGEAR  does  adjust  the  stepsize  and  order.  Briefly,  this  is  ac¬ 
complished  by  interpolating  the  past  data  w,  - ,w,  .  (using  the 

K  -  j  k  -  j 

Nordsieck  array)  to  obtain  the  values  of  w,  required  at 

k  -  j 

t^  -  jh*  where  h*  is  the  new  stepsize.  Following  a  step  of 
size  h  at  order  m,  DGEAR  attempts  to  readjust  the  stepsize  up¬ 
ward  every  m+2  steps;  this  is  done  by  estimating  the  local  trun 
cation  error  at  orders  m-1,  m,  and  m+1  and  choosing  h*  to  be 
the  largest  permitted  by  these  three,  and  the  order  m  is  reset 
accordingly.  The  estimated  truncation  error  is  effectively  com¬ 
pared  to  the  requested  local  truncation  error  bound  TOL  which  is 
specified  by  the  user. 

One  further  parameter  must  be  specified,  namely  an  initial 
stepsize  H0 .  Since  the  method  starts  out  with  a  first  order  ap¬ 
proximation  m  =  1  (no  previous  history  available) ,  H0  must  be 
small.  For  the  ID  problems  considered,  TOL  was  usually  taken  to 
be  1.  x  io  6  and  H0  to  be  1.  *  10  ^ . 


Section  5.  The  Spline  Package 

The  spline  package  contains  the  routines  for  evaluating 
splines  and  their  derivatives,  for  computing  projections  onto 
spline  subspaces,  for  computing  inner  products  of  the  form 
<DkB^  n,DkB^  n>Q»  0  <  k  <  n,  and  for  evaluating  the  "spline 
series"  Iw.B^n(x). 

a)  Evaluation  of  Modified  B-Splines  and  Their  Derivatives 

The  evaluation  of  the  B-spline  basis  functions  is  of  course 
fundamental  to  all  the  routines  in  the  spline  package.  The  B- 
spline  basis  functions  were  introduced  in  Section  3  of  Chapter  1. 
Recall  that  there  we  defined  all  the  B-splines  in  terms  of  a  singl 
function.  We  defined 

(4.21)  Bn(x)  =  ^(x-y)^'1  *  (*1)j(j)(x'j)+  1* 

Then,  the  B-splines  of  odd  degree  (even  order)  were  obtained  by 

taking  n  =  2m  (m  =  2  for  cubics,  m  =  3  for  quintics)  and  a 

partition  =  (x-}^+?  k  with  x.  =  i/N;  we  defined 
r  i  i=l-m  i 


Yi  ■  xi  -  mh 

B^,n_i (x)  =  In((x-yi)/h) ,  i  =  1-m . N+m-1. 

Finally,  the  modified  B-splines  B^  , (x)  are  computed  by  taking 
the  appropriate  linear  combinations  of  the  as  in  (1.15)- 

(1.18)  so  that  the  B^  n_^  satisfy  the  given  boundary  conditions 
F  is  the  fundamental  spline  of  order  n,  with  knots  at 
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0,1,..., n,  and  with  support  on  [0,n].  The  splines  B.  have 

i ,  n 

knots  and  the  support  of  each  iL  2m-l  ^xi-m,xi+m^ 

(or  [>Vyi  +  2m]>- 

N 

For  the  uniform  partition  A  ,  two  methods  of  evaluating  the 

B.  i (x)  may  be  used.  The  first  method  uses  the  piecewise 
l , n- 1 v  '  ' 

polynomial  representation  of  B^  n For  example,  in  the  case 
n  =  4  (cubic  splines),  we  obtain  (see  [35,  pp.  89-90]) 


(x-xi_2)3>  if  x  €  [xi.2,xil] 

h3+3h2 (x-xil) +3h(x-xi ^) 2-3 (x-xi_1) 3, 

if  x  e  txi.1,xi] , 

h3+3h2(xi+1-x)+3h(xi+1-x)2-3(xi+1-x)3, 

if  X  £  [xi,xi  +  1]  , 

Cxi  +  i ~XD 3  ,  if  X  e  txi  +  i»xi  +  2] 

0,  otherwise. 


When  each  of  the  above  expressions  is  evaluated  using  Horner's 
scheme,  this  is  a  very  efficient  method,  requiring  only  four 
multiplications.  The  derivatives  of  the  B-splines  may  be  repre¬ 
sented  similarly.  However,  a  different  representation  for  each 
n  and  for  each  derivative  is  required. 

The  algorithm  we  have  used  for  the  computations  in  Chapters 
2  and  3  may  be  used  for  splines  of  arbitrary  order  and  for  deri¬ 
vatives  of  any  order.  This  algorithm  is  based  upon  the  iterative 
formulas  satisfied  by  the  fundamental  B-spline  (see  [42,  p.  136]) 

Bn(x)  =  xBn.x(x)  +  (n-x)Bn_1(x-l) 


(4.22) 

and 
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(4.23)  DBn(x)  =  (n-1)  [Bn.1(x)  -  B^fx-l)]. 

The  iterative  formula  (4.22)  can  be  used  to  evaluate  Bn(x) 
for  any  x.  Since  we  evaluate  splines  only  when  the  values  of 
all  non-zero  splines  at  x  are  required,  as  when  computing  the 
matrix  whose  entries  are  <B^  ,,B^  .  >n  for  0  <  |i-j|  <  n, 

j  =  0,...,N,  the  algorithm  may  be  refined  to  efficiently  produce 
the  values  of  all  non-zero  splines  of  order  n  at  x. 

First  note  that  B^  (x)  =  Br  ( (x-yi  +  j  ) /h)  = 
Bn((x-yi)/h-j).  Given  n,  find  i*  such  that  n-1  (x-y^*)/h£n, 

and  define  £  =  (x-y^*)/h.  The  values  of  the  non -vanishing 
B-splines  (B^*+j  n-i(x)^j=o  then  correspond  to  the  values  of 
B n ( C - j ) ,  j  =  0 , . . . ,n-l .  To  compute  B n ( C - j ) ,  we  note  that  (4.22) 
implies 


(4.23)  BjU-i)  =  (5-i)Bj_1(C-i)  +  (i+j-OB.^^-l-i)  , 

for  j  =  2,...,n;  i  =  n-j,...,n-l. 

Now,  since  B^fS-j)  =  the  iterative  formula 

(4.23)  may  be  used  to  evaluate  Bn(£-j),  j  =  0,...,n-l.  Then, 

Bj  +  i*,n-i(x)  =  Bn(C-j),  j  =  0  , . . . ,n-l . 

If  derivatives  of  the  B^  ,  are  required,  as  when  the  ma- 

j,n-l 

7  \T  o  xt 

trix  with  entries  <D  Bv|  S,D  B?  5>Q  is  to  be  computed,  then  (4.22; 
may  be  used  in  conjunction  with  the  above  algorithm.  We  compute 
Bn_i(£-j),  j  =  0 , . . . ,n-l ,  then  compute  DBn(g-j)  =  (n-1) (If  j (£- j ) 

-  ^n-i (5* J *1) ) •  For  second  derivatives,  D2Bn(C-j)  = 

(n-1) (n-2) (Bn  2(C-j)  -  2 B^  +  B^  ^ (5- j -2) ) .  Note  that  in 

general  to  compute  D  Bn(£-j)  k  <  n-1,  we  need  only  evaluate 


i 


no 


i 


If  ^(S-j)  and  perform  k  differences  on  these  values  and  mul¬ 
tiply  by  (n-1) ! / (n-k- 1) ! .  Then,  DkB^+i* >n. 1 (x)  =  DkBn(£-;j). 

N 

Finally,  the  B.  are  computed  using  (1 . 15) - (1 . 18)  when 
J  >n 

appropriate.  We  note  that  each  of  these  modification  formulas 

N  N  N  N 

requires  only  that  Bq ,B^ ,B^_^ ,B^  be  modified.  Moreover,  these 
are  changed  only  when  x  £  [0,x2]  or  x  £  tx^ _ 2 » '  s:*-nce  t^ie 
functions  used  to  modify  vanish  outside  these  intervals.  So  we 
first  check  if  i  =  0,1,N,N-1  and  x  lies  in  one  of  the  two 
intervals.  If  so,  we  compute 


or 


N  «M  ~  N 

B?(x)  =  bJ(x)  -  ci>1B1^1(x)  -  ci>2B^2(x),  i  =  0,1 
BiU)  =  B^'(x)  -  cij3Bjj+1(x)  -  cj  }4®k+2  »  1  =  N'1  or  N» 


where  the  c^  are  the  coefficients  appearing  in  (1 . 15)  -  (1 . 18)  . 

Remark .  The  algorithm  presented  for  computing  the  unmodified 
B-splines  presented  above  is  essentially  the  one  in  [42,  p.  205] 
and  is  attributed  to  DeBoor. 

The  algorithm  to  compute  the  n  non -vanishing  B-splines  at 

2  2 
x  requires  n  +  2n-l  multiplications,  as  compared  to  n  -n 

multiplications  for  the  piecewise  polynomial  representation.  How¬ 
ever,  the  flexibility  in  computing  derivatives  and  splines  of  all 
order  in  one  algorithm  make  the  iterative  algorithm  attractive. 
Furthermore,  the  spline  computations  occur  primarily  at  level  one, 
and  so  the  slight  difference  in  efficiency  has  a  negligible  effect 
on  overall  execution  time. 

This  algorithm  appears  in  the  various  packages  as  SPQV, 

SPMVAL ,  or  SPPVAL. 


I 


Ill 


b)  Computation  of  Inner  Products 


The  computation  of  the  inner  products  is  carried  out  via  a 
composite  Gauss  -  Legendre  two-point  rule: 


(4.24) 


f1  h  K 

f(x)dx  =  j  l  f((v 

0  L  v=l 


j  -  /3/6)h)  +  f((v  -  j  +  /5/6)h) 


where  h  =  1/K,  and  K  is  chosen  to  be  a  fixed  multiple  of  N 
(usually  K  =  8N) .  The  inner  products  <B^,B^>q  are  calculated 
by  computing  xlv  =  (v  -  j  -  /3/2)h,  x2v  =  (v  -  j  +  /3/2)h,  evalua¬ 
ting  all  non-zero  splines  at  xlv>  x2v>  f°rminS  the  non-zero  pro¬ 
ducts  (xiv) (xj v)  and  accumulating  the  products  in  a  banded 
matrix  (storing  only  nonzero  subdiagonals  as  columns) .  This 
algorithm  appears  in  the  various  packages  as  SETQ ,  SETMAT, 

SETMBT . 

We  note  that  for  large  values  of  N,  it  may  be  worthwhile  to 
use  the  property  that  <B^,Bj>  =  <Bi+1»Bj+1>»  where  i  =  2,...,N-2 
for  cubic  splines  and  i  =  3,...,N-3  for  quintic  splines. 


c)  Computation  of  Projections 


The  projections  in  the  | ' l^-norm  onto  the  spline  subspaces 
are  easily  computed,  once  the  inner  products  have  been 

computed  and  stored.  Let  be  a  typical  projection  operator 


in  the 


'  m 


■norm,  m  =  0,1,  or  2.  Suppose  we  wish  to  compute  the 

'nr  aN 


projection  of  <p  onto  sk(A1N)  where  n  =  3  (k  =  0)  or  n  =  5 
(k  =  1,2,  or  3).  Then  the  projection  PN  is  given  by  P^4>  = 

_I  aiBi  where  a  =  (a0,...,aN)  is  given  by 

a  -  (QW* 


(4.25) 
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with 


IRN*).  =  <4>,Bj>m  =  <D%,DmB?>0 


Dm4 (x)  DmB^'  (x)  dx 

0 


and 


(QN)  .  .  =  <BN,BN>  =  <DmB^',DmBN>n 
^  '  ij  1  ’  j  m  l  *  j  0 


The  matrix  Q‘  is  computed  once  via  the 
stored;  thus  to  compute  the  projections, 
inner  products  <Dm4>  ,DmB^>g ,  i  =  0,...,N 
rule  (4.24]  . 


DmB^'  (x)  DmBj  (x)  dx  . 

Gauss-Legendre  rule  and 
we  need  only  compute  the 
via  the  Gauss-Legendre 


We  note  that  in  the  case  where  the  1*1  -norm  is  used,  the 

1  'm,a 

projections  are  still  given  by  (4.25);  i.e.,  the  projections  them¬ 
selves  do  not  depend  on  the  weight  a.  This  is  because  with  the 
|  •  |m  a*norrc,  a  =  (aQ^)  *aR^<j)  =  (Q^)  1R^4',  where  and  R^'  are 

as  above.  Finally,  to  compute  the  projections,  we  need  a  method 

N 

for  evaluating  £  a.B.(x). 

i  =  0  1  1 


d)  Evaluation  of  the  "Spline  Series" 


We  also  have  two  methods  for  evaluating 


(4.25) 


N 


,N 


s(x)  =  l  CAB" 

i=0  1  1,n 


i“i.n-l(x) ’ 


First  observe  that  because  of  the  local  support  of  the  B-splines, 
that  only  r.  terms  of  (4.25)  need  be  evaluated.  Thus  we  can  gen¬ 
erate  the  values  of  the  n  non-vanishing  B-splines  at  x  using 
the  previously  defined  algorithm  and  perform  the  linear  combina¬ 
tion.  A  slightly  more  efficient  method  was  pointed  out  in  [42, 
p.  193],  We  describe  this  method  briefly  here. 
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Again,  first  find  i*  such  that  n-1  <_  (x-y^*)/h  <  n,  and 
let  £  =  (x-yi*)/h.  Then, 

Six)  -  I  cJ*  j(x)  =  I  c  .  B n(E-i). 
i  =  0  1  1,n  i  =  0  1 

Define  cB  =  ci*+i,  i  =  0,...,n.  Using  (4.22),  it  is  easily 
shown  that 

n-i  n  n-1  . 

Z  cJb  (J-i)  -  I  c|b  j(E-i), 

1  =  0  1  n  i  =  1  1  n 

where  =  (£-i)c|?  +  n-^+i-lJc*?^ , 


n-1 

CiBn-2^‘^ 


i  =  . ,  .  .  .  ,n-l 


where  =  (£-i)cB  +  (n - £  +  i  - 2 ) c* . 1 , 

i  =  1 ,  .  .  .  ,n-l 


where 


(4.26) 


n-1  •_ 

I  c^B  .  (5-j), 

•  “  •  l  n-j  v  J  1  ' 


=  (C-i)cl'1  +  (n-£+i-j)c|_J 


For  j  =  n-1,  this  becomes 


S-l  ®lCC-n+D 


Thus  we  generate  the  triangular  array 
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'n- 1  » 


using  (4.26)  and  take  s(x)  =  cjj_j.  This  algorithm  requires  only 

2 

n  -n+2  operations.  To  get  s(x),  we  need  to  write  s(x)  in 

terms  of  the  unmodified  B-splines  ^  (x)  .  We  simply  find 

{c,}?:?  1  such  that 
l  l - 1  -m 


N+m-  1 


(4.27) 


i  =  l-m  Ci+m-lBi,n-l(x)  =  JQ  ciBi  ,n-l (x)  ’ 


ci+  i  =  C.  for  2  <  i  <  N-2.  The  others  are  found  using  (1.15)- 

V  ^  \r 

(1.18)  to  write  B,  „  ,  in  terms  of  the  B.  ..  The  computa- 

i  ,n -  j.  l  ,n-l 

A 

tion  of  the  coefficients  c^  need  only  be  performed  when  the  c^ 
are  changed.  For  example  in  the  case  of  cubic  splines  modified  as 
in  (1.15),  we  find  the  requirement  that 


l  c~B*  =  I  c  B * 
i= -1  11  i=0  11 


implies  o  =  ci  for  2  _<  i  <  N-2,  and 

-w  N  N 

C-1B-1  +  coBo  +  C1B1  =  C0B0  +  C1B1’ 

(recall  that  Bq  =  B^  -  4B^  and  B^  =  B^  -  Bq/4)  which  in  turn 
is  equal  to 

aW  ^  V  *  M 

'4c0B-l  +  (c0‘cl/4)B0  +  C1B1‘ 


Thus  we  find  that  c.j  =  -4cQ,  cQ  =  cQ  -  Cj/4,  c1  =  Cj,  and  simi¬ 
larly  cN1  =  Cj,  cN  =  cN  -  cn_1/4,  cn+1  =  -4cn. 
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For  the  quintic  splines,  we  find  that  corresponding  to  the 
modified  B-splines  in  (1.16)  for  simple  end  conditions,  we  obtain 


(4.28) 


£i  -  ci,  0  <  i  <  N 


C~-2  = 

12c0  - 

c2 

C~-l  = 

'3co  * 

C1 

CN+1 

1 

n 

CN-1 

CN  +  2 

12cn 

CN-2 

Corresponding  the  modified  B-splines  in  (1.19)  for  fixed  end  con¬ 
ditions,  we  obtain 


(4.29) 


c  i  ,  0  <  i  <  N 

41.2Sc0  +  32 . 5c^  +  2.25c2 

-4.12Sc0  -  2 . 2 5 c :  -  .125c2 

=  -4.125cn  -  2 . 25Cjj_  j  -  .125cN?_2 

=  41.2  5c  +  32 . 5cv  .  +  2.25cv.  , 

N  N  - 1  N  -  2 


Finally,  for  the  modified  B-splines  corresponding  to  the  canti¬ 
lever  beam  in  (1.18),  we  obtain  (4.29)  except  that 


(4.30) 


N+l 

1'5cN 

-  ,5c 

N+2 

'  3cn  ■ 

2cn-i 

Section  6.  Evaluation  of  the  Right-Hand  Side  of  the  Approximating 
System 

Since  the  most  frequently  evaluated  computations  occur  at 
level  3  where  we  evaluate  the  right-hand  side  of  the  system  (4.3) 
it  is  worthwhile  to  organize  these  computations  efficiently. 
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We  consider  ways  to  evaluate  the  right-hand  side  of  the 
system  (4.3) 

wN(t)  =  (QN)'1{KNwN(t)  +  RNF(q,t)}. 

Consider  first  the  homogeneous  part,  (Q^)  ^K^w.  Note  that 
will  have  a  block  diagonal  structure 

Qi 

Qz. 

Qs 


To  be  concrete,  let  us  restrict  ourselves  to  the  case  s  =  2, 
which  is  the  case  of  the  Euler-Bernoulli  formulation.  In  this 
case , 


qn  = 


Q?  is  (N+l)  x  (N+l) 


,  KV  is  (n+l)  x  (N+l). 


Each  of  the  submatrices  and  K?  has  a  banded  structure;  each 

will  be  7-banded  when  cubic  splines  are  used  and  11-banded  when 
quintic  splines  are  used.  Moreover,  since  each  is  a  symmetric  ma¬ 
trix,  these  can  be  stored  efficiently  by  storing  the  subdiagonals 
and  the  diagonal  as  columns  of  a  matrix  which  is  (N+l)  *  4  (or 
(N+l)  x  6  for  quintic  splines)  . 
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11^ 

Two  ways  can  be  used  to  evaluate  a  general  term  (Qp  xp, 

-IN 

where  v  =  . VN+1) .  One  way  is  to  compute  A  =  (Qp 

at  level  one  and  save  the  matrix  A.  Then  when  evaluating  the 
right-hand  side  at  level  3,  we  need  only  do  a  matrix  multipli¬ 
cation  Av.  However,  A  will  in  general  be  a  full  matrix  (no 

banded  or  sparse  structure) ,  and  so  this  multiplication  requires 
2 

(N+l)  operations. 

Another  way  to  do  this,  which  preserves  the  banded  structure, 
N  X 

is  to  factor  =  LL1  where  L  is  lower  triangular,  at  level  1 
by  the  Cholesky  algorithm,  and  store  the  factor  L  as  a  banded 
matrix  (requiring  (N+l)  x  4  or  (N+l)  x  6  locations)  .  Then 
at  level  3,  we  compute 


i)  c  ■  KNv 

ii)  L<t>  =  c 

iii)  LTx  =  5. 


Step  (i)  requires  7N-5  operations  (11N-9  for  quintics)  ,  steps  (ii) 

and  (iii)  each  require  4N-1  (6N-20  for  quintics)  to  backsolve  a 

banded  triangular  system  of  algebraic  equations.  When  done  in  this 

N  - 1  N 

manner,  the  computation  of  (Qp  Kp  requires  15N-7  opera¬ 
tions,  N  _>  4,  when  cubic  splines  are  used,  and  23N-49  operations, 

N  >  6,  when  quintic  splines  are  used.  Thus  for  large  values  of 

N  (N  >  16),  there  is  a  clear  advantage  to  the  second  approach. 

N  T 

The  Cholesky  decomposition  of  the  matrix  into  LL 

was  carried  out  by  IMSL  routine  LUDAPB,  and  the  backsolution  (ii) 
and  (iii)  by  IMSL  routine  LUELPB. 

We  can  now  compare  the  computational  efficiency  of  the  two 
approximations  discussed  for  the  Euler-Bernoulli  beam.  In  order  to 


make  a  comparison,  consider  the  case  when  y  =  0  (no  viscous 
damping) . 

For  the  quintic  spline  formulation,  we  had 

\ 

nN  _  Q‘l  0 


KN  - 


with  Q 


So  (Q^)  V  had  the  structure 


“Q21k2 


6Q21]C3 


So  with  u  =  (u1 , . . . ,uN+1) ,  v  =  (Vj , . . . ,vN+1) ,  we  compute 


cQWft)  ■  ,x  ; 

V  '  Q2 1 (aK^u  +  5K,v) 


requiring  about  34N  operations  to  evaluate  when  6  >  0  and 
about  29N  operations  when  6=0. 

For  the  cubic  spline  formulation,  we  obtained 


(QN)'1KN  =  0 


0 


-cxCA^^aJ 
0  0 


when  y,6  are  both  zero.  The  evaluation  of  (QN)_1KNfv)  requii 

\w/ 

about  30N  operations.  In  addition,  an  extra  set  of  N+l  equa- 


11 

tions  has  been  used,  compared  to  the  quintic  spline  formulation. 
Therefore,  based  upon  this  and  our  computational  experience 
there  seems  to  be  no  advantage  in  going  to  the  cubic  spline  for¬ 
mulation  . 

To  evaluate  the  non-homogeneous  term,  note  (QN) ’*R^F(q,t) 

N 

P  F(q,t).  The  projections  may  be  computed  as  in  Section  5.  If 
the  load  f(t,x;q)  can  be  written  as  g(x)h(t;q),  then  F(q,t) 

T  %r  *r 

has  the  form  (0,g)  -h(t;q)  and  the  projections  PN(0,g)  may 
be  computed  at  level  1  (once  for  any  given  N)  and  stored. 
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CHAPTER  5.  CONCLUSIONS 

We  have  presented  numerical  methods  for  parameter  identifica¬ 
tion  of  several  constant  parameters  appearing  in  the  Euler- 
Bernoulli  and  Timoshenko  equations  for  transverse  vibrations  of 
a  beam.  The  practical  utility  of  our  approach  is  supported  by  our 
computational  experience  on  a  large  number  of  examples,  a  sample 
of  which  has  been  included  here.  In  most  cases  of  interest  we  have 
given  a  complete  treatment  of  our  methods  including  numerical 
algorithms,  convergence  proofs,  and  numerical  results. 

The  fundamental  ideas  upon  which  our  convergence  results  are 
based,  involving  the  use  of  a  semigroup  theoretic  approach  for  the 
approximation  of  identification  problems  governed  by  partial  dif¬ 
ferential  equations  (distributed  parameter  systems),  first  ap¬ 
peared  in  [11].  These  methods  essentially  involve  writing  the 
initial -boundary  value  problem  as  an  abstract  equation 

(5.1)  2(t)  -^(q)z(t)  +  F  (q  ,  t) 

in  a  Hilbert  space  Z,  where  j^( q)  is  the  generator  of  a  CQ 

semigroup  in  Z,  and  approximating  the  generator  ja^(q)  by  the 
N  N  ^  N  N 

operator  srf  (q)  =  P  _of( q)P  ,  with  P  the  projection  onto  a 
finite  dimensional  subspace  ZN  spanned  by  splines.  The  result¬ 
ing  system  of  ordinary  differential  equations 

(5.2)  iN(t)  -j*N(q)zN(t)  +  FN(q,t) 
is  used  to  approximate  solutions  of  (5.1). 
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Instead  of  seeking  a  q  which  minimizes  a  cost  functional 

J(q»z(*»q))  over  mild  solutions  of  (5.1),  we  obtain  estimated 
—  — v 

values  for  q  by  seeking  a  q*  which  minimizes  the  cost  func¬ 
tional  J(q;zN(-;q))  over  solutions  zN(-;q)  to  (5.2). 

The  idea  of  estimating  solutions  to  the  identification  prob¬ 
lem  is  not  new.  The  notion  of  approximating  solutions  of  the 
partial  differential  equation  and  performing  the  optimization  on 
the  approximate  solutions  to  obtain  parameter  estimates  has  been 
used  extensively.  We  refer  the  reader  to  the  survey  articles 
12],  [26]  ,  [32] ,  and  [36]  to  see  the  variety  of  approaches  which 
have  been  used.  We  only  discuss  some  of  the  relevant  material 
here . 

We  note  that  many  investigators  have  used  finite-difference 
methods,  modal  approximations ,  and  the  methods  of  lines  to  approx 
mate  solutions  of  the  partial  differential  equation  (see  Table  1 
in  [26])  in  the  context  of  parameter  identification.  Galerkin 
methods,  to  which  our  methods  are  closely  related,  have  also  been 
used.  In  [33] ,  Galerkin  methods  are  used  for  the  heat  equation 
employing  a  basis  of  polynomials  which  satisfy  the  boundary  condi 
tions  and  which  are  orthogonalized  via  Gram  -  Schmidt . 

While  most  of  the  work  in  the  surveys  deals  with  the  estima¬ 
tion  problem  for  the  heat  equation,  several  authors  have  proposed 
numerical  methods  for  estimating  a  single  parameter  in  the  beam 
equation 


(5.3) 


ytt  =  -^xxxx- 


In  [19]  an  example  is  given  where  this  is  done  using  finite 
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differences  to  approximate  solutions  to  (5.3).  We  have  found  no 
results  in  the  literature  for  estimating  structural  or  viscous 
damping  coefficients  or  the  parameters  in  the  Timoshenko  equation. 

Cubic  splines  have  also  been  used  in  the  context  of  estima¬ 
tion  of  parameters  (see  surveys  above) .  In  [43]  they  are  used 
for  the  one-dimensional  heat  equation  to  obtain  a  "lumped  system" 
of  ordinary  differential  equations  which  approximate  the  partial 
differential  equation,  by  collocating  with  cubic  splines  in  the 
spatial  variable.  We  believe  our  approximations  using  P^_a^(q)P^ 
with  cubic  and  quintic  splines  are  new. 

Little  work  has  been  done  in  proving  convergence  of  parameter 
estimation  schemes,  and  we  have  found  no  such  theoretical  work  in 
the  literature  for  methods  comparable  to  ours.  The  literature  on 
parameter  estimation  consists  mainly  of  proposed  numerical 
methods  with  test  results  for  a  simple  example  estimating  a  single 
identification  of  many  parameters  and  have  provided  proofs  of 
convergence  for  our  methods. 

While  some  authors  have  investigated  the  identification  of 
coefficients  which  are  a  function  of  the  spatial  variable  (see 
Table  1  in  [26]),  we  have  restricted  our  attention  here  to  the 
case  of  constant  parameters.  Our  methods  do  carry  over  to  the 
case  of  spatially  varying  coefficients,  and  this  has  been  done 
for  a  special  case  of  the  convection-diffusion  equation  (see 
[5,  p.  22]  for  a  discussion).  We  have  not  treated  the  questions 
of  observability  on  identificability ,  but  have  chosen  to  empha¬ 
size  the  convergence  of  the  parameter  estimates  to  the  solution 
of  the  identification  problem. 


In  the  case  of  deterministic  data  [noise  free) ,  most  authors 


have  used  a  steepest  descent  method  to  perform  the  optimization 
(see  Table  1  in  [26]  and  other  surveys) .  It  is  well-known  that 
steepest  descent  can  exhibit  oscillatory  behavior  and  perform 
badly  in  practice  [17,  p.  18].  For  this  reason,  we  have  chosen 
the  more  robust  Levenberg-Marquardt  method. 

We  have  used  a  discrete  least  -  squares  fit-to-data  criterion 
(cost  functional) .  Many  other  type  of  cost  functionals  have  been 
used  also.  In  particular,  cost  functionals  of  the  form 


£  fT 

(5.4)  J(q,y(-,-;q))  =  l  |y(t,x.)  -  n-(t)|2dt 

i=l  Jo  1  1 


(continuous  measurements  in  time  from  a  finite  number  of  sensors) 
have  been  considered.  The  methods  we  have  proposed  extend  di¬ 
rectly  to  such  functionals.  In  fact,  ours  may  be  considered  a 
digitized  version  of  (5.4),  where  only  a  finite  number  of  time 
samples  are  recorded. 

We  have  only  used  state  observations  corresponding  to  dis¬ 
placements  y(ti  ,Xj )  to  simplify  the  discussion.  It  should  be 
clear  that  we  could  also  have  used  data  of  the  form  y  (t.,x.) 

L  X  J 

(velocity  measurements)  or  data  of  the  form  y  (t.,x.)  (from 

x  1  1 

strain  gauges)  as  well,  with  only  slight  modifications  in  the 
arguments  and  in  the  computer  codes. 

The  theory  that  we  have  developed  (following  [11])  is  the 
semidiscrete  approximations  (5.2)  (i.e.,  continuous  in  time). 

This  serves  to  decouple  the  analysis  for  the  spatial  approxima¬ 
tion  from  the  time  discretization  that  is  employed  in  practice. 
This  is  somewhat  necessary  because  we  have  used  variable-step/ 


variable  order  methods  (Gear's  package)  to  integrate  the  approxi¬ 
mating  system  of  ordinary  differential  equations  to  a  specified 
local  error  tolerance.  While  a  full  analysis  of  such  a  fully 
discretized  method  would  be  difficult,  fully  discrete  methods 
based  on  Pade  approximations  to  exp(_o^^  (q)  t)  have  been  analyzed 
for  functional  differential  equations  [39]  and  this  analysis  is 
being  carried  over  to  our  methods  for  partial  differential  equa¬ 
tions  (for  fixed  step  length  At)  . 

The  semigroup  theoretic  approach  (and  in  particular  the 
Trotter-Kato  theorem)  provides  a  simple  approach  to  obtaining 
convergence  results.  Other  methods  have  been  used  and  may  be  more 
powerful  in  obtaining  results  in  a  broader  class  of  equations. 

For  example,  [  8]  have  used  a  weak  formulation  and  Gronwall -type 
estimates  to  obtain  convergence  results  for  estimation  problems 
involving  a  class  of  non-autonomous  equations. 

We  also  note  that  in  [16]  ,  methods  to  approximate  solutions 
to  an  elliptic  equation  are  proposed,  and  an  a  priori  estimate  of 
the  error  |q^"  -  q|  is  derived.  The  Trotter-Kato  approach  does 
not  appear  to  yield  such  estimates  in  any  easy  fashion. 

In  short,  we  have  presented  a  unified  treatment  of  a  class 

of  parameter  estimation  problems  involving  certain  beam  equations. 

Our  treatment  includes  new  methods  of  approximation  (based  on  the 

N  N 

classical  approximations  P  ja^(q)P  ),  proofs  of  convergence  of  the 
approximate  identification  problem,  and  numerical  results. 
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